We propose some iterative algorithms for the generation of open-loop controls of quantum systems. The task that these controls are intended for is the transfer of a given initial state to a given final state with maximal probability. We prove that the algorithms are converging and the resulting controls are optimal.
We develop a family of monotonically convergent iterative algorithms for solving a wide class of optimal control problems in which a dynamical system interacts nonlinearly with a control. The key idea is to divide a control into identical components and to introduce auxiliary steps to update each component at every iteration step. The algorithms are proved to exhibit monotonic convergence, which is also numerically confirmed through a case study of the control of molecular orientation. The numerical results show that high-quality solutions can be obtained by using the present algorithms.
The optimal control problem for open quantum systems can be formulated as a time-dependent Lindbladian that is parameterized by a number of time-dependent control variables. Given an observable and an initial state, the goal is to tune the control variables so that the expected value of some observable with respect to the final state is maximized. In this paper, we present algorithms for solving this optimal control problem efficiently, i.e., having a poly-logarithmic dependency on the system dimension, which is exponentially faster than best-known classical algorithms. Our algorithms are hybrid, consisting of both quantum and classical components. The quantum procedure simulates time-dependent Lindblad evolution that drives the initial state to the final state, and it also provides access to the gradients of the objective function via quantum gradient estimation. The classical procedure uses the gradient information to update the control variables. At the technical level, we provide the first (to the best of our knowledge) simulation algorithm for time-dependent Lindbladians with an ℓ1-norm dependence. As an alternative, we also present a simulation algorithm in the interaction picture to improve the algorithm for the cases where the time-independent component of a Lindbladian dominates the time-dependent part. On the classical side, we heavily adapt the state-of-the-art classical optimization analysis to interface with the quantum part of our algorithms. Both the quantum simulation techniques and the classical optimization analyses might be of independent interest.
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 optimal control (QOC) provides a systematic framework for achieving high-fidelity operations in quantum systems and plays a central role in tasks such as gate synthesis, state transfer, and pulse design. Existing QOC methods broadly fall into two categories: gradient-based and gradient-free algorithms. The associated optimization landscape is often high-dimensional, non-convex, and populated by numerous local minima, making efficient gradient-free search strategies essential. To address this, we introduce a gradient-free matrix product state/tensor train (MPS/TT) sampling heuristic for discrete quantum optimal control. In our approach, the MPS defines a score function over the space of discrete control parameters, which in turn induces a sampling distribution over candidate control sequences. This distribution is iteratively refined through selection of better performing sequences and local tensor updates to bias the search toward high-performing sequences. We evaluate the method on a range of benchmark problems, including single-qubit state transfer, Bell-pair preparation, qutrit gate implementation, and open-system population transfer. Across these tasks, the method exhibits stable convergence behavior and competitive empirical performance relative to established gradient-free baselines. These results suggest that tensor network sampling offers a viable heuristic framework for discrete quantum control.
The control of quantum dynamics is often concerned with finding time-dependent optimal control fields that can take a system from an initial state to a final state to attain the desired value of an observable. This paper presents a general method for formulating monotonically convergent algorithms to iteratively improve control fields. The formulation is based on a two-point boundary-value quantum control paradigm (TBQCP) expressed as a nonlinear integral equation of the first kind arising from dynamical invariant tracking control. TBQCP is shown to be related to various existing techniques, including local control theory, the Krotov method, and optimal control theory. Several accelerated monotonic convergence schemes for iteratively computing control fields are derived based on TBQCP. Numerical simulations are compared with the Krotov method showing that the new TBQCP schemes are efficient and remain monotonically convergent over a wide range of the iteration step parameters and the control pulse lengths, which is attributable to the trap-free character of the transition probability quantum dynamics control landscape.