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 A

Optimal control of quantum-mechanical systems: Existence, numerical approximation, and applications

Anthony P. Peirce, M. Dahleh, HERSCHEL A. RABITZ

Published Jun 1, 1988• 890 citations
Publisher

Abstract

The optimal control of the path to a specified final state of a quantum-mechanical system is investigated. The problem is formulated as a minimization problem over appropriate function spaces, and the well-posedness of this problem is is established by proving the existence of an optimal solution. A Lagrange-multiplier technique is used to reduce the problem to an equivalent optimization problem and to derive necessary conditions for a minimum. These necessary conditions form the basis for a gradient iterative procedure to search for a minimum. A numerical scheme based on finite differences is used to reduce the infinite-dimensional minimization problem to an approximate finite-dimensional problem. Numerical examples are provided for final-state control of a diatomic molecule represented by a Morse potential. Within the context of this optimal control formulation, numerical results are given for the optimal pulsing strategy to demonstrate the feasibility of wave-packet control and finally to achieve a specified dissociative wave packet at a given time. The optimal external optical fields generally have a high degree of structure, including an early time period of wave-packet phase adjustment followed by a period of extensive energy deposition to achieve the imposed objective. Constraints on the form of the molecular dipole (e.g., a linear dipole) are shown to limit the accessibility (i.e., controllability) of certain types of molecular wave-packet objectives. The nontrivial structure of the optimal pulse strategies emphasizes the ultimate usefulness of an optimal-control approach to the steering of quantum systems to desired objectives.

DOI
10.1103/physreva.37.4950

Related papers

Optimal control of uncertain quantum systems

M. Dahleh, Anthony P. Peirce, HERSCHEL A. RABITZ

The design of optimal final-state controllers of quantum-mechanical systems that are insensitive to errors in the molecular Hamiltonian or to errors in the initial state of the system is considered. Control arises through the interaction of the system with an external field; the goal is optimal design of these latter fields for various physical objectives in the presence of system uncertainty. Sensitivity to modeling errors and other uncertainties in the molecular Hamiltonian is minimized by considering averaged costs for a family of Hamiltonian functions ${\mathit{H}}_{0}$(\ensuremath{\alpha}) indexed by the random variable \ensuremath{\alpha} taking values on a compact set in Euclidean space. Similarly, sensitivity of the optimal control to the initial state is minimized by viewing the initial condition as a Hilbert-space-valued random variable and considering an optimization problem with a cost functional that is averaged over the class of initial conditions. A precise formulation of the control problem is given, and its well-posedness is established. Cost propagators are defined to display the dependence of the performance index on the initial conditions explicitly, which allows analytic averaging of initial conditions.The constrained optimization problem is reduced to an unconstrained optimization problem by the introduction of Lagrange-multiplier operators. Necessary conditions for the unconstrained problem provide the basis for a gradient search for an optimal solution. Finite-difference schemes are utilized to provide a numerical approximation of the optimal control problem. Numerical examples are given for final-state control of a diatomic molecule represented by a Morse potential illustrating design for systems with initial-phase uncertainty and parametric uncertainty. The resultant insensitive controllers execute different strategies depending on the design requirements. The controller designed to be insensitive to errors in the initial phases adopts a strategy of phase imprinting during the initial stages of the control interval to compensate for a lack of knowledge of the initial phases. It is also shown that it is not possible to coerce a system from a state with completely random initial phases to a correlated state using the class of averaged controllers considered here. The controller designed to be insensitive to parametric Hamiltonian errors adopts a strategy of amplitude restraint to prevent the wave packet from taking significant excursions into the regions where the potential is uncertain. The interesting structure exhibited by the controllers in response to the different design requirements, and the superior performance of the insensitive controllers when compared with controllers designed at nominal phases and parameters, illustrate the usefulness of the cost-averaging technique for design in the presence of uncertainties.

