Local and global magnetization on the Sierpiński carpet.

The phase transition of the classical Ising model on the Sierpiński carpet, which has the fractal dimension log_{3}^{}8≈1.8927, is studied by an adapted variant of the higher-order tensor renormalization group method. The second-order phase transition is observed at the critical temperature T_{c}^{}≈1.

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.

J_{1}-J_{2} fractal studied by multirecursion tensor-network method.

We generalize a tensor-network algorithm to study the thermodynamic properties of self-similar spin lattices constructed on a square-lattice frame with two types of couplings, J_{1}^{} and J_{2}^{}, chosen to transform a regular square lattice (J_{1}^{}=J_{2}^{}) onto a fractal lattice if decreasing J_{2}^{} to zero (the fractal fully reconstructs when J_{2}^{}=0). We modified the higher-order tensor renormalization group (HOTRG) algorithm for this purpose. Single-site measurements are performed by means of so-called impurity tensors.

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.

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.

Boundary-field-driven control of discontinuous phase transitions on hyperbolic lattices.

The multistate Potts models on two-dimensional hyperbolic lattices are studied with respect to various boundary effects. The free energy is numerically calculated using the corner transfer matrix renormalization group method. We analyze phase transitions of the Potts models in the thermodynamic limit with respect to contracted boundary layers.

Free-energy analysis of spin models on hyperbolic lattice geometries.

We investigate relations between spatial properties of the free energy and the radius of Gaussian curvature of the underlying curved lattice geometries. For this purpose we derive recurrence relations for the analysis of the free energy normalized per lattice site of various multistate spin models in the thermal equilibrium on distinct non-Euclidean surface lattices of the infinite sizes. Whereas the free energy is calculated numerically by means of the corner transfer matrix renormalization group algorithm, the radius of curvature has an analytic expression.

Phase transition of the Ising model on a fractal lattice.

The phase transition of the Ising model is investigated on a planar lattice that has a fractal structure. On the lattice, the number of bonds that cross the border of a finite area is doubled when the linear size of the area is extended by a factor of 4. The free energy and the spontaneous magnetization of the system are obtained by means of the higher-order tensor renormalization group method.

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.