Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review E

Accelerated monotonic convergence of optimal control over quantum dynamics

Tak‐San Ho, HERSCHEL A. RABITZ

Published Aug 12, 2010• 44 citations
Publisher

Abstract

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.

DOI
10.1103/physreve.82.026703

Related papers

Monotonically convergent optimization in quantum control using Krotov's method

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

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.

arXiv:1008.5126•The Journal of Chemical Physics•Mar 9, 2012•219 citations

Efficient algorithms for optimal control of quantum dynamics: the Krotov method unencumbered

S. G. Schirmer, Pierre de Fouquieres

Efficient algorithms are of fundamental importance for the discovery of optimal control designs for coherent control of quantum processes. One important class of algorithms is sequential update algorithms, generally attributed to Krotov. Although widely and often successfully used, the associated theory is often involved and leaves many crucial questions unanswered, from the monotonicity and convergence of the algorithm to discretization effects, leading to the introduction of ad hoc penalty terms and suboptimal update schemes detrimental to the performance of the algorithm. We present a general framework for sequential update algorithms including specific prescriptions for efficient update rules with inexpensive dynamic search length control, taking into account discretization effects and eliminating the need for ad hoc penalty terms. The latter, despite being necessary for regularizing the problem in the limit of infinite time resolution, i.e. the continuum limit, are shown to be undesirable and unnecessary in the practically relevant case of finite time resolution. Numerical examples show that the ideas underlying many of these results extend even beyond what can be rigorously proved.

arXiv:1103.5435•New Journal of Physics•Jul 22, 2011•59 citations

On the relationship between the local tracking procedures and monotonic schemes in quantum optimal control

Julien Salomon, Gabriel Turinici

Numerical simulations of (bilinear) quantum control often rely on either monotonically convergent algorithms or tracking schemes. However, despite their mathematical simplicity, very limited intuitive understanding exists at this time to explain the former type of algorithms. Departing from the usual mathematical formalization, we present in this paper an interpretation of the monotonic algorithms as finite horizon, local in time, tracking schemes. Our purpose is not to present a new class of procedures but rather to introduce the necessary rigorous framework that supports this interpretation. As a by-product we show that at each instant, estimates of the future quality of the current control field are available and used in the optimization. When the target is expressed as reaching a prescribed final state, we also present an intuitive geometrical interpretation as the minimization of the distance between two correlated trajectories: one starting from the given initial state and the other backward in time from the target state. As an illustration, a stochastic monotonic algorithm is introduced. Numerical discretizations of the two procedures are also presented.

The Journal of Chemical Physics•Feb 14, 2006•19 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

Monotonically convergent algorithms for solving quantum optimal control problems of a dynamical system nonlinearly interacting with a control

Yukiyoshi Ohtsuki, Kazuyuki Nakagami

We develop a family of monotonically convergent iterative algorithms for solving a wide class of optimal control problems in which a dynamical system interacts nonlinearly with a control. The key idea is to divide a control into identical components and to introduce auxiliary steps to update each component at every iteration step. The algorithms are proved to exhibit monotonic convergence, which is also numerically confirmed through a case study of the control of molecular orientation. The numerical results show that high-quality solutions can be obtained by using the present algorithms.

Physical Review A•Mar 18, 2008•43 citations

A stable toolkit method in quantum control

Mohamed Belhadj, Julien Salomon, Gabriel Turinici

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.

Journal of Physics A Mathematical and Theoretical•Aug 1, 2008•8 citations
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.
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.
Efficient algorithms are of fundamental importance for the discovery of optimal control designs for coherent control of quantum processes. One important class of algorithms is sequential update algorithms, generally attributed to Krotov. Although widely and often successfully used, the associated theory is often involved and leaves many crucial questions unanswered, from the monotonicity and convergence of the algorithm to discretization effects, leading to the introduction of ad hoc penalty terms and suboptimal update schemes detrimental to the performance of the algorithm. We present a general framework for sequential update algorithms including specific prescriptions for efficient update rules with inexpensive dynamic search length control, taking into account discretization effects and eliminating the need for ad hoc penalty terms. The latter, despite being necessary for regularizing the problem in the limit of infinite time resolution, i.e. the continuum limit, are shown to be undesirable and unnecessary in the practically relevant case of finite time resolution. Numerical examples show that the ideas underlying many of these results extend even beyond what can be rigorously proved.
Numerical simulations of (bilinear) quantum control often rely on either monotonically convergent algorithms or tracking schemes. However, despite their mathematical simplicity, very limited intuitive understanding exists at this time to explain the former type of algorithms. Departing from the usual mathematical formalization, we present in this paper an interpretation of the monotonic algorithms as finite horizon, local in time, tracking schemes. Our purpose is not to present a new class of procedures but rather to introduce the necessary rigorous framework that supports this interpretation. As a by-product we show that at each instant, estimates of the future quality of the current control field are available and used in the optimization. When the target is expressed as reaching a prescribed final state, we also present an intuitive geometrical interpretation as the minimization of the distance between two correlated trajectories: one starting from the given initial state and the other backward in time from the target state. As an illustration, a stochastic monotonic algorithm is introduced. Numerical discretizations of the two procedures are also presented.
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.
We develop a family of monotonically convergent iterative algorithms for solving a wide class of optimal control problems in which a dynamical system interacts nonlinearly with a control. The key idea is to divide a control into identical components and to introduce auxiliary steps to update each component at every iteration step. The algorithms are proved to exhibit monotonic convergence, which is also numerically confirmed through a case study of the control of molecular orientation. The numerical results show that high-quality solutions can be obtained by using the present algorithms.
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.