J1-J2 fractal studied by multi-recursion tensor-network method

We generalize a tensor-network algorithm to study 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. So far, only a single local tensor and uniform extension-contraction relations have been considered in HOTRG. We introduce ten independent local tensors, each being extended and contracted by fifteen different recursion relations. We applied the Ising model to the $J_{1}^{}-J_{2}^{}$ planar fractal whose Hausdorff dimension at $J_{2}^{} = 0$ is $d^{(H)} = ln 12 / ln 4 approx 1.792$. The generalized tensor-network algorithm is applicable to a wide range of fractal patterns and is suitable for models without translational invariance.

Phase transition of the four-dimensional cross-polytope model

Thermodynamic properties of the four-dimensional cross-polytope model, the 16-cell model, which is an example of higher dimensional generalizations of the octahedron model, are studied on the square lattice. By means of the corner transfer matrix renormalization group (CTMRG) method, presence of the first-order phase transition is confirmed. The latent heat is estimated to be $L_4^{~} = 0.3172$, which is larger than that of the octahedron model $L_3^{~} = 0.0516$. The result suggests that the latent heat increases with the internal dimension $n$ when the higher-dimensional series of the cross-polytope models is considered.

Study of classical and quantum phase transitions on non-Euclidean geometries in higher dimensions

The investigation of the behavior of both classical and quantum systems on non-Euclidean surfaces near the phase transition point represents an interesting research area of modern physics. In the case of classical spin systems, a generalization of the Corner Transfer Matrix Renormalization Group algorithm has been developed and successfully applied to spin models on infinitely many regular hyperbolic lattices. In this work, we extend these studies to specific types of lattices. It is important to say that no suitable algorithms for numerical analysis of ground-states of quantum systems in similar conditions have been implemented yet. In this work, we offer a particular solution by proposing a variational numerical algorithm Tensor Product Variational Formulation, which assumes a quantum ground-state written in the form of a low-dimensional uniform tensor product state. We apply the Tensor Product Variational Formulation to three typical quantum models on a variety of regular hyperbolic lattices. The main outcomes are the following: (1) We propose an algorithm for calculation and classification of the thermodynamic properties of the Ising model on triangular-tiled hyperbolic lattices. In addition, we investigate the origin of the mean-field universality on a series of weakly curved lattices. (2) We develop the Tensor Product Variational Formulation algorithm for the numerical analysis of the ground-state of the quantum systems on the hyperbolic lattices. (3) We study quantum phase transition phenomena for the three selected spin models on various types of the hyperbolic lattices including the Bethe lattice.

Tensor Networks: Phase transition phenomena on hyperbolic and fractal geometries

One of the challenging problems in the condensed matter physics is to understand the quantum many-body systems, especially, their physical mechanisms behind. Since there are only a few complete analytical solutions of these systems, several numerical simulation methods have been proposed in recent years. Amongst all of them, the Tensor Network algorithms have become increasingly popular in recent years, especially for their adaptability to simulate strongly correlated systems. The current work focuses on the generalization of such Tensor-Network-based algorithms, which are sufficiently robust to describe critical phenomena and phase transitions of multistate spin Hamiltonians in the thermodynamic limit. We have chosen two algorithms: the Corner Transfer Matrix Renormalization Group and the Higher-Order Tensor Renormalization Group. This work, based on tensor-network analysis, opens doors for the understanding of phase transition and entanglement of the interacting systems on the non-Euclidean geometries. We focus on three main topics: A new thermodynamic model of social influence, free energy is analyzed to classify the phase transitions on an infinite set of the negatively curved geometries where a relation between the free energy and the Gaussian radius of the curvature is conjectured, a unique tensor-based algorithm is proposed to study the phase transition on fractal structures.

Entanglement-entropy study of phase transitions in six-state clock model

The Berezinskii-Kosterlitz-Thouless (BKT) transitions of the six-state clock model on the square lattice are investigated by means of the corner-transfer matrix renormalization group method. A classical analog of the entanglement entropy $S( L, T )$ is calculated for $L times L$ square system up to $L = 129$, as a function of temperature $T$. The entropy exhibits a peak at $T = T^*_{~}( L )$, where the temperature depends on both $L$ and the boundary conditions. Applying the finite-size scaling to $T^*_{~}( L )$ and assuming presence of the BKT transitions, the two distinct phase-transition temperatures are estimated to be $T_1^{~} = 0.70$ and $T_2^{~} = 0.88$. The results are in agreement with earlier studies. It should be noted that no thermodynamic functions have been used in this study.

Area-Law Study of Quantum Spin System on Hyperbolic Lattice Geometries

Magnetic properties of the transverse-field Ising model on curved (hyperbolic) lattices are studied by a tensor product variational formulation that we have generalized for this purpose. First, we identify the quantum phase transition for each hyperbolic lattice by calculating the magnetization. We study the entanglement entropy at the phase transition in order to analyze the correlations of various subsystems located at the center with the rest of the lattice. We confirm that the entanglement entropy satisfies the area law at the phase transition for fixed coordination number, i.e., it scales linearly with the increasing size of the subsystems. On the other hand, the entanglement entropy decreases as power-law with respect to the increasing coordination number.

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 by using 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 8-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.

Measurements of magnetization on the Sierpi{n}ski carpet

Phase transition of the classical Ising model on the Sierpi{n}ski carpet, which has the fractal dimension $log_3^{~} 8 approx 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_{rm c}^{~} = 1.4783(1)$. Position dependence of local functions is studied by means of impurity tensors, which are inserted at different locations on the fractal lattice. The critical exponent $beta$ associated with the local magnetization varies by two orders of magnitude, depending on lattice locations, whereas $T_{rm c}^{~}$ is not affected.

Phase Transition of the Ising Model on Fractal Lattice

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 four. The free energy and the spontaneous magnetization of the system are obtained by means of the higher-order tensor renormalization group method. The system exhibits the order-disorder phase transition, where the critical indices are different from that of the square-lattice Ising model. An exponential decay is observed in the density matrix spectrum even at the critical point. It is possible to interpret that the system is less entangled because of the fractal geometry.