Quandary

N. Anders Petersson, Stefanie Guenther

Published Jan 1, 2020• 0 citations

Abstract

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.

DOI: 10.11578/dc.20211210.3

Related papers

Quandary: An open-source C++ package for high-performance optimal control of open quantum systems

Stefanie Günther, N. Anders Petersson, Jonathan L. DuBois

Quantum optimal control can be used to shape the control pulses for realizing unitary and non-unitary transformations of quantum states. These control pulses provide the fundamental interface between the quantum compiler and the quantum hardware. Most current software for quantum optimal control (e.g. Qutip [1] or Krotov [2]) is restricted to run on shared memory platforms, limiting their applicability to smaller quantum systems, in particular if interactions with the environment are taken into account. This paper gives an overview of the open-source code Quandary, which is designed to solve quantum control problems in larger open quantum systems modelled by Lindblad’s master equation. Implemented in C++, Quandary uses the message passing paradigm for distributed memory computers that enables scalability to large numbers of compute cores. Accompanied by numerical examples, this paper presents an overview on existing theoretical developments for open optimal quantum control realizing state-to-state transfer, unitary gate optimization as well as state-preparation, and presents the numerical tools and implementation aspect as realized in Quandary, for deployment on modern high-performance computing platforms.

arXiv:2110.10310IEEE/ACM International Workshop on Quantum Computing Software (QCS)Nov 1, 202121 citations

Quandary: An open-source package for high-performance optimal control of open quantum systems

N. Anders Petersson

Quantum optimal control can be used to shape the control pulses for realizing unitary and non-unitary transformations of quantum states.These control pulses provide the fundamental interface between the quantum compiler and the quantum hardware.Most current software for quantum optimal control (e.g.Qutip [1] or Krotov [2]) is restricted to run on shared memory platforms, limiting their applicability to smaller quantum systems, in particular if interactions with the environment are taken into account.This paper gives an overview of the open source code Quandary, which models open quantum systems using Lindblad's master equation.Quandary is designed to solve optimal control problems in larger open quantum systems using a message passing paradigm, implemented in C++ for distributed memory computers that enables scalability to large numbers of compute cores.

Nov 19, 20210 citations

FOKKER–PLANCK-BASED CONTROL OF A TWO-LEVEL OPEN QUANTUM SYSTEM

Mario Annunziato, Alfio Borzı̀

The control of a two-level open quantum system subject to dissipation due to environment interaction is considered. The evolution of this system is governed by a Lindblad master equation which is augmented by a stochastic term to model the effect of time-continuous measurements. In order to control this stochastic master equation model, a Fokker–Planck control framework is investigated. Within this strategy, the control objectives are defined based on the probability density functions of the two-level stochastic process and the controls are computed as minimizers of these objectives subject to the constraints represented by the Fokker–Planck equation. This minimization problem is characterized by an optimality system including the Fokker–Planck equation and its adjoint. This optimality system is approximated by a second-order accurate, stable, conservative and positive-preserving discretization scheme. The implementation of the resulting open-loop controls is realized with a receding-horizon algorithm over a sequence of time windows. Results of numerical experiments demonstrate the effectiveness of the proposed approach.

Mathematical Models and Methods in Applied SciencesFeb 18, 20137 citations

Gate generation for open quantum systems via a monotonic algorithm with time optimization

Paulo Sérgio Pereira da Silva, Pierre Rouchon

We present a monotonic numerical algorithm including time optimization for generating quantum gates for open systems. Such systems are assumed to be governed by Lindblad master equations for the density operators on a large Hilbert-space whereas the quantum gates are relative to a sub-space of small dimension. Starting from an initial seed of the control input, this algorithm consists in the repetition of the following two steps producing a new control input: (A) backwards integration of adjoint Lindblad-Master equations (in the Heisenberg-picture) from a set of final conditions encoding the quantum gate to generate; (B) forward integration of Lindblad-Master equations in closed-loop where a Lyapunov based control produced the new control input. The numerical stability is ensured by the stability of both the open-loop adjoint backward system and the forward closed-loop system. A clock-control input can be added to the usual control input. The obtained monotonic algorithm allows then to optimise not only the shape of the control input, but also the gate time. Preliminary numerical implementations indicate that this algorithm is well suited for cat-qubit gates, where Hilbert-space dimensions (2 for the Z-gate and 4 for the CNOT-gate) are much smaller than the dimension of the physical Hilbert-space involving mainly Fock-states (typically 20 or larger for a single cat-qubit). This monotonic algorithm, based on Lyapunov control techniques, is shown to have a straightforward interpretation in terms of optimal control: its stationary conditions coincides with the first-order optimality conditions for a cost depending linearly on the final values of the quantum states.

arXiv:2403.20028Mathematical Control and Related FieldsJan 1, 20241 citation

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.05541Physical Review AMay 20, 201986 citations

最速量子操作

陽介奥平, Yosuke Okudaira

We present a general framework for finding the time-optimal quantum evolution inspired by the classical problem called the brachistochrone problem.At first, the time-optimal control problem for achieving a pure state evolution from a given initial state to a desired final state is formulated.We minimize the duration time for this transition by finding the optimal control for the Hamiltonian in an available set subject to a certain constraint dictated by physical conditions.Then, it will be extended to the unitary evolution which has greater relevance for quantum computation.This provides a systematic way to make an efficient unitary gate on a certain physical system.To demonstrate this framework, we choose the anisotropic Heisenberg model and find a general solution.Since it is not always possible to find a general solution for the fundamental equation for the time-optimal control, we introduce the simulated annealing method as a simple numerical method to show that our framework can be used for the system whose general solution cannot be found.The numerical analysis based on the Krotov iteration method is introduced and applied to the time-optimal quantum control problem to solve the problem efficiently.The optimal realization of the quantum Fourier transform for two or more qubit systems are derived.We consider two generalizations of the quantum brachistochrone for mixed state evolutions.One is the unitary evolution with ancilla system ended up with the final measurement on the ancilla.The other is the repetition of measurement or decoherence model which can be represented by the evolution with the master equation in the Lindblad form.Also, the relation between the fidelity to the target and the duration time is considerd.A.4 Quantum speed limit . . . . . . . . .

Mar 1, 20080 citations