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. B, Solid state

General- R High-Temperature Series for the Susceptibility, Second Moment, and Specific Heat of sc and fcc Ising Models with Lattice Anisotropy

Fredric Harbus, Eugene Stanley

Published Jan 1, 1973• 33 citations
Publisher

Abstract

Of particular current interest is the critical behavior of functions on crossing over from one lattice dimensionality to another. To this end, we report high-temperature series for an Ising model with lattice anisotropy-i.e., with different exchange constants for different lattice directions. The Hamiltonian is ${\mathcal{H}}_{l,\mathrm{anis}}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j}\ensuremath{\equiv}\ensuremath{-}{J}_{\mathrm{xy}}(\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}+R\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j})$ where ${s}_{i}=\ifmmode\pm\else\textpm\fi{}1$, the first summation is over all nearest-neighbor pairs in the $\mathrm{xy}$ plane, and the second sum is over pairs coupled in the $z$ direction. The susceptibility, second moment, and specific-heat series are explicitly presented for arbitrary ${J}_{\mathrm{xy}}$ and ${J}_{z}$ for the simple cubic (sc) and face-centered cubic (fcc) lattices to tenth order in inverse temperature. The general-$R$ series are essential if one wishes to study the Riedel-Wegner crossover exponent appropriate to changing lattice dimensionality, since for $R\ensuremath{\equiv}\frac{{J}_{z}}{{J}_{\mathrm{xy}}}=0$, both the sc and fcc lattices reduce to two-dimensional square lattices, while in the limit $R\ensuremath{\rightarrow}\ensuremath{\infty}$, the sc reduces to noninteracting linear chains.

DOI
10.1103/physrevb.7.365

Related papers

Partial Test of the Universality Hypothesis: The Case of Different Coupling Strengths in Different Lattice Directions

Gerald Paul, Eugene Stanley

High-temperature series expansions are used to examine the dependence of critical-point exponents upon lattice anisotropy (different interaction strengths in different directions in the lattice). The two-spin correlation function ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ is calculated to tenth order in ($\frac{1}{{k}_{B}T}$) for the Ising Hamiltonian ${\mathcal{H}}_{l \mathrm{anis}}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}\stackrel{\mathrm{xy}}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}\stackrel{z}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}$ for a wide range of anisotropy parameters $R\ensuremath{\equiv}\frac{{J}_{z}}{{J}_{\mathrm{xy}}}$ and for both the sc and fcc lattices; here the first summation is over all pairs of nearest-neighbor sites whose relative displacement vector ${\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{\mathrm{ij}}$ has no $z$ component, while the second summation is over all other pairs of nearest-neighbor sites. Hence for $R=0$, both the sc and fcc lattices reduce to the two-dimensional square lattice, while in the limit $R\ensuremath{\rightarrow}\ensuremath{\infty}$, the sc becomes a one-dimensional linear chain and the fcc becomes two noninteracting three-dimensional bcc lattices. The series for ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ are then used to obtain series of corresponding lengths for the specific heat, susceptibility, and second moment. Analysis of these series yields results consistent with the universality hypothesis of critical-point exponents. Specifically, it is found that when lattice anisotrophy is introduced, the critical-point exponents studied (the susceptibility exponent $\ensuremath{\gamma}$ and the correlation length exponent $\ensuremath{\nu}$) do not appear to change from their values for an isotropic lattice. The problem of next-nearest-neighbor interactions is treated using similar methods in Paper II of this series (and briefly discussed in this paper).

Physical review. B, Solid state•Apr 1, 1972•56 citations

Scaling with Respect to a Parameter for the Gibbs Potential and Pair Correlation Function of the S=12 Ising Model with Lattice Anisotropy

Richard Krasnow, Fredric Harbus, Luke Lokia Liu +1 more

Series for the reduced susceptibility $\overline{\ensuremath{\chi}}$, the reduced specific heat ${\overline{C}}_{H}$, and second moment ${\ensuremath{\mu}}_{2}$ of the static correlation function for the three-dimensional $S=\frac{1}{2}$ Ising model on both the simple cubic (sc) and face-centered cubic (fcc) lattices with different coupling strengths in different lattice directions have been analyzed to determine the crossover exponent $\ensuremath{\phi}$ describing the behavior of the critical temperature as a function of the anisotropy parameter $R$ in the Hamiltonian $\mathcal{H}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j}\ensuremath{\equiv}\ensuremath{-}{J}_{\mathrm{xy}}(\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}+R\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j})$. Here ${s}_{i}=+1$, the first sum is over all nearest-neighbor pairs in the $\mathrm{xy}$ plane, and the second sum is over all pairs coupled in the $z$ direction. The constant gap exponent we obtain for successive derivatives of $\overline{\ensuremath{\chi}}$ and ${\overline{C}}_{H}$ with respect to $R$ confirms the exponent predictions of scaling in the parameter $R$ for thermodynamic functions, while the results of the ${\ensuremath{\mu}}_{2}$ series confirm the exponent predictions of scaling with respect to $R$ for the two-spin correlation function. Our results agree with the predictions for $\ensuremath{\phi}$ of Abe and Suzuki, and also with rigorous relations satisfied by the exponents describing the derivatives of the various functions. Our results do not agree with previously published results, which are based on an analysis of only the susceptibility on only the sc lattice.

