Numerical time propagation of semi-linear equations of the Schrödinger type can be performed by the use of exponential integrators. The main difficulty for efficient implementation of this type of schemes lies in the evaluation of $$\varphi $$ φ -functions of a matrix argument. We develop a Chebyshev series approximation for these functions and propose a simple algorithm for the evaluation of the series coefficients. The domain of convergence of the series is consistent with the spectrum of Schrödinger type operators. This approximation is shown to be accurate and performs favorably in comparison to other state of the art methods for approximation of $$\varphi $$ φ -functions.