Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2401.07274 · Physical review. E

Corner transfer matrix renormalization group approach in the zoo of Archimedean lattices

I. V. Lukin, Andrii Sotnikov

Published Apr 17, 2024• 12 citations
arXivHTMLPDF

Abstract

We develop a new methodology to contract tensor networks within the corner transfer matrix renormalization group approach for a wide range of two-dimensional lattice geometries. We discuss contraction algorithms on the example of triangular, kagome, honeycomb, square-octagon, star, ruby, square-hexagon-dodecahedron, and dice lattices. As benchmark tests, we apply the developed method to the classical Ising model on different lattices and observe a remarkable agreement of the results with the available from the literature. The approach also shows the necessary potential to be applied to various quantum lattice models in a combination with the wave-function variational optimization schemes.

DOI
10.1103/physreve.109.045305

Related papers

Global optimization of tensor renormalization group using the corner transfer matrix

Satoshi Morita, Naoki Kawashima

A tensor network renormalization algorithm with global optimization based on the corner transfer matrix is proposed. Since the environment is updated by the corner transfer matrix renormalization group method, the forward-backward iteration is unnecessary, which is a time-consuming part of other methods with global optimization. In addition, a further approximation reducing the order of the computational cost of contraction for the calculation of the coarse-grained tensor is proposed. The computational time of our algorithm in two dimensions scales as the sixth power of the bond dimension, while the higher-order tensor renormalization group and the higher-order second renormalization group methods have the seventh power. We perform benchmark calculations in the Ising model on the square lattice and show that the time-to-solution of the proposed algorithm is faster than that of other methods.

arXiv:2009.01997•Physical review. B./Physical review. B•Jan 22, 2021•12 citations

Variational Corner Transfer Matrix Renormalization Group Method for Classical Statistical Models

X. F. Liu, Yu Fu, Weiqiang Yu +2 more

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.

arXiv:2203.17098•Chinese Physics Letters•Jun 1, 2022•13 citations

Simulation of two-dimensional quantum systems on an infinite lattice revisited: Corner transfer matrix for tensor contraction

Román Orús, Guifré Vidal

An extension of the projected entangled-pair states (PEPS) algorithm to infinite systems, known as the iPEPS algorithm, was recently proposed to compute the ground state of quantum systems on an infinite two dimensional lattice. Here we investigate a modification of the iPEPS algorithm, where the environment is computed using the corner transfer matrix renormalization group (CTMRG) method, instead of using one-dimensional transfer matrix methods as in the original proposal. We describe a variant of the CTMRG that addresses different directions of the lattice independently, and use it combined with imaginary time evolution to compute the ground state of the two dimensional quantum Ising model. Near criticality, the modified iPEPS algorithm is seen to provide a better estimation of the order parameter and correlators.

arXiv:0905.3225•Physical Review B•Sep 8, 2009•386 citations

Variational optimization of projected entangled-pair states on the triangular lattice

Jan Naumann, Jens Eisert, Philipp Schmoll

We introduce a general corner transfer matrix renormalization group algorithm tailored to projected entangled-pair states on the triangular lattice. By integrating automatic differentiation, our approach enables direct variational energy minimization on this lattice geometry. In contrast to conventional approaches that map the triangular lattice onto a square lattice with diagonal next-nearest-neighbor interactions, our native formulation yields improved variational results at the same bond dimension. This improvement stems from a more faithful and physically informed representation of the entanglement structure in the tensor network and an increased number of variational parameters. We apply our method to the antiferromagnetic nearest-neighbor Heisenberg model on the triangular and kagome lattice, and benchmark our results against previous numerical studies.

arXiv:2510.04907•Physical review. B./Physical review. B•Dec 11, 2025•1 citation

Variational optimization of projected entangled-pair states on the triangular lattice

Jan Naumann, Jens Eisert, Philipp Schmoll

We introduce a general corner transfer matrix renormalization group algorithm tailored to projected entangled-pair states on the triangular lattice. By integrating automatic differentiation, our approach enables direct variational energy minimization on this lattice geometry. In contrast to conventional approaches that map the triangular lattice onto a square lattice with diagonal next-nearest-neighbor interactions, our native formulation yields improved variational results at the same bond dimension. This improvement stems from a more faithful and physically informed representation of the entanglement structure in the tensor network and an increased number of variational parameters. We apply our method to the antiferromagnetic nearest-neighbor Heisenberg model on the triangular and kagome lattice, and benchmark our results against previous numerical studies.

•Jan 1, 2026•0 citations

Development of Two-Dimensional Tensor Network Methods for the Variational Simulation of Quantum Many-Body States

Jan Naumann

