Abstract We consider splitting methods with complex coefficients to approximate the solutions of parabolic equations. We show resorting to complex times makes it possible to construct high order composition schemes that remain applicable for non-reversible operators