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

arXiv:1104.0037

Specific heat of the Gaussian random bond Ising model on a square lattice

Olga V. Dimitrova

Published Mar 31, 2011• 0 citations
arXivHTMLPDF

Abstract

The free energy and the specific heat of the two-dimensional Gaussian random bond Ising model on a square lattice are found with high accuracy using graph expansion and analysis of continuation of the high-temperature series. At low temperatures the specific heat follows the power law $C\propto T^{3/2}$, the result confirmed independently by counting many-body states in finite size samples. The interpretation of the free energy in terms of independent droplet excitations gives the density of the two-level states, that follows a novel power law $\rho(\epsilon)\propto \sqrt{\epsilon}$ at low energies.

Related papers

Disorder-Induced Phase in the Random Bond Ising Model on the Square Lattice

K. Ziegler

The random bond Ising model in two dimensions is considered in terms of a random matrix representation. A generalization of this approach enables us to study an N → ∞ limit and the corresponding 1/N expansion. With g as the strength of the random bond fluctuations, the critical temperature T c of the pure system is split into two critical temperatures T ± = T c (1 ± c 0 exp[− c 1 / g ]) with a new phase for T − < T < T + . The specific heat, for instance, grows logarithmically (or double logarithmically in a crossover region near T ± ) with|1− T / T c | up to a bound of the order 1/ g .

Europhysics Letters (EPL)•Mar 1, 1991•16 citations

Ground-state and domain-wall energies in the spin-glass region of the two-dimensional ±J random-bond Ising model

Ronald Fisch, Alexander K. Hartmann

The statistics of the ground-state and domain-wall energies for the two-dimensional random-bond Ising model on square lattices with independent, identically distributed bonds having probability $p$ of ${J}_{ij}=\ensuremath{-}1$ and $(1\ensuremath{-}p)$ of ${J}_{ij}=+1$ are studied. We are able to consider large samples of up to ${320}^{2}$ spins by using sophisticated matching algorithms. We study $L\ifmmode\times\else\texttimes\fi{}L$ systems, but we also consider $L\ifmmode\times\else\texttimes\fi{}M$ samples, for different aspect ratios $R=L∕M$. We find that the scaling behavior of the ground-state energy and its sample-to-sample fluctuations inside the spin-glass region $({p}_{c}\ensuremath{\leqslant}p\ensuremath{\leqslant}1\ensuremath{-}{p}_{c})$ are characterized by simple scaling functions. In particular, the fluctuations exhibit a cusplike singularity at ${p}_{c}$. Inside the spin-glass region, the average domain-wall energy converges to a finite nonzero value as the sample size becomes infinite, holding $R$ fixed. Here, large finite-size effects are visible, which can be explained for all $p$ by a single exponent $\ensuremath{\omega}\ensuremath{\approx}2∕3$, provided higher-order corrections to scaling are included. Finally, we confirm the validity of aspect-ratio scaling for $R\ensuremath{\rightarrow}0$: the distribution of the domain-wall energies converges to a Gaussian for $R\ensuremath{\rightarrow}0$, although the domain walls of neighboring subsystems of size $L\ifmmode\times\else\texttimes\fi{}L$ are not independent.

arXiv:cond-mat/0603729•Physical Review B•May 10, 2007•6 citations

Specific heat of a three dimensional Ising ferromagnet above the Curie temperature. II

M F Sykes, D L Hunter, D S McKenzie +1 more

For pt. I see abstr. A37614 of 1967. High temperature series expansions for the specific heat of the Ising model of a ferromagnet are given for the face-centred cubic, body-centred cubic, and simple cubic lattices. From a numerical study it is concluded that the critical index ( alpha ) is lattice independent and that in three dimensions alpha approximately=1/8. A numerical representation of the specific heat in the range T C <or=T<or= infinity is given in each case, together with estimates of the critical energy, entropy, and free energy. Asymptotic forms for the specific heat and energy are given.

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

Shape-dependent finite-size effect of the critical two-dimensional Ising model on a triangular lattice

Xintian Wu, N. Sh. Izmailian, Wenan Guo