Physical Review A•Aug 1, 1990•147 citations

Optimally controlled quantum molecular dynamics: A perturbation formulation and the existence of multiple solutions

Metin Demіralp, HERSCHEL A. RABITZ

This work considers optimal control of quantum-mechanical systems within the framework of perturbation theory with respect to the controlling optical electric field. The control problem is expressed in terms of a cost functional including the physical objective, the penalties, and constraints. The resultant nonlinear variational equations are linearized by considering the lowest-order term in an expansion in powers of the optical-field strength. The optical field is found to satisfy a linear integral equation, and the solution may be expressed in terms of a generalized eigenvalue problem associated with the corresponding kernel. A full determination of the field is specified through the solution to the integral equation and the roots of an accompanying linearized spectral equation for a characteristic multiplier parameter. Each discrete value of the latter parameter corresponds to a particular solution to the variational equations. As a result, it is argued that under very general conditions there will be a denumerably infinite number of solutions to well-posed quantum-mechanical optimal-control problems.

Physical Review A•Feb 1, 1993•99 citations

Algebraic and computational aspects of quantum control and applications.

I. C. H. Pullen

We investigate the problem of control of quantum systems. Various notions of controllability of quantum systems are discussed and Lie algebraic techniques are employed to derive necessary and sufficient conditions for the resulting degrees of controllability. The results are employed to study the degree of controllability for various model systems. The problem of optimal control is formulated as a problem of maximisation of the expectation value of an observable of the system subject to dynamical constraints and costs. A class of iterative algorithms is discussed, and shown to converge to a solution of the Euler-Lagrange equations, which is a necessary condition for optimality of a solution. A numerical realisation of the algorithm is presented and potential implementation of the resulting control fields in the laboratory is discussed. The algorithm is applied to find optimal control fields for a variety of applications including selective and simultaneous excitation of individual quantum dots in a globally addressed ensemble, creation of superposition states, and selective excitation of certain vibrational modes of a simple molecule such as Hydrogen Fluoride. In each case we systematically explore the parameter space of the algorithm to select solutions with desirable physical features from the multitude of possible choices. As there are many possible optimal control fields which give similar yields the question of whether some are more robust than others is investigated. A number of optimal fields are analysed quantitatively in terms of robustness with regard to temporal and spectral noise affecting the fields, uncertainty in the model parameters and the effect of experimental limitations such as restrictions on the bandwidth of the fields. It is shown that whilst the solutions tend to be sensitive to uncertainty in the energy levels chosen for the model, they appear to be surprisingly tolerant to even severe bandwidth limitations.

•Jan 1, 2006•0 citations

Computational techniques for optimal control of quantum systems

Evangelos A. Coutsias, G. von Winckel

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.

•Jan 1, 2011•0 citations

Robust optimal control of quantum molecular systems in the presence of disturbances and uncertainties

Hong Zhang, HERSCHEL A. RABITZ

Recently, systematic optimal control theory has been successfully applied to design optical fields which act on a molecule to guide or transform it into a specifically prescribed state. Under realistic laboratory conditions and for real molecules, it is inevitable that there will exist disturbances and uncertainties in the control optical field, in the molecular Hamiltonian, and in the field-coupling coefficient (e.g., the dipole moment), etc. Therefore, it is important to design an optimal field which is robust to disturbances. In this paper, we provide a method through which the worst possible disturbance to the control process, under certain possible constraints, can be calculated and a corresponding optimal control field designed to be robust to the disturbances. Under disturbances of small amplitude, a perturbation solution to this problem can be applied. Illustrations are given for control of the wave-packet evolution of a diatomic molecule in the presence of disturbances.

Physical Review A•Apr 1, 1994•69 citations

A Simplified Approach to Optimally Controlled Quantum Dynamics

Jair Botina, Herschel Rabitz, N. K. Rahman

