Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review

Correlation Functions for the Two-Dimensional Ising Model

Robert Hecht

Published Jun 10, 1967• 70 citations
Publisher

Abstract

The energy-density---energy-density and energy-density---spin correlation functions are evaluated for the two-dimensional square Ising model. Results are obtained at zero magnetic field and near the critical point, i.e., for $\ensuremath{\epsilon}=(\frac{4J}{\mathrm{kT}})|(T\ensuremath{-}\frac{{T}_{c}}{{T}_{c}})|\ensuremath{\ll}1.$ The correlation functions are ${〈{E}_{R}{E}_{0}〉\ensuremath{-}{〈E〉}^{2},={(\frac{2J}{\ensuremath{\pi}})}^{2}{\ensuremath{\epsilon}}^{2}[{{K}_{1}}^{2}(\ensuremath{\epsilon}R)\ensuremath{-}{{K}_{0}}^{2}(\ensuremath{\epsilon}R)],}{〈{E}_{R}{\ensuremath{\sigma}}_{0}〉\ensuremath{-}〈E〉〈\ensuremath{\sigma}〉,=(\frac{2J}{\ensuremath{\pi}})\ensuremath{\epsilon}〈\ensuremath{\sigma}〉\ensuremath{\int}{2\ensuremath{\epsilon}R}^{\ensuremath{\infty}}\mathrm{ds} {s}^{\ensuremath{-}2}{e}^{\ensuremath{-}8},}$ where $R$ is the distance between the points in question and ${K}_{0}$ and ${K}_{1}$ are modified Bessel functions of the second kind. Very similar results hold if the coupling constants are different in the two directions. The above results agree in form with the structure predicted by the "scaling-law" approach.

DOI
10.1103/physrev.158.557

Related papers

Correlations along a Line in the Two-Dimensional Ising Model

Leo P. Kadanoff

The $2n$-spin correlation function for the two-dimensional Ising model at $T={T}_{c}$ is evaluated for the special case in which all the spins lie along a straight line, separated by many lattice constants. The resulting $2n$-spin function is simply a quotient of products of two-spin correlations. A hypothesis of reducibility of fluctuations in the critical state is introduced. This hypothesis asserts that the product of any two local fluctuating quantities in the same neighborhood of space may be effectively replaced by a finite sum of local fluctuating quantities in this neighborhood. As a result, the previously found form for the $2n$-spin function may be used to evaluate the correlation function of $n$ energy densities when all $n$ points lie on a line. The $n$-energy correlation function is simply a sum of products of two-energy correlations. The quotient form for the spin correlation plus scaling is shown to immediately imply the logarithmic specific heat.

Physical Review•Dec 10, 1969•51 citations

Spin-spin correlation functions for the two-dimensional Ising model: Exact theory in the scaling region

Tai Tsun Wu, Barry M. McCoy, Craig A. Tracy +1 more