Using a bond-propagation algorithm, we study the finite-size behavior of the critical two-dimensional Ising model on a finite triangular lattice with free boundaries in five shapes: triangular, rhomboid, trapezoid, hexagonal, and rectangular. The critical free energy, internal energy, and specific heat are calculated. The accuracy of the free energy reaches 10(-26). Based on accurate data on several finite systems with linear size up to N=2000, we extract the bulk, surface, and corner parts of the free energy, internal energy, and specific heat accurately. We confirm the conformal field theory prediction that the corner free energy is universal and find logarithmic corrections in higher-order terms in the critical free energy for the rhomboid, trapezoid, and hexagonal systems, which are absent for the triangular and rectangular systems. The logarithmic edge corrections due to edges parallel or perpendicular to the bond directions in the internal energy are found to be identical, while the logarithmic edge corrections due to corresponding edges in the specific heat are different. The corner internal energy and corner specific heat for angles π/3, π/2, and 2π/3 are obtained, as well as higher-order corrections. Comparing with the corner internal energy and corner specific heat we previously found on a rectangle of the square lattice [Phys. Rev. E 86, 041149 (2012)], we conclude that the corner internal energy and corner specific heat for the rectangular shape are not universal.

arXiv:1212.2023•Physical Review E•Feb 19, 2013•27 citations

High-temperature series expansions for the (2+1)-dimensional Ising model

Haihua He, Chris J. Hamer, J. Oitmaa

Using efficient cluster expansion methods, the known high-temperature series for the vacuum energy, specific heat, susceptibility and mass gap of the (2+1)-dimensional Ising model on the square and triangular lattices have been extended by several terms. Estimates of the critical indices demonstrate very convincing universality with the Euclidean version of the model.

Journal of Physics A Mathematical and General•May 21, 1990•110 citations

Finite Gluing of the Two-Dimensional Ising Free Energy

Yoshiki Ueoka

We test the Gluing Method on the square-lattice two-dimensional Ising model by reconstructing global free-energy generators from complementary local relations. The emphasis is kept explicit throughout: the inputs used as \emph{gluing material} are separated from quantities obtained \emph{after gluing} and from data reserved for validation. At zero magnetic field, Kramers--Wannier duality is absorbed into the coordinate $\tau=\log\sinh 2K$, so that the reduced free energy decomposes as $\Phi(\tau)=\tau/2+G(|\tau|)$. The gluing materials are then only a critical logarithmic local expansion and a low-order high-temperature outer expansion; the low-temperature sector is generated by duality rather than fitted independently. Positive flat weights preserve the prescribed critical asymptotics to every finite algebraic order. With one crossover parameter fixed from three scalar anchors, the resulting bulk free energy agrees with the exact Onsager benchmark at a maximum relative error of $1.915\times10^{-5}$ on the tested interval. Internal energy and specific heat are derived from the same glued generator. Zero-field finite-size corrections are organized separately as one-variable scaling functions multiplying $L^{-2},L^{-4},\ldots$. At nonzero field in the thermodynamic limit, the gluing materials are the established high-temperature, critical, and low-temperature local expansions of the universal Ising magnetic scaling free energy $\Phi_{\rm mag}(\eta)$, $\eta=m/|h|^{8/15}$. The low-temperature input contains the $|h|$ cusp associated with spontaneous symmetry breaking. Mapping the glued scaling free energy back to the square lattice through nonlinear scaling fields yields a single generator for the magnetic free-energy increment, magnetization, and susceptibility. The construction reproduces high-temperature linear response, the critical isotherm $M\propto |H|^{1/15}$, and the two one-sided limits $M(H\to0^\pm)=\pm M_0$ without separately fitting magnetization branches. Finally, finite size and finite field are combined in the scaling variables $x_t=Lm$ and $x_h=L^{15/8}h$. The two protected one-dimensional axes, a central mixed-response expansion, a low-temperature spectral-rounding relation, and a massive outer description are glued into a two-variable scaling surface $\mathcal Y_{\rm glue}(x_t,x_h)$. The resulting lattice free energy closes in the compact form $\Phi_L^{glue}(K,H)=\Phi_{lat}^{glue}(K,H)+\frac{\mathcal{Y}_{glue}(Lm,L^{15/8}h)- Y_{bulk}(Lm,L^{15/8}h)}{L^2}+O(L^{-4})$ where $m=m(r,H)$ and $h=h(r,H)$ are the nonlinear scaling fields and $Y_{\rm bulk}$ is the thermodynamic singular scaling surface expressed in the same variables. Thus all finite-size/finite-field information is carried by the excess $\mathcal Y_{\rm glue}-Y_{\rm bulk}$, while the thermodynamic lattice free energy remains the previously constructed $\Phi_{\rm lat}^{\rm glue}$. On eleven transfer-matrix points withheld from the construction, the blind free-energy test gives RMS absolute error $3.90\times10^{-4}$ and maximum absolute error $9.55\times10^{-4}$ in $\mathcal Y$. Differentiation of the same free energy produces finite-size magnetization and susceptibility; magnetization is especially stable on the low-temperature side, while the second derivative is more sensitive in crossover regions.

