This paper presents a framework for solving the pure-state preparation problem using numerical optimal control. As an example, we consider the case where a number of qubits are dispersively coupled to a readout cavity. We model open system quantum dynamics using the Markovian Lindblad master equation driven by external control pulses. The main result of this paper develops a basis of density matrices (a parameterization) where each basis element is a density matrix by itself. Utilizing a specific objective function, we show how an ensemble of the basis elements can be used as a single initial state throughout the optimization process—independent of the system dimension. We apply the general framework to the specific application of ground-state reset of one and two qubits coupled to a readout cavity.
This paper presents a framework for solving the pure-state preparation problem using numerical optimal control. As an example, we consider the case where a number of qubits are dispersively coupled to a readout cavity. We model open system quantum dynamics using the Markovian Lindblad master equation driven by external control pulses. The main result of this paper develops a basis of density matrices (a parameterization) where each basis element is a density matrix by itself. Utilizing a specific objective function, we show how an ensemble of the basis elements can be used as a single initial state throughout the optimization process—independent of the system dimension. We apply the general framework to the specific application of ground-state reset of one and two qubits coupled to a readout cavity.
This paper presents a numerical optimization framework for solving the ground state reset problem where a number of qubits are dispersively coupled to a readout cavity. We model the open system quantum dynamics using Lindblad's master equation, driven by external control pulses. We develop a basis of density matrices (a parameterization) where each basis element is a density matrix itself. Together with a new objective function, we show how a superposition of the basis elements can be used in such a way that the objective function can be evaluated by solving the master equation for one initial condition only - independent of the system dimension. This enables efficient ground state reset allowing for an increasing number of qubits. Numerical results demonstrate efficient optimal control for ground-state reset of one and two qubits coupled to a readout cavity.
We present a framework that combines the adjoint-state method together with reverse-time backpropagation to solve prohibitively large open-system quantum control problems. Our approach enables the optimization of arbitrary cost functions with fully general controls applied on large open quantum systems described by a Lindblad master equation. It is scalable, computationally efficient, and has a low-memory footprint. We apply this framework to optimize two inherently dissipative operations in superconducting qubits which lag behind in terms of fidelity and duration compared to other unitary operations: the dispersive readout and all-microwave reset of a transmon qubit. Our results show that while standard pulses for dispersive readout are nearly optimal, adding a transmon drive during the protocol can yield 2× improvements in fidelity and duration. We further demonstrate a 2× improvement in reset fidelity and duration through pulse shaping, indicating significant potential for enhancement in reset protocols. Our approach can readily be applied to optimize quantum controls in a vast range of applications such as reservoir engineering, autonomous quantum error correction, and leakage-reduction units.
This article provides a review of recent developments in the formulation and execution of optimal control strategies for the dynamics of quantum systems. A brief introduction to the concept of optimal control, the dynamics of of open quantum systems, and quantum information processing is followed by a presentation of recent developments regarding the two main tasks in this context: state-specific and state-independent optimal control. For the former, we present an extension of conventional theory (Pontryagin's principle) to quantum systems which undergo a non-Markovian time-evolution. Owing to its importance for the realization of quantum information processing, the main body of the review, however, is devoted to state-independent optimal control. Here, we address three different approaches: an approach which treats dissipative effects from the environment in lowest-order perturbation theory, a general method based on the time--evolution superoperator concept, as well as one based on the Kraus representation of the time-evolution superoperator. Applications which illustrate these new methods focus on single and double qubits (quantum gates) whereby the environment is modeled either within the Lindblad equation or a bath of bosons (spin-boson model). While these approaches are widely applicable, we shall focus our attention to solid-state based physical realizations, such as semiconductor- and superconductor-based systems. While an attempt is made to reference relevant and representative work throughout the community, the exposition will focus mainly on work which has emerged from our own group.
We demonstrate that the control protocols of quantum information devices can be simulated by assuming a low-rank ansatz for the density matrix. The rationale underlying this assumption is that quantum information protocols, by design, operate in a regime of nearly pure quantum states. Within the low-rank assumption, the simulation of these protocols is considerably faster than solving the full Lindblad master equation. This advantage can be used to increase the accuracy of the simulation by avoiding uncontrolled approximations and to streamline protocol optimization. We benchmark our approach on the optimization of the transmon qubit dispersive readout in a realistic transmon-resonator-filter model. With a Hilbert space dimension N = 2000 , assuming a rank as low as M = 20 , we achieve a nearly 100-fold speedup compared to full master equation integration while accurately reproducing all relevant observables. By combining the low-rank approximation with a compact pulse parametrization and gradient-free optimization, we obtain state-of-the-art readout assignment errors ɛ a ≈ 1.2 × 10 − 3 for a 40 ns readout pulse schedule, while comfortably running on a laptop and not relying on the rotating-wave approximation. Our approach is broadly applicable to most quantum control protocols, including quantum gates, state preparation, and fast reset operations. This establishes low-rank methods as a general tool for optimal control across diverse quantum platforms.
Manipulation of a quantum system requires the knowledge of how it evolves. To impose that the dynamics of a system becomes a particular target operation (for any preparation of the system), it may be more useful to have an equation of motion for the dynamics itself, rather than the state. Here we develop a Markovian master equation for the process matrix of an open system interacting with a reservoir, which resembles the Lindblad Markovian master equation. We employ this equation to introduce a scheme for optimal local coherent process control at target times and extend the Krotov technique to obtain optimal control. We illustrate utility of this framework through several quantum coherent-control scenarios, such as optimal decoherence suppression, gate simulation, and passive and indirect control of the reservoir, in all of which we aim to simulate a given terminal process at a given final time.
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.