Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA model project for reproducible papers: critical temperature for the Ising model on a square lattice
52
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this paper we present a simple, yet typical simulation in statistical physics, consisting of large scale Monte Carlo simulations followed by an involved statistical analysis of the results. The purpose is to provide an example publication to explore tools for writing reproducible papers. The simulation estimates the critical temperature where the Ising model on the square lattice becomes magnetic to be Tc /J = 2.26934(6) using a finite size scaling analysis of the crossing points of Binder cumulants. We provide a virtual machine which can be used to reproduce all figures and results.
rate research
Read More
The partition function of the square lattice Ising model on the rectangle with open boundary conditions in both directions is calculated exactly for arbitrary system size $Ltimes M$ and temperature. We start with the dimer method of Kasteleyn, McCoy & Wu, construct a highly symmetric block transfer matrix and derive a factorization of the involved determinant, effectively decomposing the free energy of the system into two parts, $F(L,M)=F_mathrm{strip}(L,M)+F_mathrm{strip}^mathrm{res}(L,M)$, where the residual part $F_mathrm{strip}^mathrm{res}(L,M)$ contains the nontrivial finite-$L$ contributions for fixed $M$. It is given by the determinant of a $frac{M}{2}times frac{M}{2}$ matrix and can be mapped onto an effective spin model with $M$ Ising spins and long-range interactions. While $F_mathrm{strip}^mathrm{res}(L,M)$ becomes exponentially small for large $L/M$ or off-critical temperatures, it leads to important finite-size effects such as the critical Casimir force near criticality. The relations to the Casimir potential and the Casimir force are discussed.
This paper deals with $tilde{chi}^{(6)}$, the six-particle contribution to the magnetic susceptibility of the square lattice Ising model. We have generated, modulo a prime, series coefficients for $tilde{chi}^{(6)}$. The length of the series is sufficient to produce the corresponding Fuchsian linear differential equation (modulo a prime). We obtain the Fuchsian linear differential equation that annihilates the depleted series $Phi^{(6)}=tilde{chi}^{(6)} - {2 over 3} tilde{chi}^{(4)} + {2 over 45} tilde{chi}^{(2)}$. The factorization of the corresponding differential operator is performed using a method of factorization modulo a prime introduced in a previous paper. The depleted differential operator is shown to have a structure similar to the corresponding operator for $tilde{chi}^{(5)}$. It splits into factors of smaller orders, with the left-most factor of order six being equivalent to the symmetric fifth power of the linear differential operator corresponding to the elliptic integral $E$. The right-most factor has a direct sum structure, and using series calculated modulo several primes, all the factors in the direct sum have been reconstructed in exact arithmetics.
We calculate very long low- and high-temperature series for the susceptibility $chi$ of the square lattice Ising model as well as very long series for the five-particle contribution $chi^{(5)}$ and six-particle contribution $chi^{(6)}$. These calculations have been made possible by the use of highly optimized polynomial time modular algorithms and a total of more than 150000 CPU hours on computer clusters. For $chi^{(5)}$ 10000 terms of the series are calculated {it modulo} a single prime, and have been used to find the linear ODE satisfied by $chi^{(5)}$ {it modulo} a prime. A diff-Pade analysis of 2000 terms series for $chi^{(5)}$ and $chi^{(6)}$ confirms to a very high degree of confidence previous conjectures about the location and strength of the singularities of the $n$-particle components of the susceptibility, up to a small set of ``additional singularities. We find the presence of singularities at $w=1/2$ for the linear ODE of $chi^{(5)}$, and $w^2= 1/8$ for the ODE of $chi^{(6)}$, which are {it not} singularities of the ``physical $chi^{(5)}$ and $chi^{(6)},$ that is to say the series-solutions of the ODEs which are analytic at $w =0$. Furthermore, analysis of the long series for $chi^{(5)}$ (and $chi^{(6)}$) combined with the corresponding long series for the full susceptibility $chi$ yields previously conjectured singularities in some $chi^{(n)}$, $n ge 7$. We also present a mechanism of resummation of the logarithmic singularities of the $chi^{(n)}$ leading to the known power-law critical behaviour occurring in the full $chi$, and perform a power spectrum analysis giving strong arguments in favor of the existence of a natural boundary for the full susceptibility $chi$.
We study the class of non-holonomic power series with integer coefficients that reduce, modulo primes, or powers of primes, to algebraic functions. In particular we try to determine whether the susceptibility of the square-lattice Ising model belongs to this class, and more broadly whether the susceptibility is a solution of a differentially algebraic equation. Initial results on Tuttes non-linear ordinary differential equation (ODE) and other simple quadratic non-linear ODEs suggest that a large set of differentially algebraic power series solutions with integer coefficients might reduce to algebraic functions modulo primes, or powers of primes. Here we give several examples of series with integer coefficients and non-zero radius of convergence that reduce to algebraic functions modulo (almost) every prime (or power of a prime). These examples satisfy differentially algebraic equations with the encoding polynomial occasionally possessing quite high degree (and thus difficult to identify even with long series). Additionally, we have extended both the high- and low-temperature Ising square-lattice susceptibility series to 5043 coefficients. We find that even this long series is insufficient to determine whether it reduces to algebraic functions modulo $3$, $5$, etc. This negative result is in contrast to the comparatively easy confirmation that the corresponding series reduce to algebraic functions modulo powers of $2$.
Based on the results published recently [J. Phys. A: Math. Theor. 50, 065201 (2017)], the universal finite-size contributions to the free energy of the square lattice Ising model on the $Ltimes M$ rectangle, with open boundary conditions in both directions, are calculated exactly in the finite-size scaling limit $L,Mtoinfty$, $Tto T_mathrm{c}$, with fixed temperature scaling variable $xpropto(T/T_mathrm{c}-1)M$ and fixed aspect ratio $rhopropto L/M$. We derive exponentially fast converging series for the related Casimir potential and Casimir force scaling functions. At the critical point $T=T_mathrm{c}$ we confirm predictions from conformal field theory by Cardy & Peschel [Nucl. Phys. B 300, 377 (1988)] and by Kleban & Vassileva [J. Phys. A: Math. Gen. 24, 3407 (1991)]. The presence of corners and the related corner free energy has dramatic impact on the Casimir scaling functions and leads to a logarithmic divergence of the Casimir potential scaling function at criticality.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
