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

arXiv:2607.14217

Analytic gradients for low-rank quantum optimal control

Leo Goutte, Vincenzo Savona

quant-phSubmitted Jul 15, 2026• Announced Jul 17, 2026• 0 citations
arXivHTMLPDF

Abstract

We introduce low-rank optimal control (LROC), a method for designing control pulses in open quantum systems whose full density-matrix simulation is prohibitively expensive. The method exploits a feature of quantum computing itself: because protocols are designed to preserve purity, the density matrix is dominated by a few pure states and admits an accurate low-rank factorization. LROC propagates only this factorized form and, by deriving the corresponding adjoint equation, obtains the gradient of any differentiable objective at the same reduced cost as the simulation, leading to a quadratic improvement in time and memory compared to the full master equation. We illustrate the breadth of the method on four superconducting-circuit tasks: preparation of a five-qubit GHZ state, a CNOT gate, qubit readout, and an error correction primitive, modeled with realistic multilevel transmons, decay, and strong drives, in each case reaching fidelities consistent with the intrinsic dissipation limits. LROC thereby extends pulse-level optimization to system sizes beyond the reach of existing gradient-based methods.

Comments
16 pages, 6 figures

Related papers

Low-rank optimal control of quantum devices

Leo Goutte, Vincenzo Savona

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.

arXiv:2508.18114•Physical Review Research•Dec 19, 2025•0 citations

Gradient-based optimal control of open quantum systems using quantum trajectories and automatic differentiation

Mohamed R. Abdelhafez, David Schuster, Jens Koch

We present a gradient-based optimal-control technique for open quantum systems that utilizes quantum trajectories to simulate the quantum dynamics during optimization. Using trajectories allows for optimizing open systems with less computational cost than the regular density matrix approaches in most realistic optimization problems. We introduce an improved-sampling algorithm which minimizes the number of trajectories needed per optimization iteration. Together with employing stochastic gradient descent techniques, this reduces the complexity of optimizing many realistic open quantum systems to the complexity encountered with closed systems. Our optimizer harnesses automatic differentiation to provide flexibility in optimization and to suit the different constraints and diverse parameter regimes of real-life experiments. We utilize the optimizer in a variety of applications to demonstrate how the use of quantum trajectories significantly reduces the computation complexity while achieving a multitude of simultaneous optimization targets. Demonstrated targets include high state-transfer fidelities despite dissipation, faster gate times, and maximization of qubit-readout fidelity while maintaining the quantum nondemolition nature of the measurement and allowing for subsequent fast resonator reset.

arXiv:1901.05541•Physical Review A•May 20, 2019•86 citations

Gradient projection method for constrained quantum control

Oleg Vasilevich Morzhin, Alexander Nikolaevich Pechen

Abstract In this work, we adopt the gradient projection method (GPM) to problems of quantum control. For general N-level closed and open quantum systems, we derive the corresponding adjoint systems and gradients of the objective functionals and provide the projection versions of the Pontryagin maximum principle and the GPM, all directly in terms of quantum objects such as evolution operator, Hamiltonians, density matrices, etc. Various forms of the GPM, including one- and two-step, are provided and compared. We formulate the GPM both for closed and open quantum systems, latter for the general case with simultaneous coherent and incoherent controls. The GPM is designed to perform local gradient based optimization in the case when bounds are imposed on the controls. The main advantage of the method is that it allows to exactly satisfy the bounds, in difference to other approaches such as adding constraints as weight to an objective. We apply the GPM to several examples including generation of one- and two-qubit gates and two-qubit Bell and Werner states for models of superconducting qubits under the constraint when controls are zero at the initial and final times, steering an open quantum system to a target density matrix for simulating action of the Werner–Holevo channel, etc.

arXiv:2411.19644•Journal of Physics A Mathematical and Theoretical•Mar 27, 2025•7 citations

Quantum Optimal Control via Semi-Automatic Differentiation

Michael H. Goerz, Sebastian Camilo Carrasco, Vladimir S. Malinovsky

We develop a framework of "semi-automatic differentiation" that combines existing gradient-based methods of quantum optimal control with automatic differentiation. The approach allows to optimize practically any computable functional and is implemented in two open source Julia packages, GRAPE.jl and Krotov.jl , part of the QuantumControl.jl framework. Our method is based on formally rewriting the optimization functional in terms of propagated states, overlaps with target states, or quantum gates. An analytical application of the chain rule then allows to separate the time propagation and the evaluation of the functional when calculating the gradient. The former can be evaluated with great efficiency via a modified GRAPE scheme. The latter is evaluated with automatic differentiation, but with a profoundly reduced complexity compared to the time propagation. Thus, our approach eliminates the prohibitive memory and runtime overhead normally associated with automatic differentiation and facilitates further advancement in quantum control by enabling the direct optimization of non-analytic functionals for quantum information and quantum metrology, especially in open quantum systems. We illustrate and benchmark the use of semi-automatic differentiation for the optimization of perfectly entangling quantum gates on superconducting qubits coupled via a shared transmission line. This includes the first direct optimization of the non-analytic gate concurrence.

