Stochastic phase-space simulation of multimode cat states via the positive-P representation

We present a comprehensive study of the transient dynamics of multimode Schrödinger cat states in dissipatively coupled resonator arrays using the positive-P phase-space method. By employing the positive-P representation, we derive the exact stochastic differential equations governing the system's dynamics, enabling the simulation of system sizes significantly larger than those accessible via direct master equation simulation. We demonstrate the utility of this method by simulating transient dynamics for networks up to N=21 sites. Furthermore, we critically examine the method's usefulness and limitations, specifically highlighting the computational instability encountered when estimating the state parity in the systems. Our results provide a pathway for scalable simulations of non-Gaussian states in large open quantum systems.