We compute exactly the spin-spin correlation functions $〈{\ensuremath{\sigma}}_{0,0}{\ensuremath{\sigma}}_{M,N}〉$ for the two-dimensional Ising model on a square lattice in zero magnetic field for $T>{T}_{c}$ and $T<{T}_{c}$. We then analyze the correlation functions in the scaling limit $T\ensuremath{\rightarrow}{T}_{c},{M}^{2}+{N}^{2}\ensuremath{\rightarrow}\ensuremath{\infty}$ such that $(T\ensuremath{-}{T}_{c})$ is fixed. In this scaling limit $〈{\ensuremath{\sigma}}_{0,0}{\ensuremath{\sigma}}_{M,N}〉={R}^{\ensuremath{-}\frac{1}{4}}{F}_{\ifmmode\pm\else\textpm\fi{}}(t)+{R}^{\ensuremath{-}\frac{5}{4}}{F}_{1\ifmmode\pm\else\textpm\fi{}}(t)+o({R}^{\ensuremath{-}\frac{5}{4}})$, where $t$ is the scaling variable $\frac{R}{\ensuremath{\xi}}$ and ${F}_{\ifmmode\pm\else\textpm\fi{}}(t)$ and ${F}_{1\ifmmode\pm\else\textpm\fi{}}(t)$ are the scaling functions ($\ensuremath{\xi}$ is the correlation length). We derive exact expressions for these scaling functions, in terms of a Painlev\'e function of the third kind and analyze both the small- and large-$t$ behavior. A table of values for ${F}_{\ifmmode\pm\else\textpm\fi{}}(t)$ (good to ten significant digits) is also given. As an application we computer the coefficients ${C}_{0\ifmmode\pm\else\textpm\fi{}}$ and ${C}_{1\ifmmode\pm\else\textpm\fi{}}$ in the expansion ${k}_{B}T\ensuremath{\chi}(T)={C}_{0\ifmmode\pm\else\textpm\fi{}}{|1\ensuremath{-}\frac{{T}_{c}}{T}|}^{\ensuremath{-}\frac{7}{4}}+{C}_{1\ifmmode\pm\else\textpm\fi{}}{|1\ensuremath{-}\frac{{T}_{c}}{T}|}^{\ensuremath{-}\frac{3}{4}}+O(1)$ of the zero-field susceptibility $\ensuremath{\chi}(T)$ as $T\ensuremath{\rightarrow}{T}_{c}^{\ifmmode\pm\else\textpm\fi{}}$.

Physical review. B, Solid state•Jan 1, 1976•694 citations

Correlation functions, free energies, and magnetizations in the two-dimensional random-field Ising model

S. L. A. de Queiroz, Robin B. Stinchcombe

Transfer-matrix methods are used to calculate spin-spin correlation functions (G), Helmholtz free energies (f) and magnetizations (m) in the two-dimensional random-field Ising model close to the zero-field bulk critical temperature T(c 0), on long strips of width L=3-18 sites, for binary field distributions. Analysis of the probability distributions of G for varying spin-spin distances R shows that describing the decay of their averaged values by effective correlation lengths is a valid procedure only for not very large R. Connections between field and correlation function distributions at high temperatures are established, yielding approximate analytical expressions for the latter, which are used for computation of the corresponding structure factor. It is shown that, for fixed R/L, the fractional widths of correlation-function distributions saturate asymptotically with L-2.2. Considering an added uniform applied field h, a connection between f(h), m(h), the Gibbs free energy g(m) and the distribution function for the uniform magnetization in a zero uniform field, P0(m), is derived and first illustrated for pure systems, and then applied for nonzero random field. From finite-size scaling and crossover arguments, coupled with numerical data, it is found that the width of P0(m) varies against (nonvanishing, but small) random-field intensity H0 as H(-3/7)(0).

arXiv:cond-mat/0105326•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Aug 28, 2001•5 citations

Two-dimensional Ising model near Tc : Approximation for small magnetic field

Barry M. McCoy, Tai Tsun Wu

The study of the equation of state of the two-dimensional Ising model is initiated by use of the recently calculated $n$-point spin-correlation functions. In the scaling region where $T\ensuremath{\rightarrow}{T}_{c}$ and $h=B{|1\ensuremath{-}\frac{{T}_{c}}{T}|}^{\ensuremath{-}\frac{15}{8}}(\frac{H}{k{T}_{c}})$ is of order 1 ($B$ is a lattice-dependent constant) the free energy below ${T}_{c}$ is of the form $f(T, h)\ensuremath{-}f(T, 0)=\mathrm{const}{(1\ensuremath{-}\frac{{T}_{c}}{T})}^{2}$ $\ensuremath{\Sigma}{n=1}^{\ensuremath{\infty}}(\frac{1}{n!}){C}_{n}^{\ensuremath{-}}{h}^{n}$ where the lattice-dependent constant is chosen so that ${C}_{1}^{\ensuremath{-}}=1$. The constant ${C}_{2}^{\ensuremath{-}}$ has been studied previously. In this paper we study ${C}_{3}^{\ensuremath{-}}$ and ${C}_{4}^{\ensuremath{-}}$ and find as a first approximation $\frac{{C}_{3}^{\ensuremath{-}}}{3!}\ensuremath{\sim}\ensuremath{-}\frac{11}{72{\ensuremath{\pi}}^{2}}$ and $\frac{{C}_{4}^{\ensuremath{-}}}{4!}\ensuremath{\sim}\frac{5}{189\ensuremath{\pi}}+\frac{59}{162{\ensuremath{\pi}}^{3}}$. These are in close agreement with the low-temperature-series values of -0.01538 and 0.0195.

