Integrals of motion as slow modes in dissipative many-body operator dynamics

We consider Lindbladian operator dynamics in many-body quantum systems with one or more integrals of motion (IOM), subject to weak local dissipation. We demonstrate that IOMs with small support become slow modes of these dynamics, in the sense that their Frobenius norm decays more slowly compared to generic operators. As a result, the eigenoperators of such Lindbladians with slowest decay rates have a large overlap with the IOMs of the underlying Hamiltonian. We demonstrate this correspondence between slow modes and IOMs numerically for a number of many-body models, and further corroborate it with perturbative arguments. These results open up a new method for the identification of IOMs, and provide insights into the dissipative many-body dynamics.