Classical simulations of time-dependent quantum systems are widely used in quantum control research. In particular, these simulations are commonly used to host iterative optimal control algorithms. This is convenient for algorithms that are too onerous to run in the loop with current-day quantum hardware, as well as for researchers without consistent access to hardware. However, if the model used to represent the system is not selected carefully, an optimized control protocol may be rendered futile when applied to hardware. We present a series of models, ordered in a hierarchy of progressive approximation, which appear in quantum control literature. The validity of each model is characterized experimentally by designing and benchmarking control protocols for an IBMQ cloud quantum device. This result demonstrates error amplification induced by the application of a first-order perturbative approximation. Furthermore, the emergence of errors that cannot be corrected by simple amplitude scaling of control pulses is demonstrated in simulation, due to an underlying mistreatment of noncomputational dynamics. Finally, an evaluation of simulated control dynamics reveals that despite the substantial variance in numerical predictions across the proposed models, the complexity of discovering local optimal control protocols appears invariant in the simple control scheme setting.
Classical simulations of time-dependent quantum systems are widely used in quantum control research. In particular, these simulations are commonly used to host iterative optimal control algorithms. This is convenient for algorithms that are too onerous to run in the loop with current-day quantum hardware, as well as for researchers without consistent access to hardware. However, if the model used to represent the system is not selected carefully, an optimized control protocol may be rendered futile when applied to hardware. We present a series of models, ordered in a hierarchy of progressive approximation, which appear in quantum control literature. The validity of each model is characterized experimentally by designing and benchmarking control protocols for an IBMQ cloud quantum device. This result demonstrates error amplification induced by the application of a first-order perturbative approximation. Furthermore, the emergence of errors that cannot be corrected by simple amplitude scaling of control pulses is demonstrated in simulation, due to an underlying mistreatment of noncomputational dynamics. Finally, an evaluation of simulated control dynamics reveals that despite the substantial variance in numerical predictions across the proposed models, the complexity of discovering local optimal control protocols appears invariant in the simple control scheme setting.
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 develop a simulator for quantum computers composed of superconducting transmon qubits. The simulation model supports an arbitrary number of transmons and resonators. Quantum gates are implemented by time-dependent pulses. Nontrivial effects such as crosstalk, leakage to non-computational states, entanglement between transmons and resonators, and control errors due to the pulses are inherently included. The time evolution of the quantum computer is obtained by solving the time-dependent Schrödinger equation. The simulation algorithm shows excellent scalability on high-performance supercomputers. We present results for the simulation of up to 16 transmons and resonators. Additionally, the model can be used to simulate environments, and we demonstrate the transition from an isolated system to an open quantum system governed by a Lindblad master equation. We also describe a procedure to extract model parameters from electromagnetic simulations or experiments. We compare simulation results to experiments on several NISQ processors of the IBM Q Experience. We find nearly perfect agreement between simulation and experiment for quantum circuits designed to probe crosstalk in transmon systems. By studying common gate metrics such as the fidelity or the diamond distance, we find that they cannot reliably predict the performance of repeated gate applications or practical quantum algorithms. As an alternative, we find that the results from two-transmon gate set tomography have an exceptional predictive power. Finally, we test a protocol from the theory of quantum error correction and fault tolerance. We find that the protocol systematically improves the performance of transmon quantum computers in the presence of characteristic control and measurement errors.
Robust control of a quantum system is essential to utilize the current noisy quantum hardware to its full potential, such as quantum algorithms. To achieve such a goal, a systematic search for an optimal control for any given experiment is essential. The design of optimal control pulses requires accurate numerical models and, therefore, accurate characterization of the system parameters. We present an online Bayesian approach for quantum characterization of qutrit systems, which automatically and systematically identifies optimal experiments that provide maximum information on the system parameters, thereby greatly reducing the number of experiments that need to be performed on the quantum testbed. Unlike most characterization protocols that provide point-estimates of the parameters, the proposed approach is able to estimate their probability distribution. The applicability of the Bayesian experimental design technique was demonstrated on test problems, where each experiment was defined by a parameterized control pulse. In addition to this, we also present an approach for iterative pulse extension, which is robust under uncertainties in transition frequencies and coherence times, and shot noise, despite being initialized with wide uninformative priors. Furthermore, we provide a mathematical proof of the theoretical identifiability of the model parameters and present conditions on the quantum state under which the parameters are identifiable. The proof and conditions for identifiability are presented for both closed and open quantum systems using the Schrödinger equation and the Lindblad master equation, respectively.
Quantum optimal control experiments and simulations have successfully manipulated the dynamics of systems ranging from atoms to biomolecules. Surprisingly, these collective works indicate that the effort (i.e., the number of algorithmic iterations) required to find an optimal control field appears to be essentially invariant to the complexity of the system. The present work explores this matter in a series of systematic optimizations of the state-to-state transition probability on model quantum systems with the number of states N ranging from 5 through 100. The optimizations occur over a landscape defined by the transition probability as a function of the control field. Previous theoretical studies on the topology of quantum control landscapes established that they should be free of suboptimal traps under reasonable physical conditions. The simulations in this work include nearly 5000 individual optimization test cases, all of which confirm this prediction by fully achieving optimal population transfer of at least 99.9 on careful attention to numerical procedures to ensure that the controls are free of constraints. Collectively, the simulation results additionally show invariance of required search effort to system dimension N. This behavior is rationalized in terms of the structural features of the underlying control landscape. The very attractive observed scaling with system complexity may be understood by considering the distance traveled on the control landscape during a search and the magnitude of the control landscape slope. Exceptions to this favorable scaling behavior can arise when the initial control field fluence is too large or when the target final state recedes from the initial state as N increases.
For paving the way to novel applications in quantum simulation, computation, and technology, increasingly large quantum systems have to be steered with high precision. It is a typical task amenable to numerical optimal control to turn the time course of pulses, i.e., piecewise constant control amplitudes, iteratively into an optimized shape. Here, we present a comparative study of optimal-control algorithms for a wide range of finite-dimensional applications. We focus on the most commonly used algorithms: GRAPE methods which update all controls concurrently, and Krotov-type methods which do so sequentially. Guidelines for their use are given and open research questions are pointed out. Moreover, we introduce a unifying algorithmic framework, DYNAMO (dynamic optimization platform), designed to provide the quantum-technology community with a convenient matlab-based tool set for optimal control. In addition, it gives researchers in optimal-control techniques a framework for benchmarking and comparing newly proposed algorithms with the state of the art. It allows a mix-and-match approach with various types of gradients, update and step-size methods as well as subspace choices. Open-source code including examples is made available at http://qlib.info.
Accurately controlling the state of a transmon qubit via applied microwave pulses is essential in successfully running quantum computers that use these qubits. Currently, control pulses are designed using classical simulations of time-dependent quantum systems that neglect any back-action of the transmon on the transmission lines that deliver the microwaves. In this work, we develop a semiclassical model for a transmon capacitively coupled to a transmission line that treats the transmon quantum mechanically and the transmission line classically. This model can be used in the future to study the detailed interactions of microwave control pulses with transmon qubits to better optimize the control pulses. We show results from a classical model of this system to highlight that the back-action of the transmon may not be negligible for common experimental parameters, emphasizing the need for more rigorous models as experimental parameters are changed to accommodate scaling in device density.