This chapter is devoted to numerical optimization methods with a focus on solution schemes that are used in the framework of quantum control problems. We first discuss classical optimization schemes that require the computation of reduced gradients and usually provide first-order convergence to a local minimum. These methods are supported by a large set of theoretical results that provide guarantee of convergence to local minima. The second section is devoted to monotonic schemes that originate from the works of Vadim Fedorovitch Krotov. While these methods remain popular in the quantum control application communities, in other frameworks the methods of choice are those discussed in Section 6.1.
This chapter is devoted to numerical optimization methods with a focus on solution schemes that are used in the framework of quantum control problems. We first discuss classical optimization schemes that require the computation of reduced gradients and usually provide first-order convergence to a local minimum. These methods are supported by a large set of theoretical results that provide guarantee of convergence to local minima. The second section is devoted to monotonic schemes that originate from the works of Vadim Fedorovitch Krotov. While these methods remain popular in the quantum control application communities, in other frameworks the methods of choice are those discussed in Section 6.1.
The non-linear optimization method developed by A. Konnov and V. Krotov [Autom. Remote Cont. (Engl. Transl.) 60, 1427 (1999)] has been used previously to extend the capabilities of optimal control theory from the linear to the non-linear Schrödinger equation [S. E. Sklarz and D. J. Tannor, Phys. Rev. A 66, 053619 (2002)]. Here we show that based on the Konnov-Krotov method, monotonically convergent algorithms are obtained for a large class of quantum control problems. It includes, in addition to nonlinear equations of motion, control problems that are characterized by non-unitary time evolution, nonlinear dependencies of the Hamiltonian on the control, time-dependent targets, and optimization functionals that depend to higher than second order on the time-evolving states. We furthermore show that the nonlinear (second order) contribution can be estimated either analytically or numerically, yielding readily applicable optimization algorithms. We demonstrate monotonic convergence for an optimization functional that is an eighth-degree polynomial in the states. For the "standard" quantum control problem of a convex final-time functional, linear equations of motion and linear dependency of the Hamiltonian on the field, the second-order contribution is not required for monotonic convergence but can be used to speed up convergence. We demonstrate this by comparing the performance of first- and second-order algorithms for two examples.
Mathematical problems of optimal control of quantum systems attract high interest in connection with both fundamental questions and existing and prospective applications. An important problem is the development of methods for constructing controls for quantum systems. One of the widely used methods is the Krotov method. The article is the first part of the review devoted to mathematical aspects of this method for optimal control of quantum systems. Beyond quantum control this method was suggested in the articles [V.F. Krotov, I.N. Feldman, Eng. Cybern. 1983. Vol. 21, 123], [A.I. Konnov, V.F. Krotov, Automat. Remote Control, 1999], and then adopted to optimal quantum control in papers by V.A. Kazakov, V.F. Krotov, R. Kosloff, C.P. Koch, D.J. Tannor, and others. The review considers this method, including cases with constraints on control spectrum and on states of the quantum system, various ways for choosing the improvement function, in particular, linear or linear-quadratic form, etc. We also discuss GRAPE (GRadient Ascent Pulse Engineering) and CRAB (Chopped Random-Basis) methods, as well as combined methods such as the method suggested by I. Maday and G. Turinici for generalizing the 1st order Krotov method and the Zhu-Rabitz method. Nowadays, various versions of the Krotov method are used as a key tool in optimal quantum control, and the review describes applications of this method to mathematical control problems for molecular dynamics, Bose-Einstein condensate, nuclear-magnetic resonance, quantum gates generation, etc. The first part of the review is devoted to closed quantum systems.
One of the most popular methods for solving numerical optimal control problems is the Krotov method, adapted for quantum control by Tannor and coworkers. The Krotov method has the following three appealing properties: (1) monotonic increase of the objective with iteration number, (2) no requirement for a line search, leading to a significant savings over gradient (first-order) methods, and (3) macrosteps at each iteration, resulting in significantly faster growth of the objective at early iterations than in gradient methods where small steps are required. The principal drawback of the Krotov method is slow convergence at later iterations, which is particularly problematic when high fidelity is desired. We show here that, near convergence, the Krotov method degenerates to a first-order gradient method. We then present a variation on the Krotov method that has all the advantages of the original Krotov method but with significantly enhanced convergence (second-order or quasi-Newton) as the optimal solution is approached. We illustrate the method by controlling the three-dimensional dynamics of the valence electron in the Na atom.
The control of matter and energy at a fundamental level will be a cornerstone of new technologies for years to come. This idea is exemplified in a distilled form by controlling the dynamics of quantum mechanical systems via a time-dependent potential. The contributions detailed within this work focus on the computational aspects of formulating and solving quantum control problems efficiently.
The accurate numerical computation of optimal controls of infinite dimensional quantum control problems is a very difficult task that requires to take into account the features of the original infinite dimensional problem. An important issue is the choice of the functional space where the minimization process is defined. A systematic comparison of L2 versus H1 based minimization shows that the choice of the appropriate functional space matters and has many consequences in the implementation of some optimization techniques.
A matrix-free cascadic BFGS algorithm is introduced in the L 2 and H1 settings and it is demonstrated that the choice of H1 over L2 results in a substantial performance and robustness increase. A comparison between optimal control resulting from function space minimization and the control obtained by minimization over Chebyshev and POD basis function coefficients is presented.
A theoretical and computational framework is presented to obtain accurate controls for fast quantum state transitions that are needed in a host of applications such as nano electronic devices and quantum computing. This method is based on a reduced Hessian Krylov-Newton scheme applied to a norm-preserving discrete model of a dipole quantum control problem. The use of second-order numerical methods for solving the control problem is justified proving existence of optimal solutions and analyzing first- and second-order optimality conditions. Criteria for the discretization of the non-convex optimization problem and for the formulation of the Hessian are given to ensure accurate gradients and a symmetric Hessian. Robustness of the Newton approach is obtained using a globalization strategy with a robust linesearch procedure. Results of numerical experiments demonstrate that the Newton approach presented in this dissertation is able to provide fast and accurate controls for high-energy state transitions.
Control of bound-to-bound and bound-to-continuum transitions in open quantum systems and vector field control of two-dimensional systems is presented. An efficient space-time spectral discretization of the time-dependent Schrodinger equation and preconditioning strategy for a fast approximate solution with Krylov methods is outlined.
The problem of finding the optimal control in numerical computer simulations of quantum control phenomena is usually addressed through the introduction of monotonically convergent algorithms that are guaranteed to improve the cost functional at each step. A recent extension of these algorithms implements a search for a control with given bounds. Within this context, this paper will present a generalisation of this algorithm to a more general class of local in time constraints. Numerical results that illustrate some of the properties of the algorithms introduced are given. KEY WORDS quantum control, bilinear control, monotonic algorithms, numerical methods, non-linear control 1
The purpose of this chapter is to revisit the quantum control problems and the solution methodologies discussed in the previous chapters presenting results of numerical experiments. In the following sections, we consider different model problems, discussed in Chapters 2, 3, and 4, and solution strategies based on the discretization schemes discussed in Chapter 5 and on the optimization methods illustrated in Chapter 6. We also provide some comparison of the strategies to provide information about which method is suitable in which context. Discussion on how to choose parameters characterizing the methods provide valuable insights for the practical implementation of quantum optimal control solvers.