We analyse in this paper the Lyapunov trajectory tracking of the Schrödinger equation for a second order coupling operator. We present a theoretical convergence result; for situations not covered by the theoretical result we propose a numerical approach that is tested and works well in practice. oui