Abstract
Quantum Optimal Control (QOC) is the field devoted to the production of external control protocols that actively guide quantum dynamics. Solutions to QOC problems were shown to constitute continuous submanifolds of control space. A solution navigation method exploiting this property to achieve secondary features in the control protocols was proposed [Larocca \textit{et al}, arXiv:1911.07105]. Originally, the technique involved the computation of the exact Hessian matrix. In this paper, we show that the navigation can be alternatively performed with a finite-difference approximation scheme, thus enabling the application of this procedure in systems where computing the exact Hessian is out of reach.