Physical review. B, Condensed matter•Nov 1, 1978•16 citations

One-dimensional random Ising ferromagnets and antiferromagnets

M. Ya. Azbel, Michael Rubinstein

The spin-spin correlation function of a one-dimensional random Ising model decays in a manner inversely proportional to the square of the partition function of the part of the system between the spins. At low temperatures and strong coupling there are logarithmic corrections to the scaling power-law dependence of the correlation length on the magnetic field. The correlation of the boundary spins as a function of the system size is evaluated analytically and numerically. At low temperatures and strong coupling this function has a crossover between different behaviors and its Fourier transform oscillates with wave vector.

Physical review. B, Condensed matter•Oct 1, 1983•5 citations

A Theory of the Second Order Phase Transitions in Spin Systems. III

Masuo Suzuki

Expressions for general correlation functions in Ising ferromagnets are derived in terms of the Lee-Yang theorem. By the use of these expressions, the asymptotic forms of correlation functions in the presence of magnetic field are analyzed, assuming those in the absence of the external field. In particular, the critical indices of spin-spin, energy-density-energy-density, energy-density-spin correlation functions are related with the critical indices of thermal susceptiblity χT0, specific heat, susceptibility χH0, and magnetization. These results agree in form with those predicted by the scaling law approach.

Progress of Theoretical Physics•Feb 1, 1968•25 citations
The energy-density---energy-density and energy-density---spin correlation functions are evaluated for the two-dimensional square Ising model. Results are obtained at zero magnetic field and near the critical point, i.e., for \ensuremathϵ=(4JkT)∣(T\ensuremath−TcTc)∣\ensuremath≪1.\ensuremath{\epsilon}=(\frac{4J}{\mathrm{kT}})|(T\ensuremath{-}\frac{{T}_{c}}{{T}_{c}})|\ensuremath{\ll}1.\ensuremathϵ=(kT4J​)∣(T\ensuremath−Tc​Tc​​)∣\ensuremath≪1. The correlation functions are 〈ERE0〉\ensuremath−〈E〉2,=(2J\ensuremathπ)2\ensuremathϵ2[K12(\ensuremathϵR)\ensuremath−K02(\ensuremathϵR)],〈ER\ensuremathσ0〉\ensuremath−〈E〉〈\ensuremathσ〉,=(2J\ensuremathπ)\ensuremathϵ〈\ensuremathσ〉\ensuremath∫2\ensuremathϵR\ensuremath∞dss\ensuremath−2e\ensuremath−8,{〈{E}_{R}{E}_{0}〉\ensuremath{-}{〈E〉}^{2},={(\frac{2J}{\ensuremath{\pi}})}^{2}{\ensuremath{\epsilon}}^{2}[{{K}_{1}}^{2}(\ensuremath{\epsilon}R)\ensuremath{-}{{K}_{0}}^{2}(\ensuremath{\epsilon}R)],}{〈{E}_{R}{\ensuremath{\sigma}}_{0}〉\ensuremath{-}〈E〉〈\ensuremath{\sigma}〉,=(\frac{2J}{\ensuremath{\pi}})\ensuremath{\epsilon}〈\ensuremath{\sigma}〉\ensuremath{\int}{2\ensuremath{\epsilon}R}^{\ensuremath{\infty}}\mathrm{ds} {s}^{\ensuremath{-}2}{e}^{\ensuremath{-}8},}〈ER​E0​〉\ensuremath−〈E〉2,=(\ensuremathπ2J​)2\ensuremathϵ2[K1​2(\ensuremathϵR)\ensuremath−K0​2(\ensuremathϵR)],〈ER​\ensuremathσ0​〉\ensuremath−〈E〉〈\ensuremathσ〉,=(\ensuremathπ2J​)\ensuremathϵ〈\ensuremathσ〉\ensuremath∫2\ensuremathϵR\ensuremath∞dss\ensuremath−2e\ensuremath−8, where RRR is the distance between the points in question and K0{K}_{0}K0​ and K1{K}_{1}K1​ are modified Bessel functions of the second kind. Very similar results hold if the coupling constants are different in the two directions. The above results agree in form with the structure predicted by the "scaling-law" approach.
The 2n2n2n-spin correlation function for the two-dimensional Ising model at T=TcT={T}_{c}T=Tc​ is evaluated for the special case in which all the spins lie along a straight line, separated by many lattice constants. The resulting 2n2n2n-spin function is simply a quotient of products of two-spin correlations. A hypothesis of reducibility of fluctuations in the critical state is introduced. This hypothesis asserts that the product of any two local fluctuating quantities in the same neighborhood of space may be effectively replaced by a finite sum of local fluctuating quantities in this neighborhood. As a result, the previously found form for the 2n2n2n-spin function may be used to evaluate the correlation function of nnn energy densities when all nnn points lie on a line. The nnn-energy correlation function is simply a sum of products of two-energy correlations. The quotient form for the spin correlation plus scaling is shown to immediately imply the logarithmic specific heat.
We compute exactly the spin-spin correlation functions $〈{\ensuremath{\sigma}}_{0,0}{\ensuremath{\sigma}}_{M,N}〉$ for the two-dimensional Ising model on a square lattice in zero magnetic field for T>TcT>{T}_{c}T>Tc​ and T<TcT<{T}_{c}T<Tc​. We then analyze the correlation functions in the scaling limit $T\ensuremath{\rightarrow}{T}_{c},{M}^{2}+{N}^{2}\ensuremath{\rightarrow}\ensuremath{\infty}$ such that $(T\ensuremath{-}{T}_{c})$ is fixed. In this scaling limit $〈{\ensuremath{\sigma}}_{0,0}{\ensuremath{\sigma}}_{M,N}〉={R}^{\ensuremath{-}\frac{1}{4}}{F}_{\ifmmode\pm\else\textpm\fi{}}(t)+{R}^{\ensuremath{-}\frac{5}{4}}{F}_{1\ifmmode\pm\else\textpm\fi{}}(t)+o({R}^{\ensuremath{-}\frac{5}{4}})$, where ttt is the scaling variable $\frac{R}{\ensuremath{\xi}}$ and ${F}_{\ifmmode\pm\else\textpm\fi{}}(t)$ and ${F}_{1\ifmmode\pm\else\textpm\fi{}}(t)$ are the scaling functions ($\ensuremath{\xi}$ is the correlation length). We derive exact expressions for these scaling functions, in terms of a Painlev\'e function of the third kind and analyze both the small- and large-ttt behavior. A table of values for ${F}_{\ifmmode\pm\else\textpm\fi{}}(t)$ (good to ten significant digits) is also given. As an application we computer the coefficients ${C}_{0\ifmmode\pm\else\textpm\fi{}}$ and ${C}_{1\ifmmode\pm\else\textpm\fi{}}$ in the expansion ${k}_{B}T\ensuremath{\chi}(T)={C}_{0\ifmmode\pm\else\textpm\fi{}}{|1\ensuremath{-}\frac{{T}_{c}}{T}|}^{\ensuremath{-}\frac{7}{4}}+{C}_{1\ifmmode\pm\else\textpm\fi{}}{|1\ensuremath{-}\frac{{T}_{c}}{T}|}^{\ensuremath{-}\frac{3}{4}}+O(1)$ of the zero-field susceptibility $\ensuremath{\chi}(T)$ as $T\ensuremath{\rightarrow}{T}_{c}^{\ifmmode\pm\else\textpm\fi{}}$.
Transfer-matrix methods are used to calculate spin-spin correlation functions (G), Helmholtz free energies (f) and magnetizations (m) in the two-dimensional random-field Ising model close to the zero-field bulk critical temperature T(c 0), on long strips of width L=3-18 sites, for binary field distributions. Analysis of the probability distributions of G for varying spin-spin distances R shows that describing the decay of their averaged values by effective correlation lengths is a valid procedure only for not very large R. Connections between field and correlation function distributions at high temperatures are established, yielding approximate analytical expressions for the latter, which are used for computation of the corresponding structure factor. It is shown that, for fixed R/L, the fractional widths of correlation-function distributions saturate asymptotically with L-2.2. Considering an added uniform applied field h, a connection between f(h), m(h), the Gibbs free energy g(m) and the distribution function for the uniform magnetization in a zero uniform field, P0(m), is derived and first illustrated for pure systems, and then applied for nonzero random field. From finite-size scaling and crossover arguments, coupled with numerical data, it is found that the width of P0(m) varies against (nonvanishing, but small) random-field intensity H0 as H(-3/7)(0).
The study of the equation of state of the two-dimensional Ising model is initiated by use of the recently calculated nnn-point spin-correlation functions. In the scaling region where $T\ensuremath{\rightarrow}{T}_{c}$ and $h=B{|1\ensuremath{-}\frac{{T}_{c}}{T}|}^{\ensuremath{-}\frac{15}{8}}(\frac{H}{k{T}_{c}})$ is of order 1 (BBB is a lattice-dependent constant) the free energy below Tc{T}_{c}Tc​ is of the form $f(T, h)\ensuremath{-}f(T, 0)=\mathrm{const}{(1\ensuremath{-}\frac{{T}_{c}}{T})}^{2}$ $\ensuremath{\Sigma}{n=1}^{\ensuremath{\infty}}(\frac{1}{n!}){C}_{n}^{\ensuremath{-}}{h}^{n}$ where the lattice-dependent constant is chosen so that ${C}_{1}^{\ensuremath{-}}=1$. The constant ${C}_{2}^{\ensuremath{-}}$ has been studied previously. In this paper we study ${C}_{3}^{\ensuremath{-}}$ and ${C}_{4}^{\ensuremath{-}}$ and find as a first approximation $\frac{{C}_{3}^{\ensuremath{-}}}{3!}\ensuremath{\sim}\ensuremath{-}\frac{11}{72{\ensuremath{\pi}}^{2}}$ and $\frac{{C}_{4}^{\ensuremath{-}}}{4!}\ensuremath{\sim}\frac{5}{189\ensuremath{\pi}}+\frac{59}{162{\ensuremath{\pi}}^{3}}$. These are in close agreement with the low-temperature-series values of -0.01538 and 0.0195.
The spin-spin correlation function of a one-dimensional random Ising model decays in a manner inversely proportional to the square of the partition function of the part of the system between the spins. At low temperatures and strong coupling there are logarithmic corrections to the scaling power-law dependence of the correlation length on the magnetic field. The correlation of the boundary spins as a function of the system size is evaluated analytically and numerically. At low temperatures and strong coupling this function has a crossover between different behaviors and its Fourier transform oscillates with wave vector.
Expressions for general correlation functions in Ising ferromagnets are derived in terms of the Lee-Yang theorem. By the use of these expressions, the asymptotic forms of correlation functions in the presence of magnetic field are analyzed, assuming those in the absence of the external field. In particular, the critical indices of spin-spin, energy-density-energy-density, energy-density-spin correlation functions are related with the critical indices of thermal susceptiblity χT0, specific heat, susceptibility χH0, and magnetization. These results agree in form with those predicted by the scaling law approach.