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:2012.09282 · Physical Review Research

Fast and differentiable simulation of driven quantum systems

Ross Shillito, Jonathan A. Gross, Agustín Di Paolo, Élie Genois, Alexandre Blais

Published Sep 20, 2021• 11 citations
arXivHTMLPDF

Abstract

The controls enacting logical operations on quantum systems are described by time-dependent Hamiltonians that often include rapid oscillations. In order to accurately capture the resulting time dynamics in numerical simulations, a very small integration time step is required, which can severely impact the simulation run time. Here, we introduce a semianalytic method based on the Dyson expansion that allows us to time-evolve driven quantum systems much faster than standard numerical integrators. This solver, which we name Dysolve, efficiently captures the effect of the highly oscillatory terms in the system Hamiltonian, significantly reducing the simulation's run time as well as its sensitivity to the time-step size. Furthermore, this solver provides the exact derivative of the time-evolution operator with respect to the drive amplitudes. This key feature allows for optimal control in the limit of strong drives and goes beyond common pulse-optimization approaches that rely on rotating-wave approximations. As an illustration of our method, we show results of the optimization of a two-qubit gate using transmon qubits in the circuit QED architecture.

DOI
10.1103/physrevresearch.3.033266

Related papers

Accurate Methods for the Analysis of Strong-Drive Effects in Parametric Gates

Alexandru Petrescu, Camille Le Calonnec, Catherine Leroux +6 more

The ability to perform fast, high-fidelity entangling gates is a requirement for a viable quantum processor. In practice, achieving fast gates often comes with the penalty of strong-drive effects that are not captured by the rotating-wave approximation. These effects can be analyzed in simulations of the gate protocol, but those are computationally costly and often hide the physics at play. Here, we show how to efficiently extract gate parameters by directly solving a Floquet eigenproblem using exact numerics and a perturbative analytical approach. As an example application of this toolkit, we study the space of parametric gates generated between two fixed-frequency transmon qubits connected by a parametrically driven coupler. Our analytical treatment, based on time-dependent Schrieffer-Wolff perturbation theory, yields closed-form expressions for gate frequencies and spurious interactions, and is valid for strong drives. From these calculations, we identify optimal regimes of operation for different types of gates including $i$swap, controlled-Z, and cnot. These analytical results are supplemented by numerical Floquet computations from which we directly extract drive-dependent gate parameters. This approach has a considerable computational advantage over full simulations of time evolutions. More generally, our combined analytical and numerical strategy allows us to characterize two-qubit gates involving parametrically driven interactions, and can be applied to gate optimization and cross-talk mitigation such as the cancelation of unwanted $ZZ$ interactions in multiqubit architectures.

arXiv:2107.02343•Physical Review Applied•Apr 3, 2023•39 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

Quantum Algorithm for Time-Dependent Hamiltonian Simulation by Permutation Expansion

Yi‐Hsiang Chen, Amir Kalev, Itay Hen

We present a quantum algorithm for the dynamical simulation of time-dependent Hamiltonians.Our method involves expanding the interaction-picture Hamiltonian as a sum of generalized permutations, which leads to an integral-free Dyson series of the time-evolution operator.Under this representation, we perform a quantum simulation for the time-evolution operator by means of the linear combination of unitaries technique.We optimize the time steps of the evolution based on the Hamiltonian's dynamical characteristics, leading to a gate count that scales with an L 1 -norm-like scaling with respect only to the norm of the interaction Hamiltonian, rather than that of the total Hamiltonian.We demonstrate that the cost of the algorithm is independent of the Hamiltonian's frequencies, implying its advantage for systems with highly oscillating components, and for time-decaying systems the cost does not scale with the total evolution time asymptotically.In addition, our algorithm retains the near optimal log(1/ )/ log log(1/ ) scaling with simulation error .

arXiv:2103.15334•PRX Quantum•Sep 9, 2021•31 citations

Quandary

N. Anders Petersson, Stefanie Guenther