Strongly correlated quantum systems are one of the most interesting and challenging phenomena in modern condensed matter physics. Unfortunately, the theoretical access even to their low energy properties is limited by the exponential growth of the Hilbert space in a many-body setting. To mitigate this computational barrier, tensor network methods are one of the most competitive approaches to classically simulate complex quantum systems. In this thesis we will discuss recent developments in the field of tensor networks to efficiently and accurately simulate two-dimensional quantum many-body states in the thermodynamic limit. Hereby, we introduce the variational update globally optimizing the wave-function, encoded as tensor network in the form of an infinite projected entangled pair state (iPEPS), in order to obtain ground states of many-body Hamiltonian models. We focus in the discussion on the well-established method to approximately contract this two-dimensional structure, the corner transfer matrix renormalization group (CTMRG), and discuss its native application on square lattice and triangular lattice based models. Furthermore, we report the development of a refined numerical algorithm, the split-CTMRG, to significantly reduce its computational complexity, pushing the limits of feasible simulations within modern hardware capabilities. Finally, the method advancements are applied to a currently widely studied model of frustrated magnetism, the Heisenberg maple-leaf antiferromagnet. We are able to find a large phase of vanishing magnetization in contradiction to previous works and, to the best of our knowledge, the lowest reported variational ground state energy at the isotropic point of the model. As an outlook, we describe a potential adaption of the previously developed methods to variationally simulate two-dimensional many-body thermal states.

•Jan 1, 2026•0 citations
We develop a new methodology to contract tensor networks within the corner transfer matrix renormalization group approach for a wide range of two-dimensional lattice geometries. We discuss contraction algorithms on the example of triangular, kagome, honeycomb, square-octagon, star, ruby, square-hexagon-dodecahedron, and dice lattices. As benchmark tests, we apply the developed method to the classical Ising model on different lattices and observe a remarkable agreement of the results with the available from the literature. The approach also shows the necessary potential to be applied to various quantum lattice models in a combination with the wave-function variational optimization schemes.
A tensor network renormalization algorithm with global optimization based on the corner transfer matrix is proposed. Since the environment is updated by the corner transfer matrix renormalization group method, the forward-backward iteration is unnecessary, which is a time-consuming part of other methods with global optimization. In addition, a further approximation reducing the order of the computational cost of contraction for the calculation of the coarse-grained tensor is proposed. The computational time of our algorithm in two dimensions scales as the sixth power of the bond dimension, while the higher-order tensor renormalization group and the higher-order second renormalization group methods have the seventh power. We perform benchmark calculations in the Ising model on the square lattice and show that the time-to-solution of the proposed algorithm is faster than that of other methods.
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.
An extension of the projected entangled-pair states (PEPS) algorithm to infinite systems, known as the iPEPS algorithm, was recently proposed to compute the ground state of quantum systems on an infinite two dimensional lattice. Here we investigate a modification of the iPEPS algorithm, where the environment is computed using the corner transfer matrix renormalization group (CTMRG) method, instead of using one-dimensional transfer matrix methods as in the original proposal. We describe a variant of the CTMRG that addresses different directions of the lattice independently, and use it combined with imaginary time evolution to compute the ground state of the two dimensional quantum Ising model. Near criticality, the modified iPEPS algorithm is seen to provide a better estimation of the order parameter and correlators.
We introduce a general corner transfer matrix renormalization group algorithm tailored to projected entangled-pair states on the triangular lattice. By integrating automatic differentiation, our approach enables direct variational energy minimization on this lattice geometry. In contrast to conventional approaches that map the triangular lattice onto a square lattice with diagonal next-nearest-neighbor interactions, our native formulation yields improved variational results at the same bond dimension. This improvement stems from a more faithful and physically informed representation of the entanglement structure in the tensor network and an increased number of variational parameters. We apply our method to the antiferromagnetic nearest-neighbor Heisenberg model on the triangular and kagome lattice, and benchmark our results against previous numerical studies.
We introduce a general corner transfer matrix renormalization group algorithm tailored to projected entangled-pair states on the triangular lattice. By integrating automatic differentiation, our approach enables direct variational energy minimization on this lattice geometry. In contrast to conventional approaches that map the triangular lattice onto a square lattice with diagonal next-nearest-neighbor interactions, our native formulation yields improved variational results at the same bond dimension. This improvement stems from a more faithful and physically informed representation of the entanglement structure in the tensor network and an increased number of variational parameters. We apply our method to the antiferromagnetic nearest-neighbor Heisenberg model on the triangular and kagome lattice, and benchmark our results against previous numerical studies.
Strongly correlated quantum systems are one of the most interesting and challenging phenomena in modern condensed matter physics. Unfortunately, the theoretical access even to their low energy properties is limited by the exponential growth of the Hilbert space in a many-body setting. To mitigate this computational barrier, tensor network methods are one of the most competitive approaches to classically simulate complex quantum systems. In this thesis we will discuss recent developments in the field of tensor networks to efficiently and accurately simulate two-dimensional quantum many-body states in the thermodynamic limit. Hereby, we introduce the variational update globally optimizing the wave-function, encoded as tensor network in the form of an infinite projected entangled pair state (iPEPS), in order to obtain ground states of many-body Hamiltonian models. We focus in the discussion on the well-established method to approximately contract this two-dimensional structure, the corner transfer matrix renormalization group (CTMRG), and discuss its native application on square lattice and triangular lattice based models. Furthermore, we report the development of a refined numerical algorithm, the split-CTMRG, to significantly reduce its computational complexity, pushing the limits of feasible simulations within modern hardware capabilities. Finally, the method advancements are applied to a currently widely studied model of frustrated magnetism, the Heisenberg maple-leaf antiferromagnet. We are able to find a large phase of vanishing magnetization in contradiction to previous works and, to the best of our knowledge, the lowest reported variational ground state energy at the isotropic point of the model. As an outlook, we describe a potential adaption of the previously developed methods to variationally simulate two-dimensional many-body thermal states.