Ising ferromagnets and antiferromagnets in an imaginary magnetic field.

We study classical Ising spin-1/2 models on a two-dimensional (2D) square lattice with ferromagnetic or antiferromagnetic nearest-neighbor interactions, under the effect of a pure imaginary magnetic field. The complex Boltzmann weights of spin configurations cannot be interpreted as a probability distribution, which prevents application of standard statistical algorithms. In this work, the mapping of the Ising spin models under consideration onto symmetric vertex models leads to real (positive or negative) Boltzmann weights.

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

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.

Finite-m scaling analysis of Berezinskii-Kosterlitz-Thouless phase transitions and entanglement spectrum for the six-state clock model.

We investigate the Berezinskii-Kosterlitz-Thouless transitions for the square-lattice six-state clock model with the corner-transfer matrix renormalization group (CTMRG). Scaling analyses for effective correlation length, magnetization, and entanglement entropy with respect to the cutoff dimension m at the fixed point of the CTMRG provide transition temperatures consistent with a variety of recent numerical studies. We also reveal that the fixed-point spectrum of the corner-transfer matrix in the critical intermediate phase of the six-state clock model is characterized by the scaling dimension consistent with the c=1 boundary conformal field theory associated with the effective Z_{6} dual sine-Gordon model.

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

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.

Critical behavior of the two-dimensional icosahedron model.

In the context of a discrete analog of the classical Heisenberg model, we investigate the critical behavior of the icosahedron model, where the interaction energy is defined as the inner product of neighboring vector spins of unit length pointing to the vertices of the icosahedron. The effective correlation length and magnetization of the model are calculated by means of the corner-transfer-matrix renormalization group (CTMRG) method. A scaling analysis with respect to the cutoff dimension m in CTMRG reveals a second-order phase transition characterized by the exponents ν=1.

Phase diagram of a truncated tetrahedral model.

Phase diagram of a discrete counterpart of the classical Heisenberg model, the truncated tetrahedral model, is analyzed on the square lattice, when the interaction is ferromagnetic. Each spin is represented by a unit vector that can point to one of the 12 vertices of the truncated tetrahedron, which is a continuous interpolation between the tetrahedron and the octahedron. Phase diagram of the model is determined by means of the statistical analog of the entanglement entropy, which is numerically calculated by the corner transfer matrix renormalization group method.

Reentrant disorder-disorder transitions in generalized multicomponent Widom-Rowlinson models.

In the lattice version of the multicomponent Widom-Rowlinson (WR) model, each site can be either empty or singly occupied by one of M different particles, all species having the same fugacity z. The only nonzero interaction potential is a nearest-neighbor hard-core exclusion between unlike particles. For Mz(d) (M).

Mean-field universality class induced by weak hyperbolic curvatures.

Order-disorder phase transition of the ferromagnetic Ising model is investigated on a series of two-dimensional lattices that have negative Gaussian curvatures. Exceptional lattice sites of coordination number seven are distributed on the triangular lattice, where the typical distance between the nearest exceptional sites is proportional to an integer parameter n. Thus, the corresponding curvature is asymptotically proportional to -n(-2).

Weak correlation effects in the Ising model on triangular-tiled hyperbolic lattices.

The Ising model is studied on a series of hyperbolic two-dimensional lattices which are formed by tessellation of triangles on negatively curved surfaces. In order to treat the hyperbolic lattices, we propose a generalization of the corner transfer matrix renormalization group method using a recursive construction of asymmetric transfer matrices. Studying the phase transition, the mean-field universality is captured by means of a precise analysis of thermodynamic functions.