arXiv:2205.15044•Quantum•Dec 7, 2022•33 citations

Structure-preserving full- and low-rank exponential methods for the optimal control of Lindblad equations

Hao Chen, Alfio Borzı̀

Optimal control of open quantum systems governed by the Lindblad master equation requires the repeated solution of forward and adjoint evolution equations. For large Hilbert spaces, these computations demand numerical methods that are structure-preserving, accurate, and computationally efficient, while remaining compatible with nonsmooth optimization techniques. In this work, a unified framework for optimal control of Lindblad equations is presented. The continuous optimality system is derived from the Pontryagin maximum principle and solved by a sequential quadratic Hamiltonian (SQH) method. Its numerical realization is based on second-order exponential midpoint propagators for the forward and adjoint Lindblad equations with time-dependent Hamiltonians. Full-rank schemes preserving the Hermitian and positive-semidefinite structure of the density matrix are developed together with low-rank formulations that substantially reduce storage requirements and computational cost. Rigorous error estimates are established for the full- and low-rank forward and adjoint propagators. Numerical experiments confirm the predicted convergence rates, demonstrate the effectiveness of the low-rank approximations, and illustrate the performance of the proposed FREM-SQH and LREM-SQH algorithms for optimal control problems with smooth and nonsmooth control costs.

arXiv:2610.09403•Oct 7, 2026•0 citations

Analytical Series Expansion for Efficient Gradient Evaluation in Multi-Qubit Optimal Control

Ashutosh Mishra, Elena Lupo, Frank K. Wilhelm +1 more

Gradient-based quantum optimal control and the theory of operator evolution are rarely discussed together. In this Letter, we bridge this gap by showing that the gradient of the time-evolved propagator involves the Heisenberg evolution of local operators. We present a unifying framework by deriving from first principles the formal solution for the gradient under an arbitrary pulse parameterization, and show that the same construction extends to derivatives of any order. We obtain a series of nested, time-independent commutators weighted by time-dependent scalar coefficients. The commutators are thus computed once, and each optimization step only updates the scalars. The method is particularly suited for simulating optimal control tasks in quantum systems with local interactions, which is a common situation in large multi-qubit platforms. In this setting, light-cone arguments and sparse-Pauli truncation heuristics may be used to keep the number of relevant terms in the series small. We compare the computational cost required for the series with the Gradient Optimization of Analytic conTrols (GOAT) method, and, focusing on the problem of preparation of a GHZ state, demonstrate more than an order of magnitude speedup for a qubit ladder and a chain geometry.