Physical review. B, Solid state•Jan 1, 1973•39 citations

High temperature series for the susceptibility of the Ising model. II. Three dimensional lattices

M F Sykes, David S. Gaunt, P D Roberts +1 more

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.

Journal of physics. A, Proceedings of the Physical Society. General•May 1, 1972•139 citations

Susceptibility scaling functions for ferromagnetic Ising films

T. Weston Capehart, Michael Ellis Fisher

The susceptibility of simple-cubic-lattice ferromagnetic Ising films of $n$ two-dimensional lattice layers is studied by extrapolation of high-temperature series expansions (to eleventh and twelfth orders) for $n=3,4,\dots{},14$ layers with periodic boundary conditions, and $n=3,4,\dots{},10$ layers with free-surface conditions. The corresponding surface susceptibility series for sc, bcc, and fcc lattices are analyzed in the light of universality hypotheses. On the basis of finite-size scaling theory explicit scaling functions are constructed describing the crossover from two-dimensional to three-dimensional critical behavior in terms of the ratio $\frac{n}{\ensuremath{\xi}(T)}$, where $\ensuremath{\xi}\ensuremath{\sim}\ensuremath{\Delta}{T}^{\ensuremath{-}\ensuremath{\nu}}$ is the bulk correlation length.

Physical review. B, Solid state•Jun 1, 1976•123 citations

Extension to order β23 of the high-temperature expansions for the spin- 12 Ising model on simple cubic and body-centered cubic lattices

P. Butera, M. Comi

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).$

arXiv:hep-lat/0006009•Physical review. B, Condensed matter•Dec 1, 2000•40 citations

High-temperature series for the susceptibility of the spin- s Ising model: Analysis of confluent singularities

William J. Camp, John P. VAN DYKE

We have extended the series for the zero-field susceptibility of the spin-$s$ Ising model through tenth order in the reduced inverse temperature $K$ on the square, triangular, simple-cubic, body-centered-cubic, and face-centered-cubic lattices. The series coefficients ${h}_{n}(s)$ are expressed as simple polynomials in $X=s(s+1)$. Using extended methods of analysis we have estimated the nature of the leading singularities on the face-centered lattice and conclude with good confidence that the susceptibility exponent ${\ensuremath{\gamma}}_{1}$ equals 1.25, independent of $s$. The exponent of the leading correction term is estimated to be ${\ensuremath{\gamma}}_{2}\ensuremath{\simeq}0.75\ifmmode\pm\else\textpm\fi{}0.08$ in good agreement with renormalization-group theory. For $s=\frac{1}{2}$ only, the amplitude of the confluent correction apparently vanishes. We have also studied the leading singularities on the triangular lattice and conclude that ${\ensuremath{\gamma}}_{1}=\frac{7}{4}$, independent of $s$, with, however, much stronger corrections than in three dimensions. These results provide a very strong corroboration of the universality hypothesis.

