© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2310.08533 · Physical review. B./Physical review. B

Efficient representation of minimally entangled typical thermal states in two dimensions via projected entangled pair states

Aritra Sinha, Marek M. Rams, Jacek Dziarmaga

Published Jan 19, 2024• 10 citations

Abstract

The minimally entangled typical thermal states (METTS) are an ensemble of pure states, equivalent to the Gibbs thermal state, designed with an efficient tensor network representation in mind. In this article, we use the projected entangled pair states (PEPS) as their representation on a two-dimensional (2D) lattice. Unlike matrix product states (MPS), which for 2D systems are limited by an exponential computational barrier in the lattice size, PEPS provides a more tractable approach. To substantiate the prowess of PEPS in modeling METTS (dubbed PEPS-METTS), we benchmark it against the purification method in the context of the 2D quantum Ising model at its critical temperature. Our analysis reveals that PEPS-METTS achieves accurate results with significantly lower bond dimensions. We further corroborate this finding in the 2D Fermi-Hubbard model. At a technical level, we introduce an efficient zipper method to obtain PEPS boundary MPS needed to compute expectation values and perform sampling. The imaginary time evolution is done with the neighborhood tensor update.

DOI: 10.1103/physrevb.109.045136

Related papers

Minimally entangled typical thermal states with auxiliary matrix-product-state bases

Chia-Min Chung, Ulrich Schollwöck

Finite temperature problems in the strong correlated systems are important but challenging tasks. Minimally entangled typical thermal states (METTS) are a powerful method in the framework of tensor network methods to simulate finite temperature systems, including Fermions and frustrated spins which introduce a sign problem in the typical Monte Carlo methods. In this work, we introduce an extension of the METTS algorithm by using a new basis, the auxiliary matrix product state. This new basis achieves the pre-summation process in the partition function, and thus improve the convergence in the Monte Carlo samplings. The method also has the advantage of simulating the grand canonical ensemble in a computationally efficient way by employing good quantum numbers. We benchmark our method on the spin-$1/2$ XXZ model on the triangular lattice, and show that the new method outperforms the original METTS as well as the purification methods at sufficiently low temperature, the usual range of applications of METTS. The new method also naturally connects the METTS method to the purification method.

arXiv:1910.03329•Oct 8, 2019•2 citations

Variational approach to projected entangled pair states at finite temperature

Piotr Czarnik, Jacek Dziarmaga

The projected entangled pair state (PEPS) ansatz can represent a thermal state in a strongly correlated system. We introduce a variational algorithm to optimize this tensor network whose essential ingredient is an auxiliary tree tensor network (TTN). Since the full tensor environment is taken into account, with increasing bond dimension the PEPS-TTN ansatz provides the exact Gibbs state. Our presentation opens with a 1D version for a matrix product state (MPS-TTN) and then generalizes to PEPS-TTN in 2D. Benchmark results in the quantum Ising model are presented.

arXiv:1503.01077•Physical Review B•Jul 29, 2015•88 citations

Projected entangled pair states at finite temperature: Imaginary time evolution with ancillas

Piotr Czarnik, Łukasz Cincio, Jacek Dziarmaga

A projected entangled pair state (PEPS) with ancillas is evolved in imaginary time. This tensor network represents a thermal state of a two-dimensional (2D) lattice quantum system. A finite-temperature phase diagram of the 2D quantum Ising model in a transverse field is obtained as a benchmark application.

arXiv:1209.0454•Physical Review B•Dec 5, 2012•100 citations

Minimally entangled typical thermal state algorithms

Stoudenmire, E. M., Steven R. White

We discuss a method based on sampling minimally entangled typical thermal states (METTS) that can simulate finite temperature quantum systems with a computational cost comparable to ground state DMRG. Detailed implementations of each step of the method are presented, along with efficient algorithms for working with matrix product states and matrix product operators. We furthermore explore how properties of METTS can reveal characteristic order and excitations of systems and discuss why METTS form an efficient basis for sampling. Finally, we explore the extent to which the average entanglement of a METTS ensemble is minimal.

•May 28, 2010•91 citations

Minimally entangled typical thermal state algorithms

E M Stoudenmire, Steven R White

We discuss a method based on sampling minimally entangled typical thermal states (METTS) that can simulate finite temperature quantum systems with a computational cost comparable to the ground state density matrix renormalization group (DMRG). Detailed implementation of each step of the method is presented, along with efficient algorithms for working with matrix product states and matrix product operators. Furthermore, we explore how the properties of METTS can reveal characteristic order and excitations of systems and discuss why METTS form an efficient basis for sampling. Finally, we explore the extent to which the average entanglement of a METTS ensemble is minimal.

arXiv:1002.1305•New Journal of Physics•May 28, 2010•300 citations

Symmetric minimally entangled typical thermal states for canonical and grand-canonical ensembles

Moritz Binder, Thomas Barthel

