Efficient numerical integration of the adiabatic master equation

The capability to simulate open quantum systems under the adiabatic condition is critical for verifying the behavior of quantum annealers. While adiabaticity is a useful condition for theoretical analysis, straightforward numerical implementations of the adiabatic master equation (AME) suffer from the need to integrate an adiabatic but highly oscillatory quantum state over a very long period of time. We show that the AME in the adiabatic frame can be split in a way that prevents large superoperator evaluations and is amenable to geometric integration methods for split time-depedent problems.