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.26867

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

Ashutosh Mishra, Elena Lupo, Frank K. Wilhelm, Alessandro Ciani

quant-phSubmitted Jul 29, 2026• Announced Jul 30, 2026• 0 citations• Revised Sep 30, 2026 (v2)
v1 · Jul 29, 2026v2 · Sep 30, 2026
arXivHTMLPDF

Abstract

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.

Comments
8 pages, 3 figures; Supplemental Material: 17 pages, 2 figures

Related papers

Fast, high precision dynamics in quantum optimal control theory

Mogens Dalgaard, Felix Motzoi

Abstract Quantum optimal control theory is becoming increasingly crucial as quantum devices become more precise, but the need to quickly optimize these systems classically remains a significant bottleneck in their operation. Here we present a new theoretical quantum control framework for much faster optimization than the state of the art by replacing standard time propagation with a product of short-time propagators, each calculated using the Magnus expansion. The derived formulas for exact series terms and their gradients, based on earlier approximate integrals in a simulation setting, allow us to subsume the high cost of calculating commutators and integrals as an initial overhead. This provides an order of magnitude speedup for quantum control optimization.

arXiv:2110.06187•Journal of Physics B Atomic Molecular and Optical Physics•Apr 2, 2022•14 citations

Lie algebra-assisted quantum simulation and quantum optimal control via high-order Magnus expansions

R. F. dos Santos, S. J. J. M. F. Kokkelmans

The evolution of a quantum system under time-dependent driving exhibits phenomena that are absent in its stationary counterpart. However, the high dimensionality and non-commutative nature of quantum dynamics make this a challenging problem. The Magnus expansion provides an analytic framework to approximate the effective dynamics on short time-scales, but computing high-order terms with existing methods is computationally expensive. We introduce a scalable approach that reduces the computational effort to depend only on the degrees of freedom defining the time-dependent control function. We focus specifically on Hamiltonians consisting of a constant drift term and a controllable term. Our method provides a polynomial expression for the Magnus expansion which can be evaluated several orders of magnitude faster than previous techniques, enabling broad applications in the realms of quantum simulation and quantum optimal control. We showcase an application of the method by designing control pulses for the 5-qubit phase gate on a neutral-atom platform utilizing Rydberg atoms.

quant-pharXiv:2512.20357•Physical Review A•Dec 24, 2025•0 citations

Analytic Filter-Function Derivatives for Quantum Optimal Control

Isabel Nha Minh Le, Julian D. Teske, Tobias Hangleiter +2 more

Autocorrelated noise appears in many solid-state qubit systems and hence needs to be taken into account when developing gate operations for quantum information processing. However, explicitly simulating this kind of noise is often less efficient than approximate methods. Here, we focus on the filter function formalism, which allows the computation of gate fidelities in the presence of autocorrelated classical noise. Hence, this formalism can be combined with optimal control algorithms to design control pulses, which optimally implement quantum gates. To enable the use of gradient-based algorithms with fast convergence, we present analytically derived filter function gradients with respect to control pulse amplitudes, and analyze the computational complexity of our results. When comparing pulse optimization using our derivatives to a gradient-free approach, we find that the gradient-based method is roughly 2 orders of magnitude faster for our test cases. We also provide a modular computational implementation compatible with quantum optimal control packages.

arXiv:2103.09126•Physical Review Applied•Feb 2, 2022•15 citations

Quantum optimal control via gradient ascent in function space and the time-bandwidth quantum speed limit

Dennis G. Lucarelli

A gradient ascent method for optimal quantum control synthesis is presented that employs a gradient derived with respect to the coefficients of a functional basis expansion of the control. Restricting the space of allowable controls to weighted sums of the Slepian sequences efficiently parametrizes the control in terms of bandwidth, resolution, and pulse duration. A bound showing minimum time evolutions scaling with the inverse of the control bandwidth [S. Lloyd and S. Montangero, Phys. Rev. Lett. 113, 010502 (2014)] is recovered and the method is shown to achieve the scaling predicted by the bound on entangling two-qubit quantum gates.

