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

Journal of Chemical Theory and Computation

Classical Optimal Control for Energy Minimization Based On Diffeomorphic Modulation under Observable-Response-Preserving Homotopy

Micheline B. Soley, Andreas Markmann, Víctor S. Batista

Published Apr 20, 2018• 4 citations
Publisher

Abstract

We introduce the so-called "Classical Optimal Control Optimization" (COCO) method for global energy minimization based on the implementation of the diffeomorphic modulation under observable-response-preserving homotopy (DMORPH) gradient algorithm. A probe particle with time-dependent mass m( t;β) and dipole μ( r, t;β) is evolved classically on the potential energy surface V( r) coupled to an electric field E( t;β), as described by the time-dependent density of states represented on a grid, or otherwise as a linear combination of Gaussians generated by the k-means clustering algorithm. Control parameters β defining m( t;β), μ( r, t;β), and E( t;β) are optimized by following the gradients of the energy with respect to β, adapting them to steer the particle toward the global minimum energy configuration. We find that the resulting COCO algorithm is capable of resolving near-degenerate states separated by large energy barriers and successfully locates the global minima of golf potentials on flat and rugged surfaces, previously explored for testing quantum annealing methodologies and the quantum optimal control optimization (QuOCO) method. Preliminary results show successful energy minimization of multidimensional Lennard-Jones clusters. Beyond the analysis of energy minimization in the specific model systems investigated, we anticipate COCO should be valuable for solving minimization problems in general, including optimization of parameters in applications to machine learning and molecular structure determination.

DOI
10.1021/acs.jctc.8b00124

Related papers

Steered Quantum Dynamics for Energy Minimization

Micheline B. Soley, Andreas Markmann, Víctor S. Batista

We introduce a quantum optimal control algorithm for energy minimization that combines the diffeomorphic modulation under observable response preserving homotopy (D-MORPH) gradient and the Broyden Fletcher Goldfarb Shanno (BFGS) iterative scheme for nonlinear optimization. An extended set of controls defining the time-dependent mass, dipole moment, and external perturbational field are optimized to find an effective Hamiltonian that steers the dynamics of the system into the global minimum without getting trapped into local minima. The algorithm is illustrated as applied to energy minimization on rugged surfaces and golf potentials comparable to those previously explored for testing quantum annealing methodologies.

The Journal of Physical Chemistry B•Aug 3, 2014•5 citations

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

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

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.

Physical Review A•Jun 1, 1988•890 citations

Observable-preserving control of quantum dynamics over a family of related systems

Adam E. Rothman, Tak-San Ho, HERSCHEL A. RABITZ

Quantum control aims at the manipulation of atomic- and molecular-scale dynamics phenomena. An important objective in this regard is the understanding of dynamical control within a family of related quantum systems. To explore this issue, diffeomorphic changes in the system Hamiltonian $\mathcal{H}(s,t)$ are introduced by scanning over a homotopy parameter $s$ and then monitoring the control field response needed to maintain the value of a specified target observable. This operation is implemented through a procedure referred to as diffeomorphic modulation under observable-response-preserving homotopy (D-MORPH). The governing D-MORPH differential equation determining the control laser field $E(s,t)$ is shown to explicitly allow for innumerable solutions, with each characterized by the choice of an arbitrary function $f(s,t)$ of $s$ and time $t$. The presence of $f(s,t)$ in the D-MORPH differential equation makes clear the origin of multiple control fields that produce the same observable objective. A stable algorithm is presented for practical execution of D-MORPH with the only criterion that the Hamiltonian $\mathcal{H}(s,t)$ permit reaching the objective over the full domain of $s$ being sampled. Both analytic and numerical examples are presented to illustrate the D-MORPH concept.

Physical Review A•Aug 25, 2005•85 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

Global Optimization by Quantum Dynamics with Energy Dissipation

Teturo Itami

Both reduction in computational load and setting parameters with ease are necessary in calculation of global optimization. For this we propose a new quantum mechanical method of the problem. Global optimization is achieved by a quantum particle with energy dissipation due to frictional force. Canonical quantization is applied to a system that consists of one mass point under friction proportional to its velocity. Realistic motion of a particle is abstracted from a resultant quantum system according to a causal interpretation. By simulation of the quantum particle motion under potential functions V(x) with local minima and one global minimum, it is clarified that the particle with its initial position at the local minimum at rest necessarily arrives at the global minimum. Only one particle is needed in the new algorithm. When we clarify mathematical structure of our algorithm, difficulties found in quantum annealing or conventional method for global optimization will thereby be solved in the new method.

IEEJ Transactions on Electronics Information and Systems•Jan 1, 2013•0 citations

Second-Order Adjoint Method for Quantum Optimal Control

Harish S. Bhat

