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

Journal of Mathematical Physics

Ising Model Spin Correlations on the Triangular Lattice. II. Fourth-Order Correlations

John Stephenson

Published Jun 1, 1966• 72 citations
Publisher

Abstract

A class of fourth-order spin correlations whose sum is closely related to the specific heat is calculated exactly by Pfaffian perturbation theory. In the particular case when the four spins lie on a lattice axis, it is shown that the reduced fourth-order correlation function has the asymptotic form 〈σ0, 0σ1, 1σk, kσk+1, k+1〉−〈σ0, 0σ1, 1〉 〈σk, kσk+1, k+1〉≏A′(e−2βk/k2),where β−1 is the (positive) mean range of order. The decay of correlations with spin separation k is symmetric about the critical point β = 0. The amplitude A′ is 1/π2 at the critical point and is 1/(2π) at all other temperatures. Correlations on ferro- and antiferromagnetic triangular and quadratic lattices are discussed in some detail.

DOI
10.1063/1.1705003

Related papers

Ising-Model Spin Correlations on the Triangular Lattice. III. Isotropic Antiferromagnetic Lattice

John Stephenson

The asymptotic behavior of the pair correlation ω2(r) = 〈σ0σr〉 between two spins at sites 0 and r on an axis of an isotropic antiferromagnetic triangular lattice is investigated with the aid of the theory of Toeplitz determinants as developed by Wu. The leading terms in the asymptotic expansion are obtained for large spin separation at fixed nonzero temperature. Evidence is presented that the zero-point behavior of the correlation is of the form ω2(r) ∼ ε0r−½ cos ⅔πr, where r = |r| is the spin separation and ε0=212(E0T)2=0.632226080…,E0T being the decay amplitude of the pair correlation at the Curie point (critical point) of an isotropic ferromagnetic triangular lattice. A special class of fourth-order correlations ω4(r) = 〈σ0σδσr σr+δ〉 − 〈σ0σδ〉 〈σrσr+δ〉 between the four spins at sites 0, δ, r, and r + δ on the same lattice axis, where δ is a lattice vector, is reconsidered. The asymptotic form of the correlation for large separation of pairs of spins r = |r| is obtained for all fixed temperatures.

Journal of Mathematical Physics•Feb 1, 1970•166 citations

Ising-Model Spin Correlations on the Triangular Lattice

John Stephenson

A Pfaffian representation, of the partition function of the triangular lattice is used to derive expressions for various two, four, and six spin correlations in terms of Pfaffians. The pair correlations along a diagonal are expressed as a Toeplitz determinant whose limiting form yields the spontaneous magnetization. At the ferromagnetic critical point the correlations decay as 1/r12 with approximately radial symmetry. At the antiferromagnetic zero point the ground state is highly degenerate—it has finite entropy—and on a given sublattice the pair correlations along a row decay as ε/r12, where ε=+ε0 on the sublattice containing the origin spin and ε≃−12ε0 on the other two sublattices. Finally, the perpendicular susceptibility, χ⊥, which depends on a finite number of correlations, is calculated; its ferromagnetic behavior is similar to that of the perpendicular susceptibilities of the quadratic and honeycomb lattices, but for an antiferromagnet χ⊥ diverges as 1/T at low temperatures.

Journal of Mathematical Physics•Aug 1, 1964•283 citations

Zero-temperature properties of quantum spin models on the triangular lattice II: the antiferromagnet

Satoshi Fujiki, Donald D. Betts

The calculation of two- and four-spin correlations of the [Formula: see text] antiferromagnet has been extended to an N = 21 site triangular lattice. By fitting a quadratic in 1/N to the nearest neighbour transverse spin correlations, we have estimated the ground-state energy per bond on the infinite lattice to be E0/3NJ = −0.2716 ± 0.005. The nearest neighbour longitudinal correlation is estimated to be [Formula: see text]. A short-range order parameter, the chirality, is defined and estimated for the infinite lattice. From the N dependence of the three sublattice or helical magnetization, the spin correlation is conjectured to decay algebraically as [Formula: see text].

Canadian Journal of Physics•Jan 1, 1987•25 citations

Theory of Critical-Point Scattering and Correlations. I. The Ising Model