A new formalism for the optimal control of quantum mechanical physical observables is presented. This approach is based on an analogous classical control technique reported previously. [1] Quantum Lagrange multiplier functions are used to preserve a chosen subset of the observable dynamics of interest. As a result, a corresponding small set of Lagrange multipliers needs to be calculated and they are only a function of time. This is a considerable simplification over traditional quantum optimal control theory. [2] The success of the new approach is based on taking advantage of the multiplicity of solutions to virtually any problem of quantum control to meet a physical objective. A family of such simplified formulations is introduced and numerically tested. Results are presented for these algorithms and compared with previous reported work on a model problem for selective unimolecular reaction induced by an external optical electric field. 1 I.

•Jan 1, 1996•18 citations
The optimal control of the path to a specified final state of a quantum-mechanical system is investigated. The problem is formulated as a minimization problem over appropriate function spaces, and the well-posedness of this problem is is established by proving the existence of an optimal solution. A Lagrange-multiplier technique is used to reduce the problem to an equivalent optimization problem and to derive necessary conditions for a minimum. These necessary conditions form the basis for a gradient iterative procedure to search for a minimum. A numerical scheme based on finite differences is used to reduce the infinite-dimensional minimization problem to an approximate finite-dimensional problem. Numerical examples are provided for final-state control of a diatomic molecule represented by a Morse potential. Within the context of this optimal control formulation, numerical results are given for the optimal pulsing strategy to demonstrate the feasibility of wave-packet control and finally to achieve a specified dissociative wave packet at a given time. The optimal external optical fields generally have a high degree of structure, including an early time period of wave-packet phase adjustment followed by a period of extensive energy deposition to achieve the imposed objective. Constraints on the form of the molecular dipole (e.g., a linear dipole) are shown to limit the accessibility (i.e., controllability) of certain types of molecular wave-packet objectives. The nontrivial structure of the optimal pulse strategies emphasizes the ultimate usefulness of an optimal-control approach to the steering of quantum systems to desired objectives.
The design of optimal final-state controllers of quantum-mechanical systems that are insensitive to errors in the molecular Hamiltonian or to errors in the initial state of the system is considered. Control arises through the interaction of the system with an external field; the goal is optimal design of these latter fields for various physical objectives in the presence of system uncertainty. Sensitivity to modeling errors and other uncertainties in the molecular Hamiltonian is minimized by considering averaged costs for a family of Hamiltonian functions H0{\mathit{H}}_{0}H0​(\ensuremath{\alpha}) indexed by the random variable \ensuremath{\alpha} taking values on a compact set in Euclidean space. Similarly, sensitivity of the optimal control to the initial state is minimized by viewing the initial condition as a Hilbert-space-valued random variable and considering an optimization problem with a cost functional that is averaged over the class of initial conditions. A precise formulation of the control problem is given, and its well-posedness is established. Cost propagators are defined to display the dependence of the performance index on the initial conditions explicitly, which allows analytic averaging of initial conditions.The constrained optimization problem is reduced to an unconstrained optimization problem by the introduction of Lagrange-multiplier operators. Necessary conditions for the unconstrained problem provide the basis for a gradient search for an optimal solution. Finite-difference schemes are utilized to provide a numerical approximation of the optimal control problem. Numerical examples are given for final-state control of a diatomic molecule represented by a Morse potential illustrating design for systems with initial-phase uncertainty and parametric uncertainty. The resultant insensitive controllers execute different strategies depending on the design requirements. The controller designed to be insensitive to errors in the initial phases adopts a strategy of phase imprinting during the initial stages of the control interval to compensate for a lack of knowledge of the initial phases. It is also shown that it is not possible to coerce a system from a state with completely random initial phases to a correlated state using the class of averaged controllers considered here. The controller designed to be insensitive to parametric Hamiltonian errors adopts a strategy of amplitude restraint to prevent the wave packet from taking significant excursions into the regions where the potential is uncertain. The interesting structure exhibited by the controllers in response to the different design requirements, and the superior performance of the insensitive controllers when compared with controllers designed at nominal phases and parameters, illustrate the usefulness of the cost-averaging technique for design in the presence of uncertainties.
This work considers optimal control of quantum-mechanical systems within the framework of perturbation theory with respect to the controlling optical electric field. The control problem is expressed in terms of a cost functional including the physical objective, the penalties, and constraints. The resultant nonlinear variational equations are linearized by considering the lowest-order term in an expansion in powers of the optical-field strength. The optical field is found to satisfy a linear integral equation, and the solution may be expressed in terms of a generalized eigenvalue problem associated with the corresponding kernel. A full determination of the field is specified through the solution to the integral equation and the roots of an accompanying linearized spectral equation for a characteristic multiplier parameter. Each discrete value of the latter parameter corresponds to a particular solution to the variational equations. As a result, it is argued that under very general conditions there will be a denumerably infinite number of solutions to well-posed quantum-mechanical optimal-control problems.
We investigate the problem of control of quantum systems. Various notions of controllability of quantum systems are discussed and Lie algebraic techniques are employed to derive necessary and sufficient conditions for the resulting degrees of controllability. The results are employed to study the degree of controllability for various model systems. The problem of optimal control is formulated as a problem of maximisation of the expectation value of an observable of the system subject to dynamical constraints and costs. A class of iterative algorithms is discussed, and shown to converge to a solution of the Euler-Lagrange equations, which is a necessary condition for optimality of a solution. A numerical realisation of the algorithm is presented and potential implementation of the resulting control fields in the laboratory is discussed. The algorithm is applied to find optimal control fields for a variety of applications including selective and simultaneous excitation of individual quantum dots in a globally addressed ensemble, creation of superposition states, and selective excitation of certain vibrational modes of a simple molecule such as Hydrogen Fluoride. In each case we systematically explore the parameter space of the algorithm to select solutions with desirable physical features from the multitude of possible choices. As there are many possible optimal control fields which give similar yields the question of whether some are more robust than others is investigated. A number of optimal fields are analysed quantitatively in terms of robustness with regard to temporal and spectral noise affecting the fields, uncertainty in the model parameters and the effect of experimental limitations such as restrictions on the bandwidth of the fields. It is shown that whilst the solutions tend to be sensitive to uncertainty in the energy levels chosen for the model, they appear to be surprisingly tolerant to even severe bandwidth limitations.
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.
Recently, systematic optimal control theory has been successfully applied to design optical fields which act on a molecule to guide or transform it into a specifically prescribed state. Under realistic laboratory conditions and for real molecules, it is inevitable that there will exist disturbances and uncertainties in the control optical field, in the molecular Hamiltonian, and in the field-coupling coefficient (e.g., the dipole moment), etc. Therefore, it is important to design an optimal field which is robust to disturbances. In this paper, we provide a method through which the worst possible disturbance to the control process, under certain possible constraints, can be calculated and a corresponding optimal control field designed to be robust to the disturbances. Under disturbances of small amplitude, a perturbation solution to this problem can be applied. Illustrations are given for control of the wave-packet evolution of a diatomic molecule in the presence of disturbances.
A new formalism for the optimal control of quantum mechanical physical observables is presented. This approach is based on an analogous classical control technique reported previously. [1] Quantum Lagrange multiplier functions are used to preserve a chosen subset of the observable dynamics of interest. As a result, a corresponding small set of Lagrange multipliers needs to be calculated and they are only a function of time. This is a considerable simplification over traditional quantum optimal control theory. [2] The success of the new approach is based on taking advantage of the multiplicity of solutions to virtually any problem of quantum control to meet a physical objective. A family of such simplified formulations is introduced and numerically tested. Results are presented for these algorithms and compared with previous reported work on a model problem for selective unimolecular reaction induced by an external optical electric field. 1 I.