We derive and implement a second-order adjoint method to compute exact gradients and Hessians for a prototypical quantum optimal control problem, that of solving for the minimal energy applied electric field that drives a molecule from a given initial state to a desired target state. For small to moderately sized systems, we demonstrate a vectorized GPU implementation of a second-order adjoint method that computes both Hessians and gradients with wall times only marginally more than those required to compute gradients via commonly used first-order adjoint methods. Pairing our secondorder adjoint method with a trust region optimizer (a type of Newton method), we show that it outperforms a first-order method, requiring significantly fewer iterations and wall time to find optimal controls for four molecular systems.

IEEE International Conference on Quantum Control, Computing and Learning (qCCL)•Jun 25, 2025•0 citations
We introduce the so-called "Classical Optimal Control Optimization" (COCO) method for global energy minimization based on the implementation of the diffeomorphic modulation under observable-response-preserving homotopy (DMORPH) gradient algorithm. A probe particle with time-dependent mass m( t;β) and dipole μ( r, t;β) is evolved classically on the potential energy surface V( r) coupled to an electric field E( t;β), as described by the time-dependent density of states represented on a grid, or otherwise as a linear combination of Gaussians generated by the k-means clustering algorithm. Control parameters β defining m( t;β), μ( r, t;β), and E( t;β) are optimized by following the gradients of the energy with respect to β, adapting them to steer the particle toward the global minimum energy configuration. We find that the resulting COCO algorithm is capable of resolving near-degenerate states separated by large energy barriers and successfully locates the global minima of golf potentials on flat and rugged surfaces, previously explored for testing quantum annealing methodologies and the quantum optimal control optimization (QuOCO) method. Preliminary results show successful energy minimization of multidimensional Lennard-Jones clusters. Beyond the analysis of energy minimization in the specific model systems investigated, we anticipate COCO should be valuable for solving minimization problems in general, including optimization of parameters in applications to machine learning and molecular structure determination.
We introduce a quantum optimal control algorithm for energy minimization that combines the diffeomorphic modulation under observable response preserving homotopy (D-MORPH) gradient and the Broyden Fletcher Goldfarb Shanno (BFGS) iterative scheme for nonlinear optimization. An extended set of controls defining the time-dependent mass, dipole moment, and external perturbational field are optimized to find an effective Hamiltonian that steers the dynamics of the system into the global minimum without getting trapped into local minima. The algorithm is illustrated as applied to energy minimization on rugged surfaces and golf potentials comparable to those previously explored for testing quantum annealing methodologies.
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.
Quantum control aims at the manipulation of atomic- and molecular-scale dynamics phenomena. An important objective in this regard is the understanding of dynamical control within a family of related quantum systems. To explore this issue, diffeomorphic changes in the system Hamiltonian H(s,t)\mathcal{H}(s,t)H(s,t) are introduced by scanning over a homotopy parameter sss and then monitoring the control field response needed to maintain the value of a specified target observable. This operation is implemented through a procedure referred to as diffeomorphic modulation under observable-response-preserving homotopy (D-MORPH). The governing D-MORPH differential equation determining the control laser field E(s,t)E(s,t)E(s,t) is shown to explicitly allow for innumerable solutions, with each characterized by the choice of an arbitrary function f(s,t)f(s,t)f(s,t) of sss and time ttt. The presence of f(s,t)f(s,t)f(s,t) in the D-MORPH differential equation makes clear the origin of multiple control fields that produce the same observable objective. A stable algorithm is presented for practical execution of D-MORPH with the only criterion that the Hamiltonian H(s,t)\mathcal{H}(s,t)H(s,t) permit reaching the objective over the full domain of sss being sampled. Both analytic and numerical examples are presented to illustrate the D-MORPH concept.
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.
Both reduction in computational load and setting parameters with ease are necessary in calculation of global optimization. For this we propose a new quantum mechanical method of the problem. Global optimization is achieved by a quantum particle with energy dissipation due to frictional force. Canonical quantization is applied to a system that consists of one mass point under friction proportional to its velocity. Realistic motion of a particle is abstracted from a resultant quantum system according to a causal interpretation. By simulation of the quantum particle motion under potential functions V(x) with local minima and one global minimum, it is clarified that the particle with its initial position at the local minimum at rest necessarily arrives at the global minimum. Only one particle is needed in the new algorithm. When we clarify mathematical structure of our algorithm, difficulties found in quantum annealing or conventional method for global optimization will thereby be solved in the new method.
We derive and implement a second-order adjoint method to compute exact gradients and Hessians for a prototypical quantum optimal control problem, that of solving for the minimal energy applied electric field that drives a molecule from a given initial state to a desired target state. For small to moderately sized systems, we demonstrate a vectorized GPU implementation of a second-order adjoint method that computes both Hessians and gradients with wall times only marginally more than those required to compute gradients via commonly used first-order adjoint methods. Pairing our secondorder adjoint method with a trust region optimizer (a type of Newton method), we show that it outperforms a first-order method, requiring significantly fewer iterations and wall time to find optimal controls for four molecular systems.