Quandary numerically simulates and optimizes the time-evolution of open quantum systems. The underlying dynamics are modelled by Lindblad's master equation, a linear ordinary differential equation (ODE) describing quantum systems interacting with the environment. Quandary solves this ODE numerically by applying a time-stepping integration scheme, and utilizes a gradient-based optimization approach to determine optimal control pulses that drive the quantum system to a desired target state. Two optimization objectives are considered: (a) Unitary gate optimization that finds controls to realize a unitary gate transformation, and (b) optimal reset that aims to drive the quantum system to the ground states. Gradient-based optimization schemes utilizing Petsc's Tao optimization package are applied to generate control pulses that minimize the respective measure. To evaluate the gradient of the objective function, the discrete adjoint method is used while leveraging techniques from Algorithmic Differentiation to produce exact and consistent gradients. To mitigate excessive execution run times, the software can be build together with the XBraid software library which provides a parallelization strategy to distribute the time-evolution of the underlying dynamics onto multiple processor.

•Jan 1, 2020•0 citations

Dynamics of a two-level system under strong driving: Quantum-gate optimization based on Floquet theory

Chunqing Deng, Feiruo Shen, Sahel Ashhab +1 more

We consider the dynamics of a two-level system (qubit) driven by strong and short resonant pulses in the framework of Floquet theory. First we derive analytical expressions for the quasienergies and Floquet states of the driven system. If the pulse amplitude varies very slowly, the system adiabatically follows the instantaneous Floquet states, which acquire dynamical phases that depend on the evolution of the quasienergies over time. The difference between the phases acquired by the two Floquet states corresponds to a qubit state rotation, generalizing the notion of Rabi oscillations to the case of large driving amplitudes. If the pulse amplitude changes very fast, the evolution is nonadiabatic, with transitions taking place between the Floquet states. We quantify and analyze the nonadiabatic transitions during the pulse by employing adiabatic perturbation theory and exact numerical simulations. We find that, for certain combinations of pulse rise and fall times and maximum driving amplitude, a destructive interference effect leads to a remarkably strong suppression of transitions between the Floquet states. This effect provides the basis of a quantum control protocol, which we name Floquet interference efficient suppression of transitions in the adiabatic basis (FIESTA), that can be used to design ultrafast high-fidelity single-qubit quantum gates.

arXiv:1605.08826•Physical Review A•Sep 27, 2016•45 citations

Time-dependent Hamiltonian Simulation of Highly Oscillatory Dynamics.

Dong An, Di Fang, Lin Lin

We propose a simple quantum algorithm for simulating highly oscillatory quantum dynamics, which does not require complicated quantum control logic for handling time-ordering operators. To our knowledge, this is the first quantum algorithm that is both insensitive to the rapid changes of the time-dependent Hamiltonian and exhibits commutator scaling. Our method can be used for efficient Hamiltonian simulation in the interaction picture. In particular, we demonstrate that for the simulation of the Schr\odinger equation, our method exhibits superconvergence and achieves a surprising second order convergence rate, of which the proof rests on a careful application of pseudo-differential calculus. Numerical results verify the effectiveness and the superconvergence property of our method.

