Optimization of the fidelity of control operations is of critical importance in the pursuit of fault tolerant quantum computation. We apply optimal control techniques to demonstrate that a single drive via the cavity in circuit quantum electrodynamics can implement a high fidelity two-qubit all-microwave gate that directly entangles the qubits via the mutual qubit-cavity couplings. This is performed by driving at one of the qubits' frequencies which generates a conditional two-qubit gate, but will also generate other spurious interactions. These optimal control techniques are used to find pulse shapes that can perform this two-qubit gate with high fidelity, robust against errors in the system parameters. The simulations were all performed using experimentally relevant parameters and constraints.
Characterizing complex quantum systems is a vital task in quantum information science. Quantum tomography, the standard tool used for this purpose, uses a well-designed measurement record to reconstruct quantum states and processes. It is, however, notoriously inefficient. Recently, the classical signal reconstruction technique known as "compressed sensing" has been ported to quantum information science to overcome this challenge: accurate tomography can be achieved with substantially fewer measurement settings, thereby greatly enhancing the efficiency of quantum tomography. Here we show that compressed sensing tomography of quantum systems is essentially guaranteed by a special property of quantum mechanics itself---that the mathematical objects that describe the system in quantum mechanics are matrices with nonnegative eigenvalues. This result has an impact on the way quantum tomography is understood and implemented. In particular, it implies that the information obtained about a quantum system through compressed sensing methods exhibits a new sense of "informational completeness." This has important consequences on the efficiency of data taking for quantum tomography, and enables us to construct informationally complete measurements that are robust to noise and modeling errors. Moreover, our result shows that one can expand the numerical tool-box used in quantum tomography and employ highly efficient algorithms developed to handle large dimensional matrices on a large dimensional Hilbert space. While we mainly present our results in the context of quantum tomography, they apply to the general case of positive semidefinite matrix recovery.
We apply the method of compressed sensing (CS) quantum process tomography (QPT) to characterize quantum gates based on superconducting Xmon and phase qubits. Using experimental data for a two-qubit controlled-Z gate, we obtain an estimate for the process matrix $\chi$ with reasonably high fidelity compared to full QPT, but using a significantly reduced set of initial states and measurement configurations. We show that the CS method still works when the amount of used data is so small that the standard QPT would have an underdetermined system of equations. We also apply the CS method to the analysis of the three-qubit Toffoli gate with numerically added noise, and similarly show that the method works well for a substantially reduced set of data. For the CS calculations we use two different bases in which the process matrix $\chi$ is approximately sparse, and show that the resulting estimates of the process matrices match each ther with reasonably high fidelity. For both two-qubit and three-qubit gates, we characterize the quantum process by not only its process matrix and fidelity, but also by the corresponding standard deviation, defined via variation of the state fidelity for different initial states.
Resource tradeoffs can often be established by solving an appropriate robust optimization problem for a variety of scenarios involving constraints on optimization variables and uncertainties. Using an approach based on sequential convex programming, we demonstrate that quantum gate transformations can be made substantially robust against uncertainties while simultaneously using limited resources of control amplitude and bandwidth. Achieving such a high degree of robustness requires a quantitative model that specifies the range and character of the uncertainties. Using a model of a controlled one-qubit system for illustrative simulations, we identify robust control fields for a universal gate set and explore the tradeoff between the worst-case gate fidelity and the field fluence. Our results demonstrate that, even for this simple model, there exist a rich variety of control design possibilities. In addition, we study the effect of noise represented by a stochastic uncertainty model.
We develop an efficient and robust approach to Hamiltonian identification for multipartite quantum systems based on the method of compressed sensing. This work demonstrates that with only O(s log(d)) experimental configurations, consisting of random local preparations and measurements, one can estimate the Hamiltonian of a d-dimensional system, provided that the Hamiltonian is nearly s-sparse in a known basis. We numerically simulate the performance of this algorithm for three- and four-body interactions in spin-coupled quantum dots and atoms in optical lattices. Furthermore, we apply the algorithm to characterize Hamiltonian fine structure and unknown system-bath interactions.
The resources required to characterise the dynamics of engineered quantum systems-such as quantum computers and quantum sensors-grow exponentially with system size. Here we adapt techniques from compressive sensing to exponentially reduce the experimental configurations required for quantum process tomography. Our method is applicable to dynamical processes that are known to be nearly-sparse in a certain basis and it can be implemented using only single-body preparations and measurements. We perform efficient, high-fidelity estimation of process matrices on an experiment attempting to implement a photonic two-qubit logic-gate. The data base is obtained under various decoherence strengths. We find that our technique is both accurate and noise robust, thus removing a key roadblock to the development and scaling of quantum technologies.
We propose and evaluate experimentally an approach to quantum process tomography that completely removes the scaling problem plaguing the standard approach. The key to this simplification is the incorporation of prior knowledge of the class of physical interactions involved in generating the dynamics, which reduces the problem to one of parameter estimation. This allows part of the problem to be tackled using efficient convex methods, which, when coupled with a constraint on some parameters allows globally optimal estimates for the Kraus operators to be determined from experimental data. Parameterising the maps provides further advantages: it allows the incorporation of mixed states of the environment as well as some initial correlation between the system and environment, both of which are common physical situations following excitation of the system away from thermal equilibrium. Although the approach is not universal, in cases where it is valid it returns a complete set of positive maps for the dynamical evolution of a quantum system at all times.
The problem of quantifying the difference between evolutions of an open quantum system (in particular, between the actual evolution of an open system and the ideal target operation on the corresponding closed system) is important in quantum control, especially in control of quantum information processing. Motivated by this problem, we develop a measure for evaluating the distance between unitary evolution operators of a composite quantum system that consists of a sub-system of interest (e.g., a quantum information processor) and environment. The main characteristic of this measure is the invariance with respect to the effect of the evolution operator on the environment, which follows from an equivalence relation that exists between unitary operators acting on the composite system, when the effect on only the sub-system of interest is considered. The invariance to the environment's transformation makes it possible to quantitatively compare the evolution of an open quantum system and its closed counterpart. The distance measure also determines the fidelity bounds of a general quantum channel (a completely positive and trace-preserving map acting on the sub-system of interest) with respect to a unitary target transformation. This measure is also independent of the initial state of the system and straightforward to numerically calculate. As an example, the measure is used in numerical simulations to evaluate fidelities of optimally controlled quantum gate operations (for one- and two-qubit systems), in the presence of a decohering environment. This example illustrates the utility of this measure for optimal control of quantum operations in the realistic case of open-system dynamics.
We consider the problem of quantum multi-parameter estimation with experimental constraints and formulate the solution in terms of a convex optimization. Specifically, we outline an efficient method to identify the optimal strategy for estimating multiple unknown parameters of a quantum process and apply this method to a realistic example. The example is two electron spin qubits coupled through the dipole and exchange interactions with unknown coupling parameters -- explicitly, the position vector relating the two qubits and the magnitude of the exchange interaction are unknown. This coupling Hamiltonian generates a unitary evolution which, when combined with arbitrary single-qubit operations, produces a universal set of quantum gates. However, the unknown parameters must be known precisely to generate high-fidelity gates. We use the Cramér-Rao bound on the variance of a point estimator to construct the optimal series of experiments to estimate these free parameters, and present a complete analysis of the optimal experimental configuration. Our method of transforming the constrained optimal parameter estimation problem into a convex optimization is powerful and widely applicable to other systems.
For an initially well designed but imperfect quantum information system, the process matrix is almost sparse in an appropriate basis. Existing theory and associated computational methods (L1-norm minimization) for reconstructing sparse signals establish conditions under which the sparse signal can be perfectly reconstructed from a very limited number of measurements (resources). Although a direct extension to quantum process tomography of the L1-norm minimization theory has not yet emerged, the numerical examples presented here, which apply L1-norm minimization to quantum process tomography, show a significant reduction in resources to achieve a desired estimation accuracy over existing methods.
We develop a theory for finding quantum error correction (QEC) procedures which are optimized for given noise channels. Our theory accounts for uncertainties in the noise channel, against which our QEC procedures are robust. We demonstrate via numerical examples that our optimized QEC procedures always achieve a higher channel fidelity than the standard error correction method, which is agnostic about the specifics of the channel. This demonstrates the importance of channel characterization before QEC procedures are applied. Our main novel finding is that in the setting of a known noise channel the recovery ancillas are redundant for optimized quantum error correction. We show this using a general rank minimization heuristic and supporting numerical calculations. Therefore, one can further improve the fidelity by utilizing all the available ancillas in the encoding block.
Maximizing the precision in estimating parameters in a quantum system subject to instrumentation constraints is cast as a convex optimization problem. We account for prior knowledge about the parameter range by developing a worst-case and average case objective for optimizing the precision. Focusing on the single parameter case, we show that the optimization problems are \em linear programs. For the average case the solution to the linear program can be expressed analytically and involves a simple search: finding the largest element in a list. An example is presented which compares what is possible under constraints against the ideal with no constraints, the Quantum Fisher Information.
This work studies the feasibility of optimal control of high-fidelity quantum gates in a model of interacting two-level particles. One particle (the qubit) serves as the quantum information processor, whose evolution is controlled by a time-dependent external field. The other particles are not directly controlled and serve as an effective environment, coupling to which is the source of decoherence. The control objective is to generate target one-qubit gates in the presence of strong environmentally-induced decoherence and under physically motivated restrictions on the control field. It is found that interactions among the environmental particles have a negligible effect on the gate fidelity and require no additional adjustment of the control field. Another interesting result is that optimally controlled quantum gates are remarkably robust to random variations in qubit-environment and inter-environment coupling strengths. These findings demonstrate the utility of optimal control for management of quantum-information systems in a very precise and specific manner, especially when the dynamics complexity is exacerbated by inherently uncertain environmental coupling.
Mar 30 2007 quant-ph
We present a semidefinite program optimization approach to quantum error correction that yields codes and recovery procedures that are robust against significant variations in the noise channel. Our approach allows us to optimize the encoding, recovery, or both, and is amenable to approximations that significantly improve computational cost while retaining fidelity. We illustrate our theory numerically for optimized 5-qubit codes, using the standard [5,1,3] code as a benchmark. Our optimized encoding and recovery yields fidelities that are uniformly higher by 1-2 orders of magnitude against random unitary weight-2 errors compared to the [5,1,3] code with standard recovery. We observe similar improvement for a 4-qubit decoherence-free subspace code.
Feb 15 2007 quant-ph
Methods of optimal control are applied to a model system of interacting two-level particles (e.g., spin-half atomic nuclei or electrons or two-level atoms) to produce high-fidelity quantum gates while simultaneously negating the detrimental effect of decoherence. One set of particles functions as the quantum information processor, whose evolution is controlled by a time-dependent external field. The other particles are not directly controlled and serve as an effective environment, coupling to which is the source of decoherence. The control objective is to generate target one- and two-qubit unitary gates in the presence of strong environmentally-induced decoherence and under physically motivated restrictions on the control field. The quantum-gate fidelity, expressed in terms of a novel state-independent distance measure, is maximized with respect to the control field using combined genetic and gradient algorithms. The resulting high-fidelity gates demonstrate the feasibility of precisely guiding the quantum evolution via optimal control, even when the system complexity is exacerbated by environmental coupling. It is found that the gate duration has an important effect on the control mechanism and resulting fidelity. An analysis of the sensitivity of the gate performance to random variations in the system parameters reveals a significant degree of robustness attained by the optimal control solutions.
Nov 21 2006 quant-ph
This paper has been withdrawn as it has been superseded by postings arXiv:quant-ph/0702147 and arXiv:0712.2935
Jun 12 2006 quant-ph
We show that the problem of designing a quantum information error correcting procedure can be cast as a bi-convex optimization problem, iterating between encoding and recovery, each being a semidefinite program. For a given encoding operator the problem is convex in the recovery operator. For a given method of recovery, the problem is convex in the encoding scheme. This allows us to derive new codes that are locally optimal. We present examples of such codes that can handle errors which are too strong for codes derived by analogy to classical error correction techniques.
Jun 09 2006 quant-ph
A distance measure is presented between two unitary propagators of quantum systems of differing dimensions along with a corresponding method of computation. A typical application is to compare the propagator of the actual (real) process with the propagator of the desired (ideal) process; the former being of a higher dimension then the latter. The proposed measure has the advantage of dealing with possibly correlated inputs, but at the expense of working on the whole space and not just the information bearing part as is usually the case, i.e., no partial trace operation is explicitly involved. It is also shown that the distance measure and an average measure of channel fidelity both depend on the size of the same matrix: as the matrix size increases, distance decreases and fidelity increases.
Dec 09 2004 quant-ph
We present a formalism for encoding the logical basis of a qubit into subspaces of multiple physical levels. The need for this multilevel encoding arises naturally in situations where the speed of quantum operations exceeds the limits imposed by the addressability of individual energy levels of the qubit physical system. A basic feature of the multilevel encoding formalism is the logical equivalence of different physical states and correspondingly, of different physical transformations. This logical equivalence is a source of a significant flexibility in designing logical operations, while the multilevel structure inherently accommodates fast and intense broadband controls thereby facilitating faster quantum operations. Another important practical advantage of multilevel encoding is the ability to maintain full quantum-computational fidelity in the presence of mixing and decoherence within encoding subspaces. The formalism is developed in detail for single-qubit operations and generalized for multiple qubits. As an illustrative example, we perform a simulation of closed-loop optimal control of single-qubit operations for a model multilevel system, and subsequently apply these operations at finite temperatures to investigate the effect of decoherence on operational fidelity.
Nov 15 2004 quant-ph
A number of problems in quantum state and system identification are addressed. Specifically, it is shown that the maximum likelihood estimation (MLE) approach, already known to apply to quantum state tomography, is also applicable to quantum process tomography (estimating the Kraus operator sum representation (OSR)), Hamiltonian parameter estimation, and the related problems of state and process (OSR) distribution estimation. Except for Hamiltonian parameter estimation, the other MLE problems are formally of the same type of convex optimization problem and therefore can be solved very efficiently to within any desired accuracy. Associated with each of these estimation problems, and the focus of the paper, is an optimal experiment design (OED) problem invoked by the Cramer-Rao Inequality: find the number of experiments to be performed in a particular system configuration to maximize estimation accuracy; a configuration being any number of combinations of sample times, hardware settings, prepared initial states, etc. We show that in all of the estimation problems, including Hamiltonian parameter estimation, the optimal experiment design can be obtained by solving a convex optimization problem. Software to solve the MLE and OED convex optimization problems is available upon request from the first author.
Mar 23 2004 quant-ph
The problem addressed is to design a detector which is maximally sensitive to specific quantum states. Here we concentrate on quantum state detection using the worst-case a posteriori probability of detection as the design criterion. This objective is equivalent to asking the question: if the detector declares that a specific state is present, what is the probability of that state actually being present? We show that maximizing this worst-case probability (maximizing the smallest possible value of this probability) is a quasiconvex optimization over the matrices of the POVM (positive operator valued measure) which characterize the measurement apparatus. We also show that with a given POVM, the optimization is quasiconvex in the matrix which characterizes the Kraus operator sum representation (OSR) in a fixed basis. We use Lagrange Duality Theory to establish the optimality conditions for both deterministic and randomized detection. We also examine the special case of detecting a single pure state. Numerical aspects of using convex optimization for quantum state detection are also discussed.