Improved contraction of finite projected entangled pair states

TL;DR

This paper introduces an improved algorithm for contracting and optimizing finite projected entangled pair states (fPEPS), enhancing efficiency and scalability for complex quantum many-body systems.

Contribution

It details contraction patterns and combines controlled bond expansion with randomized SVD, advancing the computational methods for fPEPS.

Findings

Achieved efficient contraction with minimal memory usage.

Benchmark results for 8x8 Hubbard model systems.

Supported SU(2) symmetry with high bond dimensions.

Abstract

We present an improved version of the algorithm contracting and optimizing finite projected entangled pair states (fPEPS) in conjunction with projected entangled pair operators (PEPOs). Our work has two components to it. First, we explain in detail the characteristic contraction patterns that occur in fPEPS calculations and how to slice them such that peak memory occupation remains minimal while ensuring efficient parallel computation. Second, we combine controlled bond expansion [A. Gleis, J.-W. Li, and J. von Delft, Phys. Rev. Lett. 130, 246402 (2023)] with randomized singular value decomposition [V. Rokhlin, A. Szlam, and M. Tygert, SIAM J. Matrix Anal. Appl. (2009)] and apply it throughout the fPEPS algorithm. We present benchmark results for the Hubbard model for system sizes up to 8x8 and SU(2) symmetric bond dimension of up to D = 6 for PEPS bonds and $\chi$ = 500 for the environment bonds. Finally, we comment on the state and future of the fPEPS-PEPO framework.