Chris is right about Magnus propagators being the method of choice for QM. See also this earlier post: How to efficiently exponentiate a complex (symmetric) tridiagonal matrix?
Are you sure you want this Hamiltonian matrix though?
This is non-Hermitian, which is very untypical in Markovian (=normal) QM.