IOP Publishing, Journal of Physics A: Mathematical and Theoretical, 36(41), p. 362001, 2008
DOI: 10.1088/1751-8113/41/36/362001
Full text: Download
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.