Abstract
Open quantum systems host a wide range of intriguing phenomena, yet their simulation on well-controlled quantum devices is challenging, owing to the exponential growth of the Hilbert space and the inherently non-unitary nature of the dynamics. Here we propose and experimentally demonstrate a variational quantum algorithm capable of scalable simulation of non-unitary many-body dissipative dynamics. The algorithm builds on the framework of linear combination of Hamiltonian simulation, which converts non-unitary dynamics into a weighted sum of unitary evolutions. With the further introduction of a simplified quantum circuit for loss-function evaluation, our scheme is suitable for near-term quantum hardware, with the circuit depth independent of the simulation time. We illustrate our scheme by simulating the collective dynamics of a dissipative transverse Ising model, as well as an interacting Hatano-Nelson model, on the superconducting quantum processor Wukong. Our work underlines the capability of noisy intermediate-scale quantum devices in simulating dissipative many-body dynamics and represents a step forward in exploiting their potential for solving outstanding physical problems.