Improved contraction of finite projected entangled pair states

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, with randomized singular value decomposition V. Rokhlin, A. Szlam, and M. Tygert, and apply it throughout the fPEPS algorithm. We present benchmark results for the Hubbard model for system sizes up to 8 × 8 and SU(2) symmetric bond dimension of up to D = 6 for PEPS bonds and χ = 500 for the environment bonds. Finally, we comment on the state and future of the fPEPS-PEPO framework.