Physical review. B, Solid state•Apr 1, 1975•148 citations
Of particular current interest is the critical behavior of functions on crossing over from one lattice dimensionality to another. To this end, we report high-temperature series for an Ising model with lattice anisotropy-i.e., with different exchange constants for different lattice directions. The Hamiltonian is Hl,anis=\ensuremath−Jxy\ensuremathΣ〈ij〉xysisj\ensuremath−Jz\ensuremathΣ〈ij〉zsisj\ensuremath≡\ensuremath−Jxy(\ensuremathΣ〈ij〉xysisj+R\ensuremathΣ〈ij〉zsisj){\mathcal{H}}_{l,\mathrm{anis}}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j}\ensuremath{\equiv}\ensuremath{-}{J}_{\mathrm{xy}}(\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}+R\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j})Hl,anis​=\ensuremath−Jxy​\ensuremathΣ〈ij〉xysi​sj​\ensuremath−Jz​\ensuremathΣ〈ij〉zsi​sj​\ensuremath≡\ensuremath−Jxy​(\ensuremathΣ〈ij〉xysi​sj​+R\ensuremathΣ〈ij〉zsi​sj​) where si=\ifmmode±\else\textpm\fi1{s}_{i}=\ifmmode\pm\else\textpm\fi{}1si​=\ifmmode±\else\textpm\fi1, the first summation is over all nearest-neighbor pairs in the xy\mathrm{xy}xy plane, and the second sum is over pairs coupled in the zzz direction. The susceptibility, second moment, and specific-heat series are explicitly presented for arbitrary Jxy{J}_{\mathrm{xy}}Jxy​ and Jz{J}_{z}Jz​ for the simple cubic (sc) and face-centered cubic (fcc) lattices to tenth order in inverse temperature. The general-RRR series are essential if one wishes to study the Riedel-Wegner crossover exponent appropriate to changing lattice dimensionality, since for R\ensuremath≡JzJxy=0R\ensuremath{\equiv}\frac{{J}_{z}}{{J}_{\mathrm{xy}}}=0R\ensuremath≡Jxy​Jz​​=0, both the sc and fcc lattices reduce to two-dimensional square lattices, while in the limit R\ensuremath→\ensuremath∞R\ensuremath{\rightarrow}\ensuremath{\infty}R\ensuremath→\ensuremath∞, the sc reduces to noninteracting linear chains.
High-temperature series expansions are used to examine the dependence of critical-point exponents upon lattice anisotropy (different interaction strengths in different directions in the lattice). The two-spin correlation function ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ is calculated to tenth order in (1kBT\frac{1}{{k}_{B}T}kB​T1​) for the Ising Hamiltonian ${\mathcal{H}}_{l \mathrm{anis}}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}\stackrel{\mathrm{xy}}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}\stackrel{z}{〈\mathrm{ij}〉}{S}_{i}^{z}{S}_{j}^{z}$ for a wide range of anisotropy parameters $R\ensuremath{\equiv}\frac{{J}_{z}}{{J}_{\mathrm{xy}}}$ and for both the sc and fcc lattices; here the first summation is over all pairs of nearest-neighbor sites whose relative displacement vector ${\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{\mathrm{ij}}$ has no zzz component, while the second summation is over all other pairs of nearest-neighbor sites. Hence for R=0R=0R=0, both the sc and fcc lattices reduce to the two-dimensional square lattice, while in the limit $R\ensuremath{\rightarrow}\ensuremath{\infty}$, the sc becomes a one-dimensional linear chain and the fcc becomes two noninteracting three-dimensional bcc lattices. The series for ${C}_{2}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$ are then used to obtain series of corresponding lengths for the specific heat, susceptibility, and second moment. Analysis of these series yields results consistent with the universality hypothesis of critical-point exponents. Specifically, it is found that when lattice anisotrophy is introduced, the critical-point exponents studied (the susceptibility exponent $\ensuremath{\gamma}$ and the correlation length exponent $\ensuremath{\nu}$) do not appear to change from their values for an isotropic lattice. The problem of next-nearest-neighbor interactions is treated using similar methods in Paper II of this series (and briefly discussed in this paper).
Series for the reduced susceptibility $\overline{\ensuremath{\chi}}$, the reduced specific heat C‾H{\overline{C}}_{H}CH​, and second moment ${\ensuremath{\mu}}_{2}$ of the static correlation function for the three-dimensional S=12S=\frac{1}{2}S=21​ Ising model on both the simple cubic (sc) and face-centered cubic (fcc) lattices with different coupling strengths in different lattice directions have been analyzed to determine the crossover exponent $\ensuremath{\phi}$ describing the behavior of the critical temperature as a function of the anisotropy parameter RRR in the Hamiltonian $\mathcal{H}=\ensuremath{-}{J}_{\mathrm{xy}}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}\ensuremath{-}{J}_{z}\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j}\ensuremath{\equiv}\ensuremath{-}{J}_{\mathrm{xy}}(\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{\mathrm{xy}}{s}_{i}{s}_{j}+R\ensuremath{\Sigma}{〈\mathrm{ij}〉}^{z}{s}_{i}{s}_{j})$. Here si=+1{s}_{i}=+1si​=+1, the first sum is over all nearest-neighbor pairs in the xy\mathrm{xy}xy plane, and the second sum is over all pairs coupled in the zzz direction. The constant gap exponent we obtain for successive derivatives of $\overline{\ensuremath{\chi}}$ and C‾H{\overline{C}}_{H}CH​ with respect to RRR confirms the exponent predictions of scaling in the parameter RRR for thermodynamic functions, while the results of the ${\ensuremath{\mu}}_{2}$ series confirm the exponent predictions of scaling with respect to RRR for the two-spin correlation function. Our results agree with the predictions for $\ensuremath{\phi}$ of Abe and Suzuki, and also with rigorous relations satisfied by the exponents describing the derivatives of the various functions. Our results do not agree with previously published results, which are based on an analysis of only the susceptibility on only the sc lattice.
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.
The susceptibility of simple-cubic-lattice ferromagnetic Ising films of nnn two-dimensional lattice layers is studied by extrapolation of high-temperature series expansions (to eleventh and twelfth orders) for n=3,4,…,14n=3,4,\dots{},14n=3,4,…,14 layers with periodic boundary conditions, and n=3,4,…,10n=3,4,\dots{},10n=3,4,…,10 layers with free-surface conditions. The corresponding surface susceptibility series for sc, bcc, and fcc lattices are analyzed in the light of universality hypotheses. On the basis of finite-size scaling theory explicit scaling functions are constructed describing the crossover from two-dimensional to three-dimensional critical behavior in terms of the ratio $\frac{n}{\ensuremath{\xi}(T)}$, where $\ensuremath{\xi}\ensuremath{\sim}\ensuremath{\Delta}{T}^{\ensuremath{-}\ensuremath{\nu}}$ is the bulk correlation length.
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).$
We have extended the series for the zero-field susceptibility of the spin-sss Ising model through tenth order in the reduced inverse temperature KKK on the square, triangular, simple-cubic, body-centered-cubic, and face-centered-cubic lattices. The series coefficients hn(s){h}_{n}(s)hn​(s) are expressed as simple polynomials in X=s(s+1)X=s(s+1)X=s(s+1). Using extended methods of analysis we have estimated the nature of the leading singularities on the face-centered lattice and conclude with good confidence that the susceptibility exponent ${\ensuremath{\gamma}}_{1}$ equals 1.25, independent of sss. The exponent of the leading correction term is estimated to be ${\ensuremath{\gamma}}_{2}\ensuremath{\simeq}0.75\ifmmode\pm\else\textpm\fi{}0.08$ in good agreement with renormalization-group theory. For s=12s=\frac{1}{2}s=21​ only, the amplitude of the confluent correction apparently vanishes. We have also studied the leading singularities on the triangular lattice and conclude that ${\ensuremath{\gamma}}_{1}=\frac{7}{4}$, independent of sss, with, however, much stronger corrections than in three dimensions. These results provide a very strong corroboration of the universality hypothesis.