Quantum quenches in two spatial dimensions using chain array matrix product states

We describe a method for simulating the real time evolution of extended quantum systems in two dimensions (2D). The method combines the benefits of integrability and matrix product states in one dimension to avoid several issues that hinder other applications of tensor based methods in 2D. In particular, it can be extended to infinitely long cylinders. As an example application we present results for quantum quenches in the 2D quantum [(2 + 1)-dimensional] Ising model. In quenches that cross a phase boundary we find that the return probability shows nonanalyticities in time.