High order non-unitary split-step decomposition of unitary operators

Abstract

We propose a high order numerical decomposition of exponentials of hermitean operators in terms of a product of exponentials of simple terms, following an idea which has been pioneered by M. Suzuki, however implementing it for complex coefficients. We outline a convenient fourth order formula which can be written compactly for arbitrary number of noncommuting terms in the Hamiltonian and which is superiour to the optimal formula with real coefficients, both in complexity and accuracy. We show asymptotic stability of our method for sufficiently small time step and demonstrate its efficiency and accuracy in different numerical models.