•Sep 2, 2026•0 citations
The free energy and the specific heat of the two-dimensional Gaussian random bond Ising model on a square lattice are found with high accuracy using graph expansion and analysis of continuation of the high-temperature series. At low temperatures the specific heat follows the power law C∝T3/2C\propto T^{3/2}C∝T3/2, the result confirmed independently by counting many-body states in finite size samples. The interpretation of the free energy in terms of independent droplet excitations gives the density of the two-level states, that follows a novel power law ρ(ϵ)∝ϵ\rho(\epsilon)\propto \sqrt{\epsilon}ρ(ϵ)∝ϵ​ at low energies.
The random bond Ising model in two dimensions is considered in terms of a random matrix representation. A generalization of this approach enables us to study an N → ∞ limit and the corresponding 1/N expansion. With g as the strength of the random bond fluctuations, the critical temperature T c of the pure system is split into two critical temperatures T ± = T c (1 ± c 0 exp[− c 1 / g ]) with a new phase for T − < T < T + . The specific heat, for instance, grows logarithmically (or double logarithmically in a crossover region near T ± ) with|1− T / T c | up to a bound of the order 1/ g .
The statistics of the ground-state and domain-wall energies for the two-dimensional random-bond Ising model on square lattices with independent, identically distributed bonds having probability ppp of ${J}_{ij}=\ensuremath{-}1$ and $(1\ensuremath{-}p)$ of Jij=+1{J}_{ij}=+1Jij​=+1 are studied. We are able to consider large samples of up to 3202{320}^{2}3202 spins by using sophisticated matching algorithms. We study $L\ifmmode\times\else\texttimes\fi{}L$ systems, but we also consider $L\ifmmode\times\else\texttimes\fi{}M$ samples, for different aspect ratios R=L∕MR=L∕MR=L∕M. We find that the scaling behavior of the ground-state energy and its sample-to-sample fluctuations inside the spin-glass region $({p}_{c}\ensuremath{\leqslant}p\ensuremath{\leqslant}1\ensuremath{-}{p}_{c})$ are characterized by simple scaling functions. In particular, the fluctuations exhibit a cusplike singularity at pc{p}_{c}pc​. Inside the spin-glass region, the average domain-wall energy converges to a finite nonzero value as the sample size becomes infinite, holding RRR fixed. Here, large finite-size effects are visible, which can be explained for all ppp by a single exponent $\ensuremath{\omega}\ensuremath{\approx}2∕3$, provided higher-order corrections to scaling are included. Finally, we confirm the validity of aspect-ratio scaling for $R\ensuremath{\rightarrow}0$: the distribution of the domain-wall energies converges to a Gaussian for $R\ensuremath{\rightarrow}0$, although the domain walls of neighboring subsystems of size $L\ifmmode\times\else\texttimes\fi{}L$ are not independent.
For pt. I see abstr. A37614 of 1967. High temperature series expansions for the specific heat of the Ising model of a ferromagnet are given for the face-centred cubic, body-centred cubic, and simple cubic lattices. From a numerical study it is concluded that the critical index ( alpha ) is lattice independent and that in three dimensions alpha approximately=1/8. A numerical representation of the specific heat in the range T C <or=T<or= infinity is given in each case, together with estimates of the critical energy, entropy, and free energy. Asymptotic forms for the specific heat and energy are given.
Using a bond-propagation algorithm, we study the finite-size behavior of the critical two-dimensional Ising model on a finite triangular lattice with free boundaries in five shapes: triangular, rhomboid, trapezoid, hexagonal, and rectangular. The critical free energy, internal energy, and specific heat are calculated. The accuracy of the free energy reaches 10(-26). Based on accurate data on several finite systems with linear size up to N=2000, we extract the bulk, surface, and corner parts of the free energy, internal energy, and specific heat accurately. We confirm the conformal field theory prediction that the corner free energy is universal and find logarithmic corrections in higher-order terms in the critical free energy for the rhomboid, trapezoid, and hexagonal systems, which are absent for the triangular and rectangular systems. The logarithmic edge corrections due to edges parallel or perpendicular to the bond directions in the internal energy are found to be identical, while the logarithmic edge corrections due to corresponding edges in the specific heat are different. The corner internal energy and corner specific heat for angles π/3, π/2, and 2π/3 are obtained, as well as higher-order corrections. Comparing with the corner internal energy and corner specific heat we previously found on a rectangle of the square lattice [Phys. Rev. E 86, 041149 (2012)], we conclude that the corner internal energy and corner specific heat for the rectangular shape are not universal.
We test the Gluing Method on the square-lattice two-dimensional Ising model by reconstructing global free-energy generators from complementary local relations. The emphasis is kept explicit throughout: the inputs used as \emph{gluing material} are separated from quantities obtained \emph{after gluing} and from data reserved for validation. At zero magnetic field, Kramers--Wannier duality is absorbed into the coordinate τ=log⁡sinh⁡2K\tau=\log\sinh 2Kτ=logsinh2K, so that the reduced free energy decomposes as Φ(τ)=τ/2+G(∣τ∣)\Phi(\tau)=\tau/2+G(|\tau|)Φ(τ)=τ/2+G(∣τ∣). The gluing materials are then only a critical logarithmic local expansion and a low-order high-temperature outer expansion; the low-temperature sector is generated by duality rather than fitted independently. Positive flat weights preserve the prescribed critical asymptotics to every finite algebraic order. With one crossover parameter fixed from three scalar anchors, the resulting bulk free energy agrees with the exact Onsager benchmark at a maximum relative error of 1.915×10−51.915\times10^{-5}1.915×10−5 on the tested interval. Internal energy and specific heat are derived from the same glued generator. Zero-field finite-size corrections are organized separately as one-variable scaling functions multiplying L−2,L−4,…L^{-2},L^{-4},\ldotsL−2,L−4,…. At nonzero field in the thermodynamic limit, the gluing materials are the established high-temperature, critical, and low-temperature local expansions of the universal Ising magnetic scaling free energy Φmag(η)\Phi_{\rm mag}(\eta)Φmag​(η), η=m/∣h∣8/15\eta=m/|h|^{8/15}η=m/∣h∣8/15. The low-temperature input contains the ∣h∣|h|∣h∣ cusp associated with spontaneous symmetry breaking. Mapping the glued scaling free energy back to the square lattice through nonlinear scaling fields yields a single generator for the magnetic free-energy increment, magnetization, and susceptibility. The construction reproduces high-temperature linear response, the critical isotherm M∝∣H∣1/15M\propto |H|^{1/15}M∝∣H∣1/15, and the two one-sided limits M(H→0±)=±M0M(H\to0^\pm)=\pm M_0M(H→0±)=±M0​ without separately fitting magnetization branches. Finally, finite size and finite field are combined in the scaling variables xt=Lmx_t=Lmxt​=Lm and xh=L15/8hx_h=L^{15/8}hxh​=L15/8h. The two protected one-dimensional axes, a central mixed-response expansion, a low-temperature spectral-rounding relation, and a massive outer description are glued into a two-variable scaling surface Yglue(xt,xh)\mathcal Y_{\rm glue}(x_t,x_h)Yglue​(xt​,xh​). The resulting lattice free energy closes in the compact form ΦLglue(K,H)=Φlatglue(K,H)+Yglue(Lm,L15/8h)−Ybulk(Lm,L15/8h)L2+O(L−4)\Phi_L^{glue}(K,H)=\Phi_{lat}^{glue}(K,H)+\frac{\mathcal{Y}_{glue}(Lm,L^{15/8}h)- Y_{bulk}(Lm,L^{15/8}h)}{L^2}+O(L^{-4})ΦLglue​(K,H)=Φlatglue​(K,H)+L2Yglue​(Lm,L15/8h)−Ybulk​(Lm,L15/8h)​+O(L−4) where m=m(r,H)m=m(r,H)m=m(r,H) and h=h(r,H)h=h(r,H)h=h(r,H) are the nonlinear scaling fields and YbulkY_{\rm bulk}Ybulk​ is the thermodynamic singular scaling surface expressed in the same variables. Thus all finite-size/finite-field information is carried by the excess Yglue−Ybulk\mathcal Y_{\rm glue}-Y_{\rm bulk}Yglue​−Ybulk​, while the thermodynamic lattice free energy remains the previously constructed Φlatglue\Phi_{\rm lat}^{\rm glue}Φlatglue​. On eleven transfer-matrix points withheld from the construction, the blind free-energy test gives RMS absolute error 3.90×10−43.90\times10^{-4}3.90×10−4 and maximum absolute error 9.55×10−49.55\times10^{-4}9.55×10−4 in Y\mathcal YY. Differentiation of the same free energy produces finite-size magnetization and susceptibility; magnetization is especially stable on the low-temperature side, while the second derivative is more sensitive in crossover regions.