Michael Ellis Fisher, Robert J. Burford

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={T}_{c}$ 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>~{T}_{c}$ 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) ${r}_{1}(T)$ is a slowly-varying decreasing function near ${T}_{c}$, (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 ${T}_{c}$ of 0.03 for $d=2$ and of 0.06 to 0.09 for $d=3$. Explicit formulas are given for ${\ensuremath{\kappa}}_{1}$, ${r}_{1}$, and $\ensuremath{\varphi}$ as functions of $T$. The correlations and the scattering are isotropic near ${T}_{c}$. 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>{T}_{c}$. Correspondingly, $\stackrel{^}{\ensuremath{\chi}}(\mathbf{k},T)$ exhibits a maximum for fixed k, at a temperature above ${T}_{c}$; for $d=2$ the maxima are very well marked, but for $d=3$ they are smaller and occur closer to ${T}_{c}$. The theory is compared favorably with recent neutron-scattering experiments on pure beta-brass.

Physical Review•Apr 10, 1967•753 citations

Some exact results for the anisotropic triangular lattice

R. W. Gibberd

The combinatorial approach to the triangular Ising model is shown to be considerably simplified at a given temperature TD. This enables the partition function and spin-correlation functions to be evaluated without the use of pfaffians, thus providing a simple derivation of a few of the results obtained previously by Stephenson. We confirm the results of Stephenson, that at the temperature TD, the correlation functions have the same structure as the correlations of a one-dimensional Ising lattice. The only new result is a general expression for some of the spin–spin correlation functions when the spins do not lie on the axes of the lattice.

Canadian Journal of Physics•Nov 15, 1969•4 citations

Phase transitions in the two-dimensional ferro- and antiferromagnetic potts models on a triangular lattice

Akai Kurbanovich Murtazaev, Albert Babaevich Babaev

The phase transitions in the two-dimensional ferro- and antiferromagnetic Potts models with q = 3 states of spin on a triangular lattice are studied using cluster algorithms and the classical Monte Carlo method. Systems with linear sizes L = 20–120 are considered. The method of fourth-order Binder cumulants and histogram analysis are used to discover that a second-order phase transition occurs in the ferromagnetic Potts model and a first-order phase transition takes place in the antiferromagnetic Potts model. The static critical indices of heat capacity (α), magnetic susceptibility (γ), magnetization (β), and correlation radius index (ν) are calculated for the ferromagnetic Potts model using the finite-size scaling theory.

Journal of Experimental and Theoretical Physics•Dec 1, 2012•17 citations
A class of fourth-order spin correlations whose sum is closely related to the specific heat is calculated exactly by Pfaffian perturbation theory. In the particular case when the four spins lie on a lattice axis, it is shown that the reduced fourth-order correlation function has the asymptotic form 〈σ0, 0σ1, 1σk, kσk+1, k+1〉−〈σ0, 0σ1, 1〉 〈σk, kσk+1, k+1〉≏A′(e−2βk/k2),where β−1 is the (positive) mean range of order. The decay of correlations with spin separation k is symmetric about the critical point β = 0. The amplitude A′ is 1/π2 at the critical point and is 1/(2π) at all other temperatures. Correlations on ferro- and antiferromagnetic triangular and quadratic lattices are discussed in some detail.
The asymptotic behavior of the pair correlation ω2(r) = 〈σ0σr〉 between two spins at sites 0 and r on an axis of an isotropic antiferromagnetic triangular lattice is investigated with the aid of the theory of Toeplitz determinants as developed by Wu. The leading terms in the asymptotic expansion are obtained for large spin separation at fixed nonzero temperature. Evidence is presented that the zero-point behavior of the correlation is of the form ω2(r) ∼ ε0r−½ cos ⅔πr, where r = |r| is the spin separation and ε0=212(E0T)2=0.632226080…,E0T being the decay amplitude of the pair correlation at the Curie point (critical point) of an isotropic ferromagnetic triangular lattice. A special class of fourth-order correlations ω4(r) = 〈σ0σδσr σr+δ〉 − 〈σ0σδ〉 〈σrσr+δ〉 between the four spins at sites 0, δ, r, and r + δ on the same lattice axis, where δ is a lattice vector, is reconsidered. The asymptotic form of the correlation for large separation of pairs of spins r = |r| is obtained for all fixed temperatures.
A Pfaffian representation, of the partition function of the triangular lattice is used to derive expressions for various two, four, and six spin correlations in terms of Pfaffians. The pair correlations along a diagonal are expressed as a Toeplitz determinant whose limiting form yields the spontaneous magnetization. At the ferromagnetic critical point the correlations decay as 1/r12 with approximately radial symmetry. At the antiferromagnetic zero point the ground state is highly degenerate—it has finite entropy—and on a given sublattice the pair correlations along a row decay as ε/r12, where ε=+ε0 on the sublattice containing the origin spin and ε≃−12ε0 on the other two sublattices. Finally, the perpendicular susceptibility, χ⊥, which depends on a finite number of correlations, is calculated; its ferromagnetic behavior is similar to that of the perpendicular susceptibilities of the quadratic and honeycomb lattices, but for an antiferromagnet χ⊥ diverges as 1/T at low temperatures.
The calculation of two- and four-spin correlations of the [Formula: see text] antiferromagnet has been extended to an N = 21 site triangular lattice. By fitting a quadratic in 1/N to the nearest neighbour transverse spin correlations, we have estimated the ground-state energy per bond on the infinite lattice to be E0/3NJ = −0.2716 ± 0.005. The nearest neighbour longitudinal correlation is estimated to be [Formula: see text]. A short-range order parameter, the chirality, is defined and estimated for the infinite lattice. From the N dependence of the three sublattice or helical magnetization, the spin correlation is conjectured to decay algebraically as [Formula: see text].
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=2d=2d=2, but $\ensuremath{\nu}=0.6430\ifmmode\pm\else\textpm\fi{}0.0025\ensuremath{\simeq}\frac{9}{14}$ for d=3d=3d=3. The asymptotic decay of correlation at T=TcT={T}_{c}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=2d=2d=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> TcT>~{T}_{c}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) aaa 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 qqq nearest-neighbor lattice sites, (iii) r1(T){r}_{1}(T)r1​(T) is a slowly-varying decreasing function near Tc{T}_{c}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{T}_{c}Tc​ of 0.03 for d=2d=2d=2 and of 0.06 to 0.09 for d=3d=3d=3. Explicit formulas are given for ${\ensuremath{\kappa}}_{1}$, r1{r}_{1}r1​, and $\ensuremath{\varphi}$ as functions of TTT. The correlations and the scattering are isotropic near Tc{T}_{c}Tc​. The critical scattering isotherm is curved for low kkk according to ${\stackrel{^}{\ensuremath{\chi}}}^{\ensuremath{-}1}\ensuremath{\sim}{k}^{2\ensuremath{-}\ensuremath{\eta}}$ and it intersects the isotherms for T>TcT>{T}_{c}T>Tc​. Correspondingly, $\stackrel{^}{\ensuremath{\chi}}(\mathbf{k},T)$ exhibits a maximum for fixed k, at a temperature above Tc{T}_{c}Tc​; for d=2d=2d=2 the maxima are very well marked, but for d=3d=3d=3 they are smaller and occur closer to Tc{T}_{c}Tc​. The theory is compared favorably with recent neutron-scattering experiments on pure beta-brass.
The combinatorial approach to the triangular Ising model is shown to be considerably simplified at a given temperature TD. This enables the partition function and spin-correlation functions to be evaluated without the use of pfaffians, thus providing a simple derivation of a few of the results obtained previously by Stephenson. We confirm the results of Stephenson, that at the temperature TD, the correlation functions have the same structure as the correlations of a one-dimensional Ising lattice. The only new result is a general expression for some of the spin–spin correlation functions when the spins do not lie on the axes of the lattice.
The phase transitions in the two-dimensional ferro- and antiferromagnetic Potts models with q = 3 states of spin on a triangular lattice are studied using cluster algorithms and the classical Monte Carlo method. Systems with linear sizes L = 20–120 are considered. The method of fourth-order Binder cumulants and histogram analysis are used to discover that a second-order phase transition occurs in the ferromagnetic Potts model and a first-order phase transition takes place in the antiferromagnetic Potts model. The static critical indices of heat capacity (α), magnetic susceptibility (γ), magnetization (β), and correlation radius index (ν) are calculated for the ferromagnetic Potts model using the finite-size scaling theory.