We present a new variational method, based on the matrix product operator (MPO) ansatz, for finding the steady state of dissipative quantum chains governed by master equations of the Lindblad form. Instead of requiring an accurate representation of the system evolution until the stationary state is attained, the algorithm directly targets the final state, thus allowing for a faster convergence when the steady state is a MPO with small bond dimension. Our numerical simulations for several

more »

... tive spin models over a wide range of parameters illustrate the performance of the method and show that indeed the stationary state is often well described by a MPO of very moderate dimensions.