A stable toolkit method in quantum control

Recently the 'toolkit' discretization introduced to accelerate the numerical resolution of the time-dependent Schrödinger equation arising in quantum optimal control problems demonstrated good results on a large range of models. However, when coupling this class of methods with the so-called monotonically convergent algorithms, numerical instabilities affect the convergence of the discretized scheme. We present an adaptation of the 'toolkit' method which preserves the monotonicity of the procedure. The theoretical properties of the new algorithm are illustrated by numerical simulations.