arXiv:1708.06115 · Physical review. E

Original electric-vertex formulation of the symmetric eight-vertex model on the square lattice is fully nonuniversal

Roman Krčmár, Ladislav Šamaj

Published Jan 9, 2018• 3 citations

Abstract

The partition function of the symmetric (zero electric field) eight-vertex model on a square lattice can be formulated either in the original "electric" vertex format or in an equivalent "magnetic" Ising-spin format. In this paper, both electric and magnetic versions of the model are studied numerically by using the corner transfer matrix renormalization-group method which provides reliable data. The emphasis is put on the calculation of four specific critical exponents, related by two scaling relations, and of the central charge. The numerical method is first tested in the magnetic format, the obtained dependencies of critical exponents on the model's parameters agree with Baxter's exact solution, and weak universality is confirmed within the accuracy of the method due to the finite size of the system. In particular, the critical exponents η and δ are constant as required by weak universality. On the other hand, in the electric format, analytic formulas based on the scaling relations are derived for the critical exponents η_{e} and δ_{e} which agree with our numerical data. These exponents depend on the model's parameters which is evidence for the full nonuniversality of the symmetric eight-vertex model in the original electric formulation.

DOI: 10.1103/physreve.97.012108

Related papers

Critical properties of the eight-vertex model in a field

Roman Krčmár, Ladislav Šamaj

The general eight-vertex model on a square lattice is studied numerically by using the Corner Transfer Matrix Renormalization Group method. The method is tested on the symmetric (zero-field) version of the model, the obtained dependence of critical exponents on model's parameters is in agreement with Baxter's exact solution and weak universality is verified with a high accuracy. It was suggested long time ago that the symmetric eight-vertex model is a special exceptional case and in the presence of external fields the eight-vertex model falls into the Ising universality class. We confirm numerically this conjecture in a subspace of vertex weights, except for two specific combinations of vertical and horizontal fields for which the system still exhibits weak universality.

arXiv:1610.08657Europhysics Letters (EPL)Sep 1, 20166 citations

Full nonuniversality of the symmetric 16-vertex model on the square lattice

Eva Pospíšilová, Roman Krčmár, Andrej Gendiar +1 more

We consider the symmetric two-state 16-vertex model on the square lattice whose vertex weights are invariant under any permutation of adjacent edge states. The vertex-weight parameters are restricted to a critical manifold which is self-dual under the gauge transformation. The critical properties of the model are studied numerically with the Corner Transfer Matrix Renormalization Group method. Accuracy of the method is tested on two exactly solvable cases: the Ising model and a specific version of the Baxter eight-vertex model in a zero field that belong to different universality classes. Numerical results show that the two exactly solvable cases are connected by a line of critical points with the polarization as the order parameter. There are numerical indications that critical exponents vary continuously along this line in such a way that the weak universality hypothesis is violated.

arXiv:2002.02757Physical review. EJul 13, 20200 citations

Modified universality hypothesis for the eight-vertex model

J. D. Gunton, Th. Niemeijer

This paper is a comment on recent papers by Kadanoff and Wegner, and Barber and Baxter, concerning the critical exponents of the eight-vertex model. It is shown that the critical exponents of this model as derived by these authors can be obtained from the square Ising lattice critical exponents by a renormalization of the temperature only. A new "universality" hypothesis is proposed which is consistent with these results.

Physical review. B, Solid stateJan 1, 197517 citations

Transfer matrix and Monte Carlo tests of critical exponents in lattice models

J. Kaupužs

The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems of square geometry with periodic boundaries oriented either along <10> or along <11> direction have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction proportional to 1/L^0.25 of a very small magnitude also has been detected in agreement with predictions of our recently developed GFD (grouping of Feynman diagrams) theory. A refined analysis of the recent MC data for 3D Ising, phi^4, and XY lattice models has been performed. It includes an analysis of the partition function zeros of 3D Ising model, an estimation of the correction-to-scaling exponent omega from the Binder cumulant data near criticality, as well as a study of the effective critical exponent eta and the effective amplitudes in the asymptotic expansion of susceptibility at the critical point. In all cases a refined analysis is consistent with our (GFD) asymptotic values of the critical exponents (nu=2/3, omega=1/2, eta=1/8 for 3D Ising model and omega=5/9 for 3D XY model), while the actually accepted "conventional" exponents are, in fact, effective exponents which are valid for approximation of the finite--size scaling behavior of not too large systems.

arXiv:cond-mat/0201221Jan 14, 20020 citations

Generalized codes and their application to Ising models with four-spin interactions including the eight-vertex model

David S. Gaunt

Generalized codes are defined for the simple quadratic lattice and it is shown how they may be used to derive high magnetic field or low temperature expansions for two Ising-type models with four-spin interactions. One of these models, namely Baxter's eight-vertex model, has critical exponents which are known to depend on x, the ratio of the strengths of the four-spin to two-spin interactions. The series so derived are analysed to yield estimates of the critical exponent delta for the magnetization as a function of magnetic field at T c . The results for both models are consistent with a constant value delta =15 independent of x as predicted by the scaling laws.

Journal of Physics A Mathematical Nuclear and GeneralSep 1, 19746 citations

The ice and Potts models in statistical mechanics

Stewart Beaton Kelland

This thesis deals with two two-dimensional lattice models in statistical mechanics. The first of these is the zero-field 20-vertex model on a triangular lattice. This model is the analogue to the ice-type, zero-field six-vertex model on a square lattice. We solve this problem exactly for a restricted set of vertex weights. The other model is the q-component Potts model on a square lattice. This model is equivalent to a staggered ice-type model on a square lattice. Using a variational method we estimate the critical exponent ß for several values of q. In part I we outline the fundamental statistical mechanical ideas needed, including the definition of the partition function and its relation to the thermodynamic quantities. Mention is also made of the existence and uniqueness of the thermodynamic limit for these models. The 20-vertex model is discussed in part II. The transfer matrix approach is adopted and a Bethe ansatz is used for the eigenvector. Unlike the square lattice ice-type model, the triangular lattice 20-vertex model is only soluble by this method for a restricted set of vertex weights. In particular, the ten parameters determining the vertex weights must be expressed in terms of four independent parameters. This results in only one type of non-trivial temperature independent solution. This solution exhibits critical behaviour similar to the square lattice KDP model. The general solution has three states: disordered, antiferroelectric and ferroelectric. Which state the model is in is determined by one of the above four parameters. The transition from the disordered to the antiferroelectric state is of infinite-order (i.e the free energy has an essential singularity), and from the disordered to the ferroelectric state is of finite-order. We also compare our results to the free-fermion solution of the same model and note where the two methods overlap. In part III the Potts model is considered. The model is defined and shown to be equivalent to a staggered ice-type model. A variational principle for the eigenvalues of the transfer matrix of the staggered ice model is determined. Using the zero-temperature solution as a ^uide a simple form for a sequence of approximations to the eigenvector, below T , is developed. By using the variational principle we arrive at a set of matrix equations determining the eigenvector at a given level of approximation. By noting some of the symmetries of the transfer matrix we are able to separate some of the dependences of the matrix elements and consequently simplify the original set of equations. We are also able to express the percolation probability ( the analogue of the spontaneous magnetization for general values of q ) of the Potts model as a simple function of these matrix elements. At each level of approximation there is a set of algebraic equations to solve. We have used an electronic computer to solve these equations for the first five levels of approximation, and for several values of q between zero and four. With our results we estimate the critical exponent 3 for the Potts model. These estimates agree reasonably well with previous work and indicate 3=0.145 ±.017 for q=l.

Bulletin of the Australian Mathematical SocietyOct 1, 19761 citation