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

Perturbative formulation of optimal control for two-photon transitions: Reduction to an eigenvalue problem

Valery V. Dubov, HERSCHEL A. RABITZ

Published Jul 1, 1996• 9 citations
Publisher

Abstract

An alternative perturbative formulation is introduced to reduce the task of designing optimal controls of molecular motion. The method provides multiple optimal control fields as eigenvectors corresponding to the positive eigenvalues of a matrix {${\mathit{A}}_{\mathit{mn}}$}. Analytical equations for the matrix elements ${\mathit{A}}_{\mathit{mn}}$ are derived. Numerical examples are given for simple three- and four-state molecular systems. The method is shown to be computationally simple and efficient in providing good quality control results. \textcopyright{} 1996 The American Physical Society.

DOI
10.1103/physreva.54.710

Related papers

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

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

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

Optimal control of coupled quantum systems based on the first-order Magnus expansion: Application to multiple dipole-dipole-coupled molecular rotors

Andrew Ma, Alicia B Magann, Tak‐San Ho +1 more

This paper presents a method for performing approximate optimal control simulations for quantum systems with multiple coupled degrees of freedom. The time evolution is simulated using the first-order Magnus expansion in the interaction picture, where the couplings between different degrees of freedom are treated as the perturbation. A numerical implementation procedure is presented that leverages upon pairwise couplings and the separability of the zeroth-order time evolution operator to achieve a reduced computational cost, which is analyzed with respect to the number of degrees of freedom. The formulation is compatible with gradient-free methods to optimize the control field, and a stochastic hill climbing algorithm is adopted for this purpose. As illustrations, optimal control simulations are performed for systems of two and three dipole-dipole-coupled molecular rotors under the influence of a control field. For the two-rotor system, the field is optimized to achieve either orientation or entanglement objectives. For the three-rotor system, the field is optimized either to orient all three rotors in the same direction or to orient one rotor in a particular direction while the other two rotors point in the opposite direction.

arXiv:1811.12785•Physical Review A•Jul 30, 2020•8 citations

Optimal control of quantum systems: a projection approach

Chao-Jung Cheng, Chi-Chuan Hwang, Teh‐Lu Liao +1 more

This paper considers the optimal control of quantum systems. The controlled quantum systems are described by the probability-density-matrix-based Liouville–von Neumann equation. Using projection operators, the states of the quantum system are decomposed into two sub-spaces, namely the 'main state' space and the 'remaining state' space. Since the control energy is limited, a solution for optimizing the external control force is proposed in which the main state is brought to the desired main state at a certain target time, while the population of the remaining state is simultaneously suppressed in order to diminish its effects on the final population of the main state. The optimization problem is formulated by maximizing a general cost functional of states and control force. An efficient algorithm is developed to solve the optimization problem. Finally, using the hydrogen fluoride (HF) molecular population transfer problem as an illustrative example, the effectiveness of the proposed scheme for a quantum system initially in a mixed state or in a pure state is investigated through numerical simulations.

Journal of Physics A Mathematical and General•Jan 14, 2005•2 citations

Reduction of the molecular hamiltonian matrix using quantum community detection

Susan M. Mniszewski, Pavel A. Dub, Sergei Tretiak +3 more

Quantum chemistry is interested in calculating ground and excited states of molecular systems by solving the electronic Schrödinger equation. The exact numerical solution of this equation, frequently represented as an eigenvalue problem, remains unfeasible for most molecules and requires approximate methods. In this paper we introduce the use of Quantum Community Detection performed using the D-Wave quantum annealer to reduce the molecular Hamiltonian matrix in Slater determinant basis without chemical knowledge. Given a molecule represented by a matrix of Slater determinants, the connectivity between Slater determinants (as off-diagonal elements) is viewed as a graph adjacency matrix for determining multiple communities based on modularity maximization. A gauge metric based on perturbation theory is used to determine the lowest energy cluster. This cluster or sub-matrix of Slater determinants is used to calculate approximate ground state and excited state energies within chemical accuracy. The details of this method are described along with demonstrating its performance across multiple molecules of interest and bond dissociation cases. These examples provide proof-of-principle results for approximate solution of the electronic structure problem using quantum computing. This approach is general and shows potential to reduce the computational complexity of post-Hartree-Fock methods as future advances in quantum hardware become available.

Scientific Reports•Feb 18, 2021•18 citations
An alternative perturbative formulation is introduced to reduce the task of designing optimal controls of molecular motion. The method provides multiple optimal control fields as eigenvectors corresponding to the positive eigenvalues of a matrix {Amn{\mathit{A}}_{\mathit{mn}}Amn​}. Analytical equations for the matrix elements Amn{\mathit{A}}_{\mathit{mn}}Amn​ are derived. Numerical examples are given for simple three- and four-state molecular systems. The method is shown to be computationally simple and efficient in providing good quality control results. \textcopyright{} 1996 The American Physical Society.
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.
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 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.
This paper presents a method for performing approximate optimal control simulations for quantum systems with multiple coupled degrees of freedom. The time evolution is simulated using the first-order Magnus expansion in the interaction picture, where the couplings between different degrees of freedom are treated as the perturbation. A numerical implementation procedure is presented that leverages upon pairwise couplings and the separability of the zeroth-order time evolution operator to achieve a reduced computational cost, which is analyzed with respect to the number of degrees of freedom. The formulation is compatible with gradient-free methods to optimize the control field, and a stochastic hill climbing algorithm is adopted for this purpose. As illustrations, optimal control simulations are performed for systems of two and three dipole-dipole-coupled molecular rotors under the influence of a control field. For the two-rotor system, the field is optimized to achieve either orientation or entanglement objectives. For the three-rotor system, the field is optimized either to orient all three rotors in the same direction or to orient one rotor in a particular direction while the other two rotors point in the opposite direction.
This paper considers the optimal control of quantum systems. The controlled quantum systems are described by the probability-density-matrix-based Liouville–von Neumann equation. Using projection operators, the states of the quantum system are decomposed into two sub-spaces, namely the 'main state' space and the 'remaining state' space. Since the control energy is limited, a solution for optimizing the external control force is proposed in which the main state is brought to the desired main state at a certain target time, while the population of the remaining state is simultaneously suppressed in order to diminish its effects on the final population of the main state. The optimization problem is formulated by maximizing a general cost functional of states and control force. An efficient algorithm is developed to solve the optimization problem. Finally, using the hydrogen fluoride (HF) molecular population transfer problem as an illustrative example, the effectiveness of the proposed scheme for a quantum system initially in a mixed state or in a pure state is investigated through numerical simulations.
Quantum chemistry is interested in calculating ground and excited states of molecular systems by solving the electronic Schrödinger equation. The exact numerical solution of this equation, frequently represented as an eigenvalue problem, remains unfeasible for most molecules and requires approximate methods. In this paper we introduce the use of Quantum Community Detection performed using the D-Wave quantum annealer to reduce the molecular Hamiltonian matrix in Slater determinant basis without chemical knowledge. Given a molecule represented by a matrix of Slater determinants, the connectivity between Slater determinants (as off-diagonal elements) is viewed as a graph adjacency matrix for determining multiple communities based on modularity maximization. A gauge metric based on perturbation theory is used to determine the lowest energy cluster. This cluster or sub-matrix of Slater determinants is used to calculate approximate ground state and excited state energies within chemical accuracy. The details of this method are described along with demonstrating its performance across multiple molecules of interest and bond dissociation cases. These examples provide proof-of-principle results for approximate solution of the electronic structure problem using quantum computing. This approach is general and shows potential to reduce the computational complexity of post-Hartree-Fock methods as future advances in quantum hardware become available.