© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1008.5126 · The Journal of Chemical Physics

Monotonically convergent optimization in quantum control using Krotov's method

Daniel M. Reich, Mamadou Ndong, Christiane P. Koch

Published Mar 9, 2012• 219 citations

Abstract

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.

DOI: 10.1063/1.3691827

Related papers

Accelerated monotonic convergence of optimal control over quantum dynamics

Tak‐San Ho, HERSCHEL A. RABITZ

The control of quantum dynamics is often concerned with finding time-dependent optimal control fields that can take a system from an initial state to a final state to attain the desired value of an observable. This paper presents a general method for formulating monotonically convergent algorithms to iteratively improve control fields. The formulation is based on a two-point boundary-value quantum control paradigm (TBQCP) expressed as a nonlinear integral equation of the first kind arising from dynamical invariant tracking control. TBQCP is shown to be related to various existing techniques, including local control theory, the Krotov method, and optimal control theory. Several accelerated monotonic convergence schemes for iteratively computing control fields are derived based on TBQCP. Numerical simulations are compared with the Krotov method showing that the new TBQCP schemes are efficient and remain monotonically convergent over a wide range of the iteration step parameters and the control pulse lengths, which is attributable to the trap-free character of the transition probability quantum dynamics control landscape.

Physical Review E•Aug 12, 2010•44 citations

Optimal control with accelerated convergence: Combining the Krotov and quasi-Newton methods

Reuven Eitan, Michael D. Mundt, David J. Tannor

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.

Physical Review A•May 24, 2011•106 citations

Convergence of the time-discretized monotonic schemes

Julien Salomon

Many numerical simulations in (bilinear) quantum control use the monotonically convergent Krotov algorithms (introduced by Tannor et al. [Time Dependent Quantum Molecular Dynamics (1992) 347–360]), Zhu and Rabitz [J. Chem. Phys. (1998) 385–391] or their unified form described in Maday and Turinici [J. Chem. Phys. (2003) 8191–8196]. In Maday et al. [Num. Math. (2006) 323–338], a time discretization which preserves the property of monotonicity has been presented. This paper introduces a proof of the convergence of these schemes and some results regarding their rate of convergence.

ESAIM Mathematical Modelling and Numerical Analysis•Jan 1, 2007•20 citations

Krotov Method in Optimal Quantum Control. I. Closed Systems

Oleg Vasilevich Morzhin, Alexander Nikolaevich Pechen

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.

•Sep 25, 2018•0 citations

New formulations of monotonically convergent quantum control algorithms

Yvon Maday, Gabriel Turinici

Most of the numerical simulation in quantum (bilinear) control have used one of the monotonically convergent algorithms of Krotov (introduced by Tannor et al.) or of Zhu and Rabitz. However, until now no explicit relationship has been revealed between the two algorithms in order to understand their common properties. Within this framework, we propose in this paper a unified formulation that comprises both algorithms and that extends to a new class of monotonically convergent algorithms. Numerical results show that the newly derived algorithms behave as well as (and sometimes better than) the well-known algorithms cited above.

The Journal of Chemical Physics•Apr 23, 2003•193 citations

Monotonically convergent algorithms for locally constraint quantum controls

Gabriel Turinici

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

•Feb 16, 2005•0 citations

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.

The control of quantum dynamics is often concerned with finding time-dependent optimal control fields that can take a system from an initial state to a final state to attain the desired value of an observable. This paper presents a general method for formulating monotonically convergent algorithms to iteratively improve control fields. The formulation is based on a two-point boundary-value quantum control paradigm (TBQCP) expressed as a nonlinear integral equation of the first kind arising from dynamical invariant tracking control. TBQCP is shown to be related to various existing techniques, including local control theory, the Krotov method, and optimal control theory. Several accelerated monotonic convergence schemes for iteratively computing control fields are derived based on TBQCP. Numerical simulations are compared with the Krotov method showing that the new TBQCP schemes are efficient and remain monotonically convergent over a wide range of the iteration step parameters and the control pulse lengths, which is attributable to the trap-free character of the transition probability quantum dynamics control landscape.

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.

Many numerical simulations in (bilinear) quantum control use the monotonically convergent Krotov algorithms (introduced by Tannor et al. [Time Dependent Quantum Molecular Dynamics (1992) 347–360]), Zhu and Rabitz [J. Chem. Phys. (1998) 385–391] or their unified form described in Maday and Turinici [J. Chem. Phys. (2003) 8191–8196]. In Maday et al. [Num. Math. (2006) 323–338], a time discretization which preserves the property of monotonicity has been presented. This paper introduces a proof of the convergence of these schemes and some results regarding their rate of convergence.

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.

Most of the numerical simulation in quantum (bilinear) control have used one of the monotonically convergent algorithms of Krotov (introduced by Tannor et al.) or of Zhu and Rabitz. However, until now no explicit relationship has been revealed between the two algorithms in order to understand their common properties. Within this framework, we propose in this paper a unified formulation that comprises both algorithms and that extends to a new class of monotonically convergent algorithms. Numerical results show that the newly derived algorithms behave as well as (and sometimes better than) the well-known algorithms cited above.

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