•Nov 7, 2021•1 citation
The controls enacting logical operations on quantum systems are described by time-dependent Hamiltonians that often include rapid oscillations. In order to accurately capture the resulting time dynamics in numerical simulations, a very small integration time step is required, which can severely impact the simulation run time. Here, we introduce a semianalytic method based on the Dyson expansion that allows us to time-evolve driven quantum systems much faster than standard numerical integrators. This solver, which we name Dysolve, efficiently captures the effect of the highly oscillatory terms in the system Hamiltonian, significantly reducing the simulation's run time as well as its sensitivity to the time-step size. Furthermore, this solver provides the exact derivative of the time-evolution operator with respect to the drive amplitudes. This key feature allows for optimal control in the limit of strong drives and goes beyond common pulse-optimization approaches that rely on rotating-wave approximations. As an illustration of our method, we show results of the optimization of a two-qubit gate using transmon qubits in the circuit QED architecture.
The ability to perform fast, high-fidelity entangling gates is a requirement for a viable quantum processor. In practice, achieving fast gates often comes with the penalty of strong-drive effects that are not captured by the rotating-wave approximation. These effects can be analyzed in simulations of the gate protocol, but those are computationally costly and often hide the physics at play. Here, we show how to efficiently extract gate parameters by directly solving a Floquet eigenproblem using exact numerics and a perturbative analytical approach. As an example application of this toolkit, we study the space of parametric gates generated between two fixed-frequency transmon qubits connected by a parametrically driven coupler. Our analytical treatment, based on time-dependent Schrieffer-Wolff perturbation theory, yields closed-form expressions for gate frequencies and spurious interactions, and is valid for strong drives. From these calculations, we identify optimal regimes of operation for different types of gates including iiiswap, controlled-Z, and cnot. These analytical results are supplemented by numerical Floquet computations from which we directly extract drive-dependent gate parameters. This approach has a considerable computational advantage over full simulations of time evolutions. More generally, our combined analytical and numerical strategy allows us to characterize two-qubit gates involving parametrically driven interactions, and can be applied to gate optimization and cross-talk mitigation such as the cancelation of unwanted ZZZZZZ interactions in multiqubit architectures.
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.
We present a quantum algorithm for the dynamical simulation of time-dependent Hamiltonians.Our method involves expanding the interaction-picture Hamiltonian as a sum of generalized permutations, which leads to an integral-free Dyson series of the time-evolution operator.Under this representation, we perform a quantum simulation for the time-evolution operator by means of the linear combination of unitaries technique.We optimize the time steps of the evolution based on the Hamiltonian's dynamical characteristics, leading to a gate count that scales with an L 1 -norm-like scaling with respect only to the norm of the interaction Hamiltonian, rather than that of the total Hamiltonian.We demonstrate that the cost of the algorithm is independent of the Hamiltonian's frequencies, implying its advantage for systems with highly oscillating components, and for time-decaying systems the cost does not scale with the total evolution time asymptotically.In addition, our algorithm retains the near optimal log(1/ )/ log log(1/ ) scaling with simulation error .
Quandary numerically simulates and optimizes the time-evolution of open quantum systems. The underlying dynamics are modelled by Lindblad's master equation, a linear ordinary differential equation (ODE) describing quantum systems interacting with the environment. Quandary solves this ODE numerically by applying a time-stepping integration scheme, and utilizes a gradient-based optimization approach to determine optimal control pulses that drive the quantum system to a desired target state. Two optimization objectives are considered: (a) Unitary gate optimization that finds controls to realize a unitary gate transformation, and (b) optimal reset that aims to drive the quantum system to the ground states. Gradient-based optimization schemes utilizing Petsc's Tao optimization package are applied to generate control pulses that minimize the respective measure. To evaluate the gradient of the objective function, the discrete adjoint method is used while leveraging techniques from Algorithmic Differentiation to produce exact and consistent gradients. To mitigate excessive execution run times, the software can be build together with the XBraid software library which provides a parallelization strategy to distribute the time-evolution of the underlying dynamics onto multiple processor.
We consider the dynamics of a two-level system (qubit) driven by strong and short resonant pulses in the framework of Floquet theory. First we derive analytical expressions for the quasienergies and Floquet states of the driven system. If the pulse amplitude varies very slowly, the system adiabatically follows the instantaneous Floquet states, which acquire dynamical phases that depend on the evolution of the quasienergies over time. The difference between the phases acquired by the two Floquet states corresponds to a qubit state rotation, generalizing the notion of Rabi oscillations to the case of large driving amplitudes. If the pulse amplitude changes very fast, the evolution is nonadiabatic, with transitions taking place between the Floquet states. We quantify and analyze the nonadiabatic transitions during the pulse by employing adiabatic perturbation theory and exact numerical simulations. We find that, for certain combinations of pulse rise and fall times and maximum driving amplitude, a destructive interference effect leads to a remarkably strong suppression of transitions between the Floquet states. This effect provides the basis of a quantum control protocol, which we name Floquet interference efficient suppression of transitions in the adiabatic basis (FIESTA), that can be used to design ultrafast high-fidelity single-qubit quantum gates.
We propose a simple quantum algorithm for simulating highly oscillatory quantum dynamics, which does not require complicated quantum control logic for handling time-ordering operators. To our knowledge, this is the first quantum algorithm that is both insensitive to the rapid changes of the time-dependent Hamiltonian and exhibits commutator scaling. Our method can be used for efficient Hamiltonian simulation in the interaction picture. In particular, we demonstrate that for the simulation of the Schr\odinger equation, our method exhibits superconvergence and achieves a surprising second order convergence rate, of which the proof rests on a careful application of pseudo-differential calculus. Numerical results verify the effectiveness and the superconvergence property of our method.