Non-unitarity maximizing unraveling of open quantum dynamics

The dynamics of many-body quantum states in open systems is commonly numerically simulated by unraveling the density matrix into pure-state trajectories. In this work, we introduce a new unraveling strategy that can adaptively minimize the averaged entanglement in the trajectory states. This enables a more efficient classical representation of trajectories using matrix product decompositions. Our new approach is denoted non-unitarity maximizing unraveling (NUMU). It relies on the idea that adaptively maximizing the averaged non-unitarity of a set of Kraus operators leads to a more efficient trajectory entanglement destruction. Compared to other adaptive entanglement lowering algorithms, NUMU is computationally inexpensive. We demonstrate its utility in large-scale simulations with random quantum circuits. NUMU lowers runtimes in practical calculations, and it also provides new insight on the question of classical simulability of quantum dynamics. We show that for the quantum circuits considered here, unraveling methods are much less efficient than full matrix product density operator simulations, hinting to a still large potential for finding more advanced adaptive unraveling schemes.