P18: Tensor network algorithms
For many fundamental models of strongly correlated quantum many-body systems, a thorough understanding of their physical properties has still not been achieved. This holds in particular true for strongly correlated fermion systems, described e.g. by the Hubbard model, and for many basic models of quantum magnetism, in particular for systems with (geometric) frustration. In recent years, a series of new algorithmic approaches to such systems have been put forward, inspired to a large degree by insights gained from quantum information science on the entanglement properties of quantum many-body states. These methods provide a new approach to study basic models of strongly correlated quantum many-body systems. The basic ingredient to these methods are complete sets of quantum states that are defined in terms of networks of locally contracted tensors, which generalize the previously introduced matrix product states (MPS) (1),(2) that also underlie the density matrix renormalization group method used in P5 and P12. The MPS have thus been employed already very successfully to study lattice-based quantum many-body systems in one spatial dimension (3). The more general class of states, also called projected entangled pair states (PEPS) (4),(5), allow for an efficient representation of the entanglement properties of correlated quantum
systems also beyond one spatial dimension. Currently, they are intensively employed to study two-dimensional quantum systems. In the first part of this project, we will assess the reliability of PEPS-based methods in treating specific two-dimensional frustrated quantum spin systems (6),(7) and quantum phase transitions(8). For this purpose, a native version of the PEPS algorithms will be implemented, based on the itensor library.
In particular, we will perform a detailed comparison of the PEPS-based results with previous unbiased quantum Monte Carlo results that exist for these specific systems(9),(10). Of special interest here will be the non-variational aspects, that enter the iterative tensor-contraction. We also plan to examine these aspects from an analytical perspective. These studies will allow us to judge the robustness of these new approaches. It will also be important to explore in details the scaling behavior of the numerical data with the various control parameters of the numerical algorithms in order to provide robust means of extrapolations. In a second step, we plan to apply also optimization methods to the tensor network states employing Monte Carlo sampling, and extend these methods to the treatment of
fermion systems (11)-(15).
The central goal of this project thus is to establish tensor network based algorithms as a reliable approach to the study of complex quantum many-body systems, in particular for studying ground state properties such as e.g. in project P2.
(1) S. Ostlund and S. Rommer, Phys. Rev. Lett. 75, 3537 (1995)
(2) F. Verstraete, D. Porras and J. I. Cirac, Phys. Rev. Lett. 93 227205 (2004)
(3) U. Schollwöck, Ann. of Phys. 326, 96 (2011)
(4) F. Verstraete and J. I. Cirac, cond-mat/0407066 (2004)
(5) F. Verstraete, V. Murg and J. I. Cirac, Adv. Phys. 57, 143 (2008)
(6) Z. Y. Meng and S. Wessel, Phys. Rev. B 78, 224416 (2008)
(7) S. Wessel and M. Troyer, Phys. Rev. Lett. 95, 127205 (2005)
(8) L. Fritz, R.L. Doretto, S. Wessel, S. Wenzel, S. Burdin and M. Vojta, Phys. Rev. B 83, 174416 (2011)
(9) F. Alet, S. Wessel and M. Troyer, Phys. Rev. E 71, 036706 (2005)
(10) O. F. Syljuåsen and A. W. Sandvik, Phys. Rev. E 66, 046701 (2002)
(11) N. Schuch, M. M. Wolf, F. Verstraete and J. I. Cirac, Phys. Rev. Lett. 100, 040501 (2008)
(12) F. Mezzacapo, N. Schuch, M. Boninsegni and J. I. Cirac, New J. Phys. 11, 083026 (2009)
(13) L. Wang, I. Pižorn and F. Verstraete, Phys. Rev. B 83, 134421 (2011)
(14) A. W. Sandvik and G. Vidal, Phys. Rev. Lett. 99, 220602 (2007)
(15) C. P. Chou, F. Pollmann and T. K. Lee, Phys. Rev. B 86, 041105 (2012)