Spectral functions with infinite projected entangled-pair states

Infinite projected entangled-pair states (iPEPS) provide a powerful tool to study two-dimensional strongly correlated systems directly in the thermodynamic limit. In this work, we extend the iPEPS toolbox by a method to efficiently evaluate nonequal time two-point correlators, enabling the computation of spectral functions. It is based on an iPEPS ansatz of the ground state in a large unit cell, with an operator applied in the center of the cell, which is evolved in real-time using the fast-full update method. At every time step, the two-point correlators within a cell are computed based on the corner transfer matrix renormalization group method. Benchmark results for the two-dimensional transverse field Ising model show that the main features of the dynamical structure factor can already be reproduced at relatively small bond dimensions and unit cell sizes. The results for the magnon dispersion are found to be in good agreement with previous data obtained with the iPEPS excitation ansatz.