The simulation of out-of-equilibrium dissipative quantum many body systems is a problem of fundamental interest to a number of fields in physics, ranging from condensed matter to cosmology. For unitary systems, tensor network methods have proved successful and extending these to open systems is a natural avenue for study. In particular, an important question concerns the possibility of approximating the critical dynamics of non-equilibrium systems with tensor networks. Here, we investigate this by performing numerical simulations of a paradigmatic quantum non-equilibrium system with an absorbing state: the quantum contact process. We consider the application of matrix product states and the time-evolving block decimation algorithm to simulate the time-evolution of the quantum contact process at criticality. In the Lindblad formalism, we find that the Heisenberg picture can be used to improve the accuracy of simulations over the Schrödinger approach, which can be understood by considering the evolution of operator-space entanglement. Furthermore, we also consider a quantum trajectories approach, which we find can reproduce the expected universal behaviour of key observables for a significantly longer time than direct simulation of the average state. These improved results provide further evidence that the universality class of the quantum contact process is not directed percolation, which is the class of the classical contact process.Schrödinger picture and Heisenberg picture. In the QCP, this asymmetry is particularly pronounced, and the Heisenberg picture does not display an absorbing state. It is then of interest to investigate any difference in performance between simulations in the two.Finally, the critical physics of the 1d QCP is similar to that of the CCP, in the sense that key observables display power-law behaviour but with different exponents. This means that, at criticality, different numerical methods or approaches to the dynamics can be compared by their ability to reproduce the expected power-laws. Furthermore, since the universality class of the QCP is currently debated, [7,9,11], it is of considerable interest in its own right to make estimates of critical exponents, comparing these with previous estimates and known cases.To simulate the non-equilibrium dynamics of the QCP, we apply matrix product states (MPSs) and the timeevolving block-decimation (TEBD) algorithm [13][14][15]. This algorithm belongs to a more general class of tensor network (TN) methods, well established for the simulation of closed quantum systems in 1d, which have also been applied to dissipative quantum systems previously in a number of cases [16][17][18][19][20][21][22][23][24]. In the context of studying dissipative quantum dynamics, a key question for TN methods is whether different approaches, such as quantum trajectories (QTs) as opposed to the Lindblad master equation, can lead to substantially different accuracies. This question has been explored previously in [16,17] and it has been suggested that in highentanglement ...