  Infinite projected entangled pair states (iPEPS), the tensor network ansatz
for two-dimensional systems in the thermodynamic limit, already provide
excellent results on ground-state quantities using either imaginary-time
evolution or variational optimisation. Here, we show (i) the feasibility of
real-time evolution in iPEPS to simulate the dynamics of an infinite system
after a global quench and (ii) the application of disorder-averaging to obtain
translationally invariant systems in the presence of disorder. To illustrate
the approach, we study the short-time dynamics of the square lattice Heisenberg
model in the presence of a bi-valued disorder field.