Based on the density matrix renormalization group (DMRG), strongly correlated quantum many-body systems at finite temperatures can be simulated by sampling over a certain class of pure matrix product states (MPS) called minimally entangled typical thermal states (METTS). When a system features symmetries, these can be utilized to substantially reduce MPS computation costs. It is conceptually straightforward to simulate canonical ensembles using symmetric METTS. In practice, it is important to alternate between different symmetric collapse bases to decrease autocorrelations in the Markov chain of METTS. To this purpose, we introduce symmetric Fourier and Haar-random block bases that are efficiently mixing. We also show how grand-canonical ensembles can be simulated efficiently with symmetric METTS. We demonstrate these approaches for spin-$1/2$ $XXZ$ chains and discuss how the choice of the collapse bases influences autocorrelations as well as the distribution of measurement values and, hence, convergence speeds.

arXiv:1701.03872•Physical review. B./Physical review. B•May 22, 2017•25 citations

The minimally entangled typical thermal states (METTS) are an ensemble of pure states, equivalent to the Gibbs thermal state, designed with an efficient tensor network representation in mind. In this article, we use the projected entangled pair states (PEPS) as their representation on a two-dimensional (2D) lattice. Unlike matrix product states (MPS), which for 2D systems are limited by an exponential computational barrier in the lattice size, PEPS provides a more tractable approach. To substantiate the prowess of PEPS in modeling METTS (dubbed PEPS-METTS), we benchmark it against the purification method in the context of the 2D quantum Ising model at its critical temperature. Our analysis reveals that PEPS-METTS achieves accurate results with significantly lower bond dimensions. We further corroborate this finding in the 2D Fermi-Hubbard model. At a technical level, we introduce an efficient zipper method to obtain PEPS boundary MPS needed to compute expectation values and perform sampling. The imaginary time evolution is done with the neighborhood tensor update.

Finite temperature problems in the strong correlated systems are important but challenging tasks. Minimally entangled typical thermal states (METTS) are a powerful method in the framework of tensor network methods to simulate finite temperature systems, including Fermions and frustrated spins which introduce a sign problem in the typical Monte Carlo methods. In this work, we introduce an extension of the METTS algorithm by using a new basis, the auxiliary matrix product state. This new basis achieves the pre-summation process in the partition function, and thus improve the convergence in the Monte Carlo samplings. The method also has the advantage of simulating the grand canonical ensemble in a computationally efficient way by employing good quantum numbers. We benchmark our method on the spin-

1/2

XXZ model on the triangular lattice, and show that the new method outperforms the original METTS as well as the purification methods at sufficiently low temperature, the usual range of applications of METTS. The new method also naturally connects the METTS method to the purification method.

The projected entangled pair state (PEPS) ansatz can represent a thermal state in a strongly correlated system. We introduce a variational algorithm to optimize this tensor network whose essential ingredient is an auxiliary tree tensor network (TTN). Since the full tensor environment is taken into account, with increasing bond dimension the PEPS-TTN ansatz provides the exact Gibbs state. Our presentation opens with a 1D version for a matrix product state (MPS-TTN) and then generalizes to PEPS-TTN in 2D. Benchmark results in the quantum Ising model are presented.

We discuss a method based on sampling minimally entangled typical thermal states (METTS) that can simulate finite temperature quantum systems with a computational cost comparable to ground state DMRG. Detailed implementations of each step of the method are presented, along with efficient algorithms for working with matrix product states and matrix product operators. We furthermore explore how properties of METTS can reveal characteristic order and excitations of systems and discuss why METTS form an efficient basis for sampling. Finally, we explore the extent to which the average entanglement of a METTS ensemble is minimal.

We discuss a method based on sampling minimally entangled typical thermal states (METTS) that can simulate finite temperature quantum systems with a computational cost comparable to the ground state density matrix renormalization group (DMRG). Detailed implementation of each step of the method is presented, along with efficient algorithms for working with matrix product states and matrix product operators. Furthermore, we explore how the properties of METTS can reveal characteristic order and excitations of systems and discuss why METTS form an efficient basis for sampling. Finally, we explore the extent to which the average entanglement of a METTS ensemble is minimal.

Based on the density matrix renormalization group (DMRG), strongly correlated quantum many-body systems at finite temperatures can be simulated by sampling over a certain class of pure matrix product states (MPS) called minimally entangled typical thermal states (METTS). When a system features symmetries, these can be utilized to substantially reduce MPS computation costs. It is conceptually straightforward to simulate canonical ensembles using symmetric METTS. In practice, it is important to alternate between different symmetric collapse bases to decrease autocorrelations in the Markov chain of METTS. To this purpose, we introduce symmetric Fourier and Haar-random block bases that are efficiently mixing. We also show how grand-canonical ensembles can be simulated efficiently with symmetric METTS. We demonstrate these approaches for spin-

1/2

XXZ

chains and discuss how the choice of the collapse bases influences autocorrelations as well as the distribution of measurement values and, hence, convergence speeds.