arXiv:1611.00188•Physical Review A•Jun 29, 2018•55 citations

Approximate dynamics leading to more optimal control: Efficient exact derivatives

Jesper Hasseriis Mohr Jensen, Frederik Møller, Jens Jakob Sørensen +1 more

Accurate derivatives are important for locally traversing and converging efficiently in quantum optimization landscapes. By deriving analytically exact control derivatives (gradient and Hessian) for unitary control tasks, we show here that the computational feasibility of meeting this accuracy requirement depends on the choice of propagation scheme and problem representation. Even when exact propagation is sufficiently cheap it is, perhaps surprisingly, much more efficient to optimize the (appropriately) approximate propagators: Approximations in the dynamics are traded off for significant complexity reductions in the exact derivative calculations. Importantly, past the initial analytical considerations, only standard numerical techniques are explicitly required with straightforward application to realistic systems. These results are numerically verified for two concrete problems of increasing Hilbert space dimensionality. The best schemes obtain unit fidelity to machine precision, whereas the results for other schemes are separated consistently by orders of magnitude in computation time and in worst-case ten orders of magnitude in achievable fidelity. Since these gaps continually increase with system size and complexity, this methodology allows numerically efficient optimization of very high-dimensional dynamics, e.g., in many-body contexts, operating in the high-fidelity regime.

arXiv:2005.09943•Physical Review A•Jun 21, 2021•19 citations

Pulse engineering via projection of response functions at infinite nonlinear order

Lia Kley, Ludwig Mathey

Optimal control problems arise in a wide range of scientific disciplines, but the corresponding optimization algorithms often display a strong dependence on hyperparameters that significantly influence performance and convergence. For the optimal implementation of quantum algorithms, these challenges are further amplified by high-dimensional control landscapes and the need for high-fidelity operations. Here, we propose an algorithm for optimal control problems in quantum computing to efficiently generate high-fidelity control protocols for multi-qubit systems in a hyperparameter and gradient free manner. The method, referred to as Pulse Engineering via Projection of response functions at infinite nonlinear order (PEPRino), leverages the framework of response theory to navigate the control landscape to find high-fidelity implementations. This is achieved by determining the control landscape via response functions to infinite order, efficiently evaluated by resummation in terms of the first and second order response function. To demonstrate the approach, we apply it to quantum systems consisting of two and three qubits for the optimal implementation of the Quantum Fourier Transform (QFT). We benchmark the proposed algorithm against the Chopped Random Basis (CRAB) algorithm utilizing the Nelder-Mead method, focusing on the 2-qubit scenario. The results indicate faster convergence regarding iteration steps and computational time, highlighting the advantages of our approach.

