Gradient methods for variational optimization of projected entangled-pair states

The authors present a new method for variationally optimizing projected entangled-pair states (PEPS) in the thermodynamic limit based on the global energy functional and novel tensor contraction strategies. It is successfully benchmarked using Ising and Heisenberg Hamiltonians. The method is completely independent from previous PEPS implementations and improves on the variational energies and order parameters as compared to other state-of-the-art PEPS simulations. The algorithm presented here should prove valuable in simulating ground-state properties of strongly correlated quantum spin systems, can be used with other variational ansatzes, and can be applied to tensor network states in general.