Numerical simulation of quantum nonequilibrium phase transitions without finite size effects

Classical (1+1)-dimensional (D) cellular automata, as for instance Domany-Kinzel cellular automata, are paradigmatic systems for the study of nonequilibrium phenomena. Such systems evolve in discrete time steps, and are thus free of time-discretization errors. Moreover, they display nonequilibrium phase transitions which can be studied by simulating the evolution of an initial seed. At any finite time, this has support only on a finite light cone. Thus, essentially numerically exact simulations free of finite-size errors or boundary effects are possible, leading to high-accuracy estimates of critical exponents. Here, we show how similar advantages can be gained in the quantum regime: The many-body critical dynamics occurring in (1+1)D quantum cellular automata with an absorbing state can be studied directly on an infinite lattice when starting from seed initial conditions. This can be achieved efficiently by simulating the dynamics of an associated one-dimensional, nonunitary quantum cellular automaton using tensor networks. We apply our method to a model introduced recently and find accurate values for universal exponents, suggesting that this approach can be a powerful tool for precisely studying nonequilibrium universal physics in quantum systems.