We formulate the corner transfer matrix equations, and the corner transfer matrix renormalization group built on them, as a set of local equations on general planar graphs. These equations depend only on the local subgraphs of the planar graph and are glued together consistently to describe the whole graph. Solving them yields message tensors and corner matrices, which can be used to construct an MPS environment for any connected subgraph. We test the proposed equations on the ferromagnetic Ising model, on a variety of infinite lattices (including hyperbolic ones) and on finite graphs, and find systematic improvement over belief propagation.
We formulate the corner transfer matrix equations, and the corner transfer matrix renormalization group built on them, as a set of local equations on general planar graphs. These equations depend only on the local subgraphs of the planar graph and are glued together consistently to describe the whole graph. Solving them yields message tensors and corner matrices, which can be used to construct an MPS environment for any connected subgraph. We test the proposed equations on the ferromagnetic Ising model, on a variety of infinite lattices (including hyperbolic ones) and on finite graphs, and find systematic improvement over belief propagation.
In the context of tensor network states, we for the first time reformulate the corner transfer matrix renormalization group (CTMRG) method into a variational bilevel optimization algorithm. The solution of the optimization problem corresponds to the fixed-point environment pursued in the conventional CTMRG method, from which the partition function of a classical statistical model, represented by an infinite tensor network, can be efficiently evaluated. The validity of this variational idea is demonstrated by the high-precision calculation of the residual entropy of the dimer model, and is further verified by investigating several typical phase transitions in classical spin models, where the obtained critical points and critical exponents all agree with the best known results in literature. Its extension to three-dimensional tensor networks or quantum lattice models is straightforward, as also discussed briefly.
We revisit the corner transfer matrix renormalization group (CTMRG) method of Nishino and Okunishi for contracting two-dimensional (2D) tensor networks and demonstrate that its performance can be substantially improved by determining the tensors using an eigenvalue solver as opposed to the power method used in CTMRG. We also generalize the variational uniform matrix product state (VUMPS) ansatz for diagonalizing 1D quantum Hamiltonians to the case of 2D transfer matrices and discuss similarities with the corner methods. These two new algorithms will be crucial to improving the performance of variational infinite projected entangled pair state (PEPS) methods.
We propose a tensor-network-based algorithm to study the classical Ising model on an infinitely large hyperbolic lattice with a regular 3D tesselation of identical dodecahedra. We reformulate the corner transfer matrix renormalization group (CTMRG) algorithm from 2D to 3D to reproduce the known results on the cubic lattice. We subsequently generalize the CTMRG to a hyperbolic lattice with dodecahedral cells, which is an infinite-dimensional lattice. We analyze the spontaneous magnetization, von Neumann entropy, and correlation length to find a continuous noncritical phase transition on the dodecahedral lattice. We estimate the phase-transition temperature and find the magnetic critical exponents β=0.4999 and δ=3.007, which confirm the mean-field universality class, in accord with predictions from Monte Carlo and high-temperature series expansions. The algorithm can be applied to arbitrary multistate spin models.
We introduce "matrix product belief propagation" (MP-BP) as a controlled method for the approximate contraction of two-dimensional tensor networks and graphical models, which generically leads to a quadratic reduction of errors relative to existing methods: at similar computational effort, the number of accurate digits is asymptotically doubled. For infinite systems, MP-BP can be understood as a practical implementation of Baxter's corner-transfer-matrix method; for reflection-symmetric networks, it further coincides with the boundary matrix-product-state and corner-transfer-matrix renormalization group algorithms as conventionally implemented. In the limit of unit matrix-product rank, MP-BP is equivalent to belief propagation, and more generally, can be understood as a generalization of BP in which the "messages" live on surfaces and are assumed to have a matrix-product factorization. We further extend the quadratic improvement to the computation of expectation values, and numerically demonstrate the qualitatively improved convergence for a variety of classical and quantum 2D tensor network problems.
The concept of the corner transfer matrix (CTM) was first discovered by Baxter in 1968, when he derived a set of variational matrix equations, which for finite matrix sizes, could be numerically solved to obtain a sequence of approximations for the statistical mechanical properties of a system of monomers and dimers on a rectangular lattice [1]. It was not until 1978, however, in a seminal paper entitled "Variational Approximations for Square Lattice Models in Statistical Mechanics" [4], that Baxter outlined his CTM variational method, which brought to light the potential power of the former objects to obtain numerics and series expansions for unsolved models and to calculate the order parameter of solved ones. Subsequent numerical work led to the realisation that the method, though general, was not very efficient; and increasing efficiency required making model specific modifications, which restricted the transferability of the resulting algorithm to other models. The CTM variational method was thus not widely adopted, despite holding much promise. More recently, however, Nishino and Okunishi discovered that White's density matrix renormalisation group (DMRG) algorithm [55, 56] could be efficiently extended to study two-dimensional classical lattice models, if the density matrix was approximated by Baxters CTMs. Numerical tests of their algorithm, the corner transfer matrix renormalisation group (CTMRG) method, were met with much success. Notable among these is its implementation within the infinite projected entangled-pair states (iPEPS) algorithm of Orus and Vidal to simulate the ground state of infinite two-dimensional quantum systems. In this thesis we review CTM derived variational methods: Baxter's original CTM iterative method, and developments of the CTM concept within the DMRG algorithm by Nishino and Okunishi, which led to the CTMRG method to calculate critical phenomena of classical and quantum systems. This will begin with elucidating the theoretical formalisms of both methods in two dimensions and their generalisations to three dimensions; followed by review of important application, namely, within the iPEPS algorithm of Orus and Vidal to numerically study infinite planar quantum systems.
The field of tensor networks, kicked off in 1992 by Steve White's invention of the spectacularly successful density matrix renormalization group (DMRG) algorithm, has exploded in popularity in recent years. Tensor networks are poised to play a role in helping us solve some of the greatest open physics problems of our time, such as understanding the nature of high-temperature superconductivity and illuminating a theory of quantum gravity. DMRG and extensions based on a class of variational states known as tensor network states have been indispensable tools in helping us understand both numerically and theoretically the properties of complicated classical and quantum many-body systems. However, practical challenges to these techniques still remain, and algorithmic developments are needed before tensor network algorithms can be applied to more physics problems. In this thesis we present a variety of recent advancements to tensor network algorithms. First we describe a DMRG-like algorithm for noninteracting fermions. Noninteracting fermions, naturally being gapless and therefore having high levels of entanglement, are actually a challenging setting for standard DMRG algorithms, and we believe this new algorithm can help with tensor network calculations in that setting. Next we explain a new algorithm called the variational uniform matrix product state (VUMPS) algorithm that is a DMRG-like algorithm that works directly in the thermodynamic limit, improving upon currently available MPS-based methods for studying infinite 1D and quasi-1D quantum many-body systems. Finally, we describe a variety of improvements to algorithms for contracting 2D tensor networks, a common problem in tensor network algorithms, for example for studying 2D classical statistical mechanics problems and 2D quantum many-body problems with projected entangled pair states (PEPS). One is a new variant of the corner transfer matrix renormalization group (CTMRG) algorithm of Nishino and Okunishi that improves the numerical stability for contracting asymmetric two-dimensional tensor networks compared to the most commonly used method. Another is the application of the VUMPS algorithm to contracting 2D tensor networks. The last is a new alternative to CTMRG, where the tensors are solved for with eigenvalue equations instead of a power method, which we call the fixed point corner method (FPCM). We present results showing the transfer matrix VUMPS algorithm and FPCM significantly improve upon the convergence time of CTMRG. We expect these algorithms will play an important role in expanding the set of 2D classical and 2D quantum many-body problems that can be addressed with tensor networks.