quant-pharXiv:2607.24725•Jul 28, 2026•0 citations
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.
Abstract Quantum optimal control theory is becoming increasingly crucial as quantum devices become more precise, but the need to quickly optimize these systems classically remains a significant bottleneck in their operation. Here we present a new theoretical quantum control framework for much faster optimization than the state of the art by replacing standard time propagation with a product of short-time propagators, each calculated using the Magnus expansion. The derived formulas for exact series terms and their gradients, based on earlier approximate integrals in a simulation setting, allow us to subsume the high cost of calculating commutators and integrals as an initial overhead. This provides an order of magnitude speedup for quantum control optimization.
The evolution of a quantum system under time-dependent driving exhibits phenomena that are absent in its stationary counterpart. However, the high dimensionality and non-commutative nature of quantum dynamics make this a challenging problem. The Magnus expansion provides an analytic framework to approximate the effective dynamics on short time-scales, but computing high-order terms with existing methods is computationally expensive. We introduce a scalable approach that reduces the computational effort to depend only on the degrees of freedom defining the time-dependent control function. We focus specifically on Hamiltonians consisting of a constant drift term and a controllable term. Our method provides a polynomial expression for the Magnus expansion which can be evaluated several orders of magnitude faster than previous techniques, enabling broad applications in the realms of quantum simulation and quantum optimal control. We showcase an application of the method by designing control pulses for the 5-qubit phase gate on a neutral-atom platform utilizing Rydberg atoms.
Autocorrelated noise appears in many solid-state qubit systems and hence needs to be taken into account when developing gate operations for quantum information processing. However, explicitly simulating this kind of noise is often less efficient than approximate methods. Here, we focus on the filter function formalism, which allows the computation of gate fidelities in the presence of autocorrelated classical noise. Hence, this formalism can be combined with optimal control algorithms to design control pulses, which optimally implement quantum gates. To enable the use of gradient-based algorithms with fast convergence, we present analytically derived filter function gradients with respect to control pulse amplitudes, and analyze the computational complexity of our results. When comparing pulse optimization using our derivatives to a gradient-free approach, we find that the gradient-based method is roughly 2 orders of magnitude faster for our test cases. We also provide a modular computational implementation compatible with quantum optimal control packages.
A gradient ascent method for optimal quantum control synthesis is presented that employs a gradient derived with respect to the coefficients of a functional basis expansion of the control. Restricting the space of allowable controls to weighted sums of the Slepian sequences efficiently parametrizes the control in terms of bandwidth, resolution, and pulse duration. A bound showing minimum time evolutions scaling with the inverse of the control bandwidth [S. Lloyd and S. Montangero, Phys. Rev. Lett. 113, 010502 (2014)] is recovered and the method is shown to achieve the scaling predicted by the bound on entangling two-qubit quantum gates.
Accurate derivatives are important for locally traversing and converging efficiently in quantum optimization landscapes. By deriving analytically exact control derivatives (gradient and Hessian) for unitary control tasks, we show here that the computational feasibility of meeting this accuracy requirement depends on the choice of propagation scheme and problem representation. Even when exact propagation is sufficiently cheap it is, perhaps surprisingly, much more efficient to optimize the (appropriately) approximate propagators: Approximations in the dynamics are traded off for significant complexity reductions in the exact derivative calculations. Importantly, past the initial analytical considerations, only standard numerical techniques are explicitly required with straightforward application to realistic systems. These results are numerically verified for two concrete problems of increasing Hilbert space dimensionality. The best schemes obtain unit fidelity to machine precision, whereas the results for other schemes are separated consistently by orders of magnitude in computation time and in worst-case ten orders of magnitude in achievable fidelity. Since these gaps continually increase with system size and complexity, this methodology allows numerically efficient optimization of very high-dimensional dynamics, e.g., in many-body contexts, operating in the high-fidelity regime.
Optimal control problems arise in a wide range of scientific disciplines, but the corresponding optimization algorithms often display a strong dependence on hyperparameters that significantly influence performance and convergence. For the optimal implementation of quantum algorithms, these challenges are further amplified by high-dimensional control landscapes and the need for high-fidelity operations. Here, we propose an algorithm for optimal control problems in quantum computing to efficiently generate high-fidelity control protocols for multi-qubit systems in a hyperparameter and gradient free manner. The method, referred to as Pulse Engineering via Projection of response functions at infinite nonlinear order (PEPRino), leverages the framework of response theory to navigate the control landscape to find high-fidelity implementations. This is achieved by determining the control landscape via response functions to infinite order, efficiently evaluated by resummation in terms of the first and second order response function. To demonstrate the approach, we apply it to quantum systems consisting of two and three qubits for the optimal implementation of the Quantum Fourier Transform (QFT). We benchmark the proposed algorithm against the Chopped Random Basis (CRAB) algorithm utilizing the Nelder-Mead method, focusing on the 2-qubit scenario. The results indicate faster convergence regarding iteration steps and computational time, highlighting the advantages of our approach.