Strang splitting schemes for N-level Bloch models

Abstract: We extend to the N-level Bloch model the splitting scheme which use exact numerical solutions of sub-equations. These exact solutions involve matrix exponentials which we want to avoid to calculate at each time step. We use Newton interpolation to reduce the computational cost. The resulting scheme is nonstandard and preserves all qualitative properties of the Bloch equations. We show numerical simulations to compare this approach with a few other schemes.