quant-pharXiv:2607.26867•Jul 30, 2026•0 citations
We introduce low-rank optimal control (LROC), a method for designing control pulses in open quantum systems whose full density-matrix simulation is prohibitively expensive. The method exploits a feature of quantum computing itself: because protocols are designed to preserve purity, the density matrix is dominated by a few pure states and admits an accurate low-rank factorization. LROC propagates only this factorized form and, by deriving the corresponding adjoint equation, obtains the gradient of any differentiable objective at the same reduced cost as the simulation, leading to a quadratic improvement in time and memory compared to the full master equation. We illustrate the breadth of the method on four superconducting-circuit tasks: preparation of a five-qubit GHZ state, a CNOT gate, qubit readout, and an error correction primitive, modeled with realistic multilevel transmons, decay, and strong drives, in each case reaching fidelities consistent with the intrinsic dissipation limits. LROC thereby extends pulse-level optimization to system sizes beyond the reach of existing gradient-based methods.
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.
We present a gradient-based optimal-control technique for open quantum systems that utilizes quantum trajectories to simulate the quantum dynamics during optimization. Using trajectories allows for optimizing open systems with less computational cost than the regular density matrix approaches in most realistic optimization problems. We introduce an improved-sampling algorithm which minimizes the number of trajectories needed per optimization iteration. Together with employing stochastic gradient descent techniques, this reduces the complexity of optimizing many realistic open quantum systems to the complexity encountered with closed systems. Our optimizer harnesses automatic differentiation to provide flexibility in optimization and to suit the different constraints and diverse parameter regimes of real-life experiments. We utilize the optimizer in a variety of applications to demonstrate how the use of quantum trajectories significantly reduces the computation complexity while achieving a multitude of simultaneous optimization targets. Demonstrated targets include high state-transfer fidelities despite dissipation, faster gate times, and maximization of qubit-readout fidelity while maintaining the quantum nondemolition nature of the measurement and allowing for subsequent fast resonator reset.
Abstract In this work, we adopt the gradient projection method (GPM) to problems of quantum control. For general N-level closed and open quantum systems, we derive the corresponding adjoint systems and gradients of the objective functionals and provide the projection versions of the Pontryagin maximum principle and the GPM, all directly in terms of quantum objects such as evolution operator, Hamiltonians, density matrices, etc. Various forms of the GPM, including one- and two-step, are provided and compared. We formulate the GPM both for closed and open quantum systems, latter for the general case with simultaneous coherent and incoherent controls. The GPM is designed to perform local gradient based optimization in the case when bounds are imposed on the controls. The main advantage of the method is that it allows to exactly satisfy the bounds, in difference to other approaches such as adding constraints as weight to an objective. We apply the GPM to several examples including generation of one- and two-qubit gates and two-qubit Bell and Werner states for models of superconducting qubits under the constraint when controls are zero at the initial and final times, steering an open quantum system to a target density matrix for simulating action of the Werner–Holevo channel, etc.
We develop a framework of "semi-automatic differentiation" that combines existing gradient-based methods of quantum optimal control with automatic differentiation. The approach allows to optimize practically any computable functional and is implemented in two open source Julia packages, GRAPE.jl and Krotov.jl , part of the QuantumControl.jl framework. Our method is based on formally rewriting the optimization functional in terms of propagated states, overlaps with target states, or quantum gates. An analytical application of the chain rule then allows to separate the time propagation and the evaluation of the functional when calculating the gradient. The former can be evaluated with great efficiency via a modified GRAPE scheme. The latter is evaluated with automatic differentiation, but with a profoundly reduced complexity compared to the time propagation. Thus, our approach eliminates the prohibitive memory and runtime overhead normally associated with automatic differentiation and facilitates further advancement in quantum control by enabling the direct optimization of non-analytic functionals for quantum information and quantum metrology, especially in open quantum systems. We illustrate and benchmark the use of semi-automatic differentiation for the optimization of perfectly entangling quantum gates on superconducting qubits coupled via a shared transmission line. This includes the first direct optimization of the non-analytic gate concurrence.
Optimal control of open quantum systems governed by the Lindblad master equation requires the repeated solution of forward and adjoint evolution equations. For large Hilbert spaces, these computations demand numerical methods that are structure-preserving, accurate, and computationally efficient, while remaining compatible with nonsmooth optimization techniques. In this work, a unified framework for optimal control of Lindblad equations is presented. The continuous optimality system is derived from the Pontryagin maximum principle and solved by a sequential quadratic Hamiltonian (SQH) method. Its numerical realization is based on second-order exponential midpoint propagators for the forward and adjoint Lindblad equations with time-dependent Hamiltonians. Full-rank schemes preserving the Hermitian and positive-semidefinite structure of the density matrix are developed together with low-rank formulations that substantially reduce storage requirements and computational cost. Rigorous error estimates are established for the full- and low-rank forward and adjoint propagators. Numerical experiments confirm the predicted convergence rates, demonstrate the effectiveness of the low-rank approximations, and illustrate the performance of the proposed FREM-SQH and LREM-SQH algorithms for optimal control problems with smooth and nonsmooth control costs.
Gradient-based quantum optimal control and the theory of operator evolution are rarely discussed together. In this Letter, we bridge this gap by showing that the gradient of the time-evolved propagator involves the Heisenberg evolution of local operators. We present a unifying framework by deriving from first principles the formal solution for the gradient under an arbitrary pulse parameterization, and show that the same construction extends to derivatives of any order. We obtain a series of nested, time-independent commutators weighted by time-dependent scalar coefficients. The commutators are thus computed once, and each optimization step only updates the scalars. The method is particularly suited for simulating optimal control tasks in quantum systems with local interactions, which is a common situation in large multi-qubit platforms. In this setting, light-cone arguments and sparse-Pauli truncation heuristics may be used to keep the number of relevant terms in the series small. We compare the computational cost required for the series with the Gradient Optimization of Analytic conTrols (GOAT) method, and, focusing on the problem of preparation of a GHZ state, demonstrate more than an order of magnitude speedup for a qubit ladder and a chain geometry.