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

The Journal of Chemical Physics

A hybrid local/global optimal control algorithm for dissipative systems with time-dependent targets: Formulation and application to relaxing adsorbates

Stephanie Beyvers, Peter Saalfrank

Published Feb 19, 2008• 10 citations
Publisher

Abstract

Open-system quantum optimal control theory for optical control of the dynamics of a quantum system in contact with a dissipative bath is used here for explicitly time-dependent target operators, O(t). Global and local control strategies are combined in a novel algorithm by defining a set of time slices, into which the total control time is subdivided. The optimization then proceeds locally forward in time from subinterval to subinterval, while within each subinterval global control theory is used with iterative forward-backward propagation. The subintervals are connected by appropriate boundary conditions. In the present paper, all operators are represented in the basis of the eigenstates of the field-free system Hamiltonian. The algorithm is first applied to and its computational performance tested for a two-level system with energy and phase relaxation, and later extended to a many-level model. Model parameters are chosen to represent the IR pulse excitation of the adsorbate-surface stretch mode of vibrationally relaxing CO on a Cu(100) surface. Various time-dependent targets are formulated to achieve (i) population inversion, (ii) the creation of a wavepacket, and (iii) overtone excitation by "ladder climbing."

DOI
10.1063/1.2830709

Related papers

Free-time and fixed end-point optimal control theory in dissipative media: Application to entanglement generation and maintenance

K. Mishima, Koichi Yamashita

We develop monotonically convergent free-time and fixed end-point optimal control theory (OCT) in the density-matrix representation to deal with quantum systems showing dissipation. Our theory is more general and flexible for tailoring optimal laser pulses in order to control quantum dynamics with dissipation than the conventional fixed-time and fixed end-point OCT in that the optimal temporal duration of laser pulses can also be optimized exactly. To show the usefulness of our theory, it is applied to the generation and maintenance of the vibrational entanglement of carbon monoxide adsorbed on the copper (100) surface, CO/Cu(100). We demonstrate the numerical results and clarify how to combat vibrational decoherence as much as possible by the tailored shapes of the optimal laser pulses. It is expected that our theory will be general enough to be applied to a variety of dissipative quantum dynamics systems because the decoherence is one of the quantum phenomena sensitive to the temporal duration of the quantum dynamics.

The Journal of Chemical Physics•Jul 7, 2009•16 citations

Using quantum optimal control to drive intramolecular vibrational redistribution and to perform quantum computing

Ludovic Santos, Nathalie Vaeck

Quantum optimal control theory is applied to find optimal pulses for controlling the motion of an ion and a molecule for two different applications. Those optimal pulses enable the control of the dynamics of the system by driving the atom or the molecule from an initial state to desired states.The evolution equations obtained by means of the quantum optimal control theory are resolved iteratively using a monotonic convergent algorithm. A number of simulation parameters are varied in order to get the optimal pulses including the duration of the pulses, the time step of the time grid, a penalty factor that limits the maximal intensity of the fields, and a guess pulse which is used to start the optimal control.The optimal pulses obtained for each application are analyzed by Fourier transform, and also by looking at the time evolution of the populations that they generate in the system.The first application is the preparation of specific vibrational states of acetylene that are usually not reachable from the ground state, and that would remain unpopulated by usual spectroscopy. Relevant state energies and transition dipole moments are extracted from the experimental literature and especially from the global acetylene Hamiltonian conferring an uncommon precision to the control simulation. The control starts from the ground state. The target states belongs to the polyad Ns=1, Nr=5 of acetylene which includes two vibrational dark states and one vibrational bright state. First, the simulation is performed with the Schrödinger equation and in a second step, with the Liouville--von Neumann equation, as mixed states are prepared. Indeed, the control starts from a Boltzmann distribution of population in the rotational levels of the vibrational ground state chosen in order to simulate an experimental condition. But the distribution is truncated to limit the computational effort. One of the dark states appears to be a potential target for a realistic experimental investigation because the average population of the Rabi oscillation remains high and decoherence is expected to be weak. The optimal pulses obtained have a high fidelity, have a spectrum with well-resolved peak frequencies, and their experimental feasibility seems achievable within the current abilities of experimental laboratories.The second application is to propose an experimental realization of a microscopic physical device able to simulate quantum dynamics. The idea is to use the motional states of a Cd^+ ion trapped in an anharmonic potential to realize a quantum dynamics simulator of a single-particle Schrödinger equation. In this way, the motional states store the information and the optimal pulse manipulates this information to realize operations. In the present case, the simulated dynamics was the propagation of a wave packet in a harmonic potential. Starting from an initial quantum state, the pulse acts on the system to modify the motional states of the ion in such a way that the final superposition of motional states corresponds to the results of the dynamics. This simulation is performed with the Liouville--von Neumann equation and also with the Lindblad equation as dissipation is included to test the robustness of the pulse against perturbations of the potential. The optimal pulses that are obtained have a high fidelity which shows that the ion trap system has correctly realized the quantum dynamics simulation. The optimal pulses are valid for any initial condition if the potential of the simulation or the mass of the propagated wave packet is unchanged.

•Nov 28, 2017•0 citations

OPTIMALLY CONTROLLED VIBRATIONAL POPULATION TRANSFER IN A DIATOMIC QUANTUM SYSTEM

Praveen Kumar, Sitansh Sharma, Harjinder Singh

A time-dependent formulation of quantum control is employed to investigate optimally controlled vibrational population transfer in a diatomic quantum system. The problem of finding the optimal laser field needed to achieve a specific quantum transition from an initial state to the desired target goal is formulated using an iterative method and the conjugate gradient method (CGM). The time-dependent Schrödinger equation is solved with interaction of laser radiation with matter included within a dipole approximation in the Hamiltonian. Appropriate boundary conditions are chosen for the evolution problem. The control objective is chosen as the value of transition probability from an initial state to a target state. A comparison is made between the results obtained using the iterative method and the CGM for optimization. Finally, quantum bits are encoded using the vibrational states of the diatomic in the regime of low-vibrational excitation.

Journal of Theoretical and Computational Chemistry•Feb 1, 2009•11 citations

Dynamics of high-dimensional quantum systems coupled to a harmonic bath. General theory and implementation via multiconfigurational wave packets and truncated hierarchical equations for the mean-fields

David Picconi

Modeling the dynamics of a quantum system coupled to a dissipative environment becomes particularly challenging when the system's dimensionality is too high to permit the computation of its eigenstates. This problem is addressed by introducing an eigenstate-free formalism, where the open quantum system is represented as a mixture of high-dimensional, time-dependent wave packets governed by coupled Schrödinger equations, while the environment is described by a multi-component quantum master equation. An efficient computational implementation of this formalism is presented, employing a variational mixed Gaussian/multiconfigurational time-dependent Hartree (G-MCTDH) ansatz for the wave packets and propagating the environment dynamics via hierarchical equations, truncated at the first or second level of the hierarchy. The effectiveness of the proposed methodology is demonstrated on a 61-dimensional model of phonon-driven vibrational relaxation of an adsorbate. G-MCTDH calculations on 4- and 10-dimensional reduced models, combined with truncated hierarchical equations for the mean fields, nearly quantitatively replicate the full-dimensional quantum dynamical results on vibrational relaxation while significantly reducing the computational time. This approach thus offers a promising quantum dynamical method for modeling complex system-bath interactions, where a large number of degrees of freedom must be explicitly considered.

The Journal of Chemical Physics•Oct 25, 2024•7 citations

Near-time-optimal control for quantum systems

Qiming Chen, Re-Bing Wu, Tian-Ming Zhang +1 more

For a quantum system controlled by an external field, time-optimal control is referred to as the shortest-time-duration control that can still permit maximizing an objective function $J$, which is especially a desirable goal for engineering quantum dynamics against decoherence effects. However, since rigorously finding a time-optimal control is usually very difficult and in many circumstances the control is only required to be sufficiently short and precise, one can design algorithms seeking such suboptimal control solutions for much reduced computational effort. In this paper, we propose an iterative algorithm for finding near-time-optimal control in a high level set (i.e., the set of controls that achieves the same value of $J$) that can be arbitrarily close to the global optima. The algorithm proceeds seeking to decrease the time duration $T$ while the value of $J$ remains invariant, until $J$ leaves the level-set value; the deviation of $J$ due to numerical errors is corrected by gradient climbing that brings the search back to the level-set $J$ value. Since the level set is very close to the maximum value of $J$, the resulting control solution is nearly time optimal with manageable precision. Numerical examples demonstrate the effectiveness and general applicability of the algorithm.

Physical Review A•Dec 16, 2015•20 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
Open-system quantum optimal control theory for optical control of the dynamics of a quantum system in contact with a dissipative bath is used here for explicitly time-dependent target operators, O(t). Global and local control strategies are combined in a novel algorithm by defining a set of time slices, into which the total control time is subdivided. The optimization then proceeds locally forward in time from subinterval to subinterval, while within each subinterval global control theory is used with iterative forward-backward propagation. The subintervals are connected by appropriate boundary conditions. In the present paper, all operators are represented in the basis of the eigenstates of the field-free system Hamiltonian. The algorithm is first applied to and its computational performance tested for a two-level system with energy and phase relaxation, and later extended to a many-level model. Model parameters are chosen to represent the IR pulse excitation of the adsorbate-surface stretch mode of vibrationally relaxing CO on a Cu(100) surface. Various time-dependent targets are formulated to achieve (i) population inversion, (ii) the creation of a wavepacket, and (iii) overtone excitation by "ladder climbing."
We develop monotonically convergent free-time and fixed end-point optimal control theory (OCT) in the density-matrix representation to deal with quantum systems showing dissipation. Our theory is more general and flexible for tailoring optimal laser pulses in order to control quantum dynamics with dissipation than the conventional fixed-time and fixed end-point OCT in that the optimal temporal duration of laser pulses can also be optimized exactly. To show the usefulness of our theory, it is applied to the generation and maintenance of the vibrational entanglement of carbon monoxide adsorbed on the copper (100) surface, CO/Cu(100). We demonstrate the numerical results and clarify how to combat vibrational decoherence as much as possible by the tailored shapes of the optimal laser pulses. It is expected that our theory will be general enough to be applied to a variety of dissipative quantum dynamics systems because the decoherence is one of the quantum phenomena sensitive to the temporal duration of the quantum dynamics.
Quantum optimal control theory is applied to find optimal pulses for controlling the motion of an ion and a molecule for two different applications. Those optimal pulses enable the control of the dynamics of the system by driving the atom or the molecule from an initial state to desired states.The evolution equations obtained by means of the quantum optimal control theory are resolved iteratively using a monotonic convergent algorithm. A number of simulation parameters are varied in order to get the optimal pulses including the duration of the pulses, the time step of the time grid, a penalty factor that limits the maximal intensity of the fields, and a guess pulse which is used to start the optimal control.The optimal pulses obtained for each application are analyzed by Fourier transform, and also by looking at the time evolution of the populations that they generate in the system.The first application is the preparation of specific vibrational states of acetylene that are usually not reachable from the ground state, and that would remain unpopulated by usual spectroscopy. Relevant state energies and transition dipole moments are extracted from the experimental literature and especially from the global acetylene Hamiltonian conferring an uncommon precision to the control simulation. The control starts from the ground state. The target states belongs to the polyad Ns=1, Nr=5 of acetylene which includes two vibrational dark states and one vibrational bright state. First, the simulation is performed with the Schrödinger equation and in a second step, with the Liouville--von Neumann equation, as mixed states are prepared. Indeed, the control starts from a Boltzmann distribution of population in the rotational levels of the vibrational ground state chosen in order to simulate an experimental condition. But the distribution is truncated to limit the computational effort. One of the dark states appears to be a potential target for a realistic experimental investigation because the average population of the Rabi oscillation remains high and decoherence is expected to be weak. The optimal pulses obtained have a high fidelity, have a spectrum with well-resolved peak frequencies, and their experimental feasibility seems achievable within the current abilities of experimental laboratories.The second application is to propose an experimental realization of a microscopic physical device able to simulate quantum dynamics. The idea is to use the motional states of a Cd^+ ion trapped in an anharmonic potential to realize a quantum dynamics simulator of a single-particle Schrödinger equation. In this way, the motional states store the information and the optimal pulse manipulates this information to realize operations. In the present case, the simulated dynamics was the propagation of a wave packet in a harmonic potential. Starting from an initial quantum state, the pulse acts on the system to modify the motional states of the ion in such a way that the final superposition of motional states corresponds to the results of the dynamics. This simulation is performed with the Liouville--von Neumann equation and also with the Lindblad equation as dissipation is included to test the robustness of the pulse against perturbations of the potential. The optimal pulses that are obtained have a high fidelity which shows that the ion trap system has correctly realized the quantum dynamics simulation. The optimal pulses are valid for any initial condition if the potential of the simulation or the mass of the propagated wave packet is unchanged.
A time-dependent formulation of quantum control is employed to investigate optimally controlled vibrational population transfer in a diatomic quantum system. The problem of finding the optimal laser field needed to achieve a specific quantum transition from an initial state to the desired target goal is formulated using an iterative method and the conjugate gradient method (CGM). The time-dependent Schrödinger equation is solved with interaction of laser radiation with matter included within a dipole approximation in the Hamiltonian. Appropriate boundary conditions are chosen for the evolution problem. The control objective is chosen as the value of transition probability from an initial state to a target state. A comparison is made between the results obtained using the iterative method and the CGM for optimization. Finally, quantum bits are encoded using the vibrational states of the diatomic in the regime of low-vibrational excitation.
Modeling the dynamics of a quantum system coupled to a dissipative environment becomes particularly challenging when the system's dimensionality is too high to permit the computation of its eigenstates. This problem is addressed by introducing an eigenstate-free formalism, where the open quantum system is represented as a mixture of high-dimensional, time-dependent wave packets governed by coupled Schrödinger equations, while the environment is described by a multi-component quantum master equation. An efficient computational implementation of this formalism is presented, employing a variational mixed Gaussian/multiconfigurational time-dependent Hartree (G-MCTDH) ansatz for the wave packets and propagating the environment dynamics via hierarchical equations, truncated at the first or second level of the hierarchy. The effectiveness of the proposed methodology is demonstrated on a 61-dimensional model of phonon-driven vibrational relaxation of an adsorbate. G-MCTDH calculations on 4- and 10-dimensional reduced models, combined with truncated hierarchical equations for the mean fields, nearly quantitatively replicate the full-dimensional quantum dynamical results on vibrational relaxation while significantly reducing the computational time. This approach thus offers a promising quantum dynamical method for modeling complex system-bath interactions, where a large number of degrees of freedom must be explicitly considered.
For a quantum system controlled by an external field, time-optimal control is referred to as the shortest-time-duration control that can still permit maximizing an objective function JJJ, which is especially a desirable goal for engineering quantum dynamics against decoherence effects. However, since rigorously finding a time-optimal control is usually very difficult and in many circumstances the control is only required to be sufficiently short and precise, one can design algorithms seeking such suboptimal control solutions for much reduced computational effort. In this paper, we propose an iterative algorithm for finding near-time-optimal control in a high level set (i.e., the set of controls that achieves the same value of JJJ) that can be arbitrarily close to the global optima. The algorithm proceeds seeking to decrease the time duration TTT while the value of JJJ remains invariant, until JJJ leaves the level-set value; the deviation of JJJ due to numerical errors is corrected by gradient climbing that brings the search back to the level-set JJJ value. Since the level set is very close to the maximum value of JJJ, the resulting control solution is nearly time optimal with manageable precision. Numerical examples demonstrate the effectiveness and general applicability of the algorithm.
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.