Simple pulses for elimination of leakage in weakly nonlinear qubits

Simple pulses for elimination of leakage in weakly nonlinear qubits


In realizations of quantum computing, a two-level system (qubit) is often singled out from the many levels of an anharmonic oscillator. In these cases, simple qubit control fails on short time scales because of coupling to leakage levels. We provide an easy to implement analytic formula that inhibits this leakage from any single-control analog or pixelated pulse. It is based on adding a second control that is proportional to the time-derivative of the first. For realistic parameters of superconducting qubits, this strategy reduces the error by an order of magnitude relative to the state of the art, all based on smooth and feasible pulse shapes. These results show that even weak anharmonicity is sufficient and in general not a limiting factor for implementing quantum gates.

03.67.Lx, 02.30.Y, 85.25.-j

Quantum information processing devices are paving new inroads into the understanding of coherent processes and their applicability to high-complexity computation. The fundamental building blocks for these devices are quantum bits (qubits), which are quantum two-level systems. Naturally occurring qubits such as spin particles and photons are difficult to both isolate and control. Fortunately, effective qubits can also be manufactured from an anharmonic multi-level system. These systems can be understood as a generalized atom, and all share the common problem of leakage out of the qubit subspace if the control bandwidth is comparable to the anharmonicity. Examples of such systems are superconducting qubits [1-15], optical lattices Maneshi:2008a (), and quantum dots Schliemann:2008a ().

In superconducting qubits Clarke:2008a (), researchers have demonstrated single qubit gates Lucero:2008a (); Chow:2008a (), two qubit gates Yamamoto:2003b (); Plantenberg:2007a (), and the use of a quantum bus Majer:2007a (); Sillanpaa:2007a (), suggesting the possibility of a scalable technology. Decoherence has dropped dramatically Martinis:2005a (); Houck:2008a (); Bertet:2005a () from the original designs Nakamura:1999a (), due to both device and operation improvement Martinis:2005a (); Koch:2007a (); Vion:2002a (). To date, the single-qubit accuracy of such systems has been in the range of 98% Lucero:2008a ()- 99.3% Chow:2008a (), which is a major achievement, yet not sufficient for quantum computing.

For large scale computation, gate errors are required below a fault-tolerance threshold, conservatively estimated at Knill:2005a (). In superconducting qubits such as the Phase Martinis:2002a (); Martinis:2005a (); Lucero:2008a () and the Transmon Koch:2007a (); Schreier:2008a (); Houck:2008a (); Chow:2008a () qubits, long pulses are limited by decoherence and short pulses by excitations out of the qubit space Chow:2008a (); Steffen:2003a (). Thus, one prevalent approach to achieve high fidelity gates is to increase the coherence time; alternatively, the approach we will take is to reduce the gate time by using optimized control pulses that avoid leakage.

Thus far, trying to minimize gate times through better control has been met with limited success. Little improvement has been achieved beyond what is possible with a Gaussian waveform Bauer:1984a (); Steffen:2003a (). Optimal control theory with only amplitude modulation has required pulses which have sharp features Rebentrost:2008a (); Safaei:2008a (); Jirari:2005a (). For example, the pulses in Steffen:2003a (); Rebentrost:2008a () require 3 pulse steps and 2 waiting periods, which requires currently impractical pulse shaping hardware. Furthermore, these pulses are limited to a minimal time, given by the inverse of the anharmonicity.

We propose to solve this problem using two quadrature controls for effective phase modulation, which allows us to achieve instantaneous cancellation of leakage using slowly varying controls. We present an analytical approach to accomplish this called Derivative Removal by Adiabatic Gate (DRAG). Then, we map the full potential of the technique using numerical refinement and consider the case when decoherence is included.

I The system

We consider the lowest three levels of a driven slightly anharmonic energy spectrum with only nearest level coupling; the first two levels comprise the qubit and the third level accounts for leakage, provided its population is minimized during the entire pulse. That is, in the lab frame, the Hamiltonian can be written as


where is the projector for the level, is the lowering operator, and the corresponding transition energies are denoted , with the ground energy set to zero. The intrinsic anharmonicity of the system is . The parameter weighs the relative strength of the to transition, and without loss of generality we will take and . Driving and control of the system is represented by


which is a single frequency carrier with two independent quadrature controls, and ; is the time taken for one gate operation.

Making the rotating wave approximation (RWA) and moving to the interaction frame with respect to the drive frequency , the Hamiltonian can be written as


where and . and are the detunings of the transitions with respect to the drive frequency . When the drive frequency is resonant with the qubit frequency (), is the anharmonicity.

In the ideal case, , the Hamiltonian and its evolution reduce to that of a qubit. Thus qubit operations,


can be trivially driven by the controls, Eq. (2), since the effect is , which is a unitary that acts only in the qubit subspace, and the global phases are irrelevant Rebentrost:2008a (). For example, the NOT gate (which without loss of generality we choose to focus on) is implemented by simply setting , . Then, the time-ordering operator can be dropped and any controls that satisfy lead to the desired gate .

However, in the more general case when is non-zero, the pulse can result in large unwanted leakage into the third level. Fourier analysis can show how this leakage is dictated by the anharmonicity Warren:1984a (); Steffen:2003a (). The standard approach to diminish this error is to use Gaussian () Bauer:1984a () or tangential () envelope shaping


where is the standard deviation for and the rise-time parameter for . is chosen such that the correct amount of rotation is implemented (e.g. for a NOT) while enforces the pulse to start and end at zero. The small frequency response bandwidth of the Gaussian envelope ensures little excitation at the leakage transition frequency. Additionally, Gaussians and tangential pulses have relatively simple pulse-shaping hardware requirements. This being said, these pulses still suffer large errors at short gate times, as shown by the solid blue line in Fig. 1, where the error is at a gate time of  ns (vertical dashed line), and scales as . To quantify this error, we use the gate fidelity averaged over all input states existing in the qubit Hilbert space, which using an argument similar to Bowdrey:2002a () gives


where is the actual process in the three dimensionional Hilbert space, are the six axial states on the Bloch sphere, and is defined in Eq. (4).

Figure 1: Error vs. pulse duration for Gaussian (=) (blue solid), Gaussian (=) with DRAG (red dot-dashed), tangential (=-1.4) with DRAG (green dashed), and GRAPE pulse (purple crosses with dotted guiding line). The standard Gaussian and DRAG pulses are continuous, while the GRAPE pulse has a pixel size of 1ns. Parameters are , MHz) (typical for the Phase Martinis:2002a (); Martinis:2005a (); Lucero:2008a () and Transmon Koch:2007a (); Schreier:2008a (); Houck:2008a (); Chow:2008a () qubits). The top axis is set to dimensionless units, .

Ii Derivative Removal by Adiabatic Gate (DRAG)

To go beyond these simple pulses we introduce an adiabatic transformation that allows us to work entirely in the qubit subspace. This transformation is


where is the adiabaticity parameter. For the gates to be identical in both frames we require , such that by the end of the pulse it has no net effect. This is achieved by setting . In this frame, to first order in , the effective Hamiltonian reads


Here, we see we can remove the imaginary inertial term and the ac-Stark shift (phase) error with the second quadrature and drive detuning set to


Thus, we eliminate the leakage (second line of Eq. (8)) to order , which in this frame comes about from the induced transition 1. This error can also be removed by adding the term to the transformation . This results in a higher order error being added to the transition, , which can easily be compensated by rescaling . Increasing the precision of the transformation up to 5th order results in


These pulses remove most of the leakage error, as demonstrated in Fig. 1, often with only the first order of corrections. The green dashed (optimized tangential) and red dot-dashed (optimized Gaussian) lines in Fig. 1 show that DRAG reduces the gate error by orders of magnitude (e.g. 4 orders at 6ns). In fact, we verified with randomized testing that DRAG is largely impervious to envelope shape, provided it has the right area and starts and ends at 0. Fig. 2A shows an optimized Gaussian.

If cannot be controlled in real time, phase ramping can be applied Patt:1992a (). It amounts to transforming to the frame given by transformation where . Then, alternatively, we can pick , , and .

Figure 2: A. 6ns analog Gaussian () with DRAG; B. 4ns digital GRAPE with no detuning and pixels at 1ns. In both cases, the blue solid line is the control and the green dashed is the control. The red dot-dashed line in A represents explicit detuning, . Other parameters are the same as in Fig. 1.

Iii Numerical optimization (GRAPE)

We apply optimal control theory to our system by employing the gradient ascent numerical optimization algorithm, GRAPE Khaneja:2005a (). We try to minimize the gate error, , with as defined in Eq. (6). This is equivalent (for dissipation-free evolution) to maximizing


where is the unitary evolution found by GRAPE Rebentrost:2008a (). GRAPE works by time-slicing the evolution into discrete time pixels of constant amplitudes and , as in Fig. 2B. The gradient of with respect to each control at each time step provides a direction in search space. That is, each iteration of the algorithm increments the controls at each time step of the pulse proportionally to the gradient with respect to the control at that time (including the sign). This unique search direction allows for fast optimization, particularly with few pixels, and the algorithm quickly converges to an optimal solution () for all gate times (see Fig. 3).

Restricting GRAPE to 1ns pixels (which is consistent with current experimental limitations 2) and setting the initial condition for the quadrature to be the Gaussian pulse, Eq. (5), the algorithm quickly converges to the optimal solution. After optimization, the shape of the Gaussian in the quadrature is largely unchanged. However the quadrature changes from 0 to a shape similar to the derivative of the Gaussian, see Fig. 2B. GRAPE essentially takes the idea of DRAG to infinite order and discretizes the pulse. The shapes of the pulses are smooth with no sudden rises or falls in the control amplitudes. The purple dotted line in Fig. 1 plots the gate error as a function of gate time for GRAPE. The extra optimization increases the performance over DRAG.

The two quadrature control GRAPE pulse outperforms its one quadrature control version Rebentrost:2008a (). Fig. 3 shows the minimal pulse duration to obtain gate errors smaller than vs the pixel width for both one (red solid) and two (blue dashed) controls. Here, we see that at small pixel widths (continuous limit) the one control saturates at a time of Rebentrost:2008a () (black dashed line), whereas the two controls allows for arbitrarily small pulse times. This is because, in the one control case, phase cancelation of the leakage is done through phase accumulation by the natural precession of the third level, whereas, for two controls, it can be done instantaneously. The insert of Fig. 3 shows such a continuous two control pulse; this pulse clearly displays a combination of DRAG and composite features 3. At the larger pixels (more experimentally realizable), the two controls offer a substantial improvement over one control. For example, at 1ns pixels, there is a factor of 3.5 improvement.

Figure 3: Pulse duration vs. pixel width for 1 control (blue solid) and 2 controls (red dashed). Pulse duration is calculated as the minimal time to bring the error down to . The insert shows the two control GRAPE result for (blue solid is , and green dashed is ). Other parameters are the same as in Fig. 1.

Iv Decoherence

In realistic physical systems, decoherence cannot be neglected. Here, we include this effect by simulating the evolution by a master equation for each optimal pulse. The master equation we use is


where is the density matrix, and represent relaxation and pure dephasing times respectively, and is the damping superoperator, defined as


For the simulation, we used parameters consistent with the Transmon qubit and as such take pure dephasing to be zero. The results are shown in Fig. 4. Fig. 4A plots the gate error () as a function of gate time for s. Here, we see there is an optimal value for the gate time (when the error due to decoherence is about the same as the error due to leakage) and for both the DRAG (green dashed) and the GRAPE (purple dotted) pulse, this is much less then the standard Gaussian pulse (blue solid) with an optimal error of and respectively. Fig. 4B plots this minimum as a function of where we observe that at current experimental values s Houck:2008a () DRAG and GRAPE still out preform the simple Gaussian by about an order of magnitude.

Figure 4: A. error vs. pulse duration for s and B. minimum error vs. for Gaussian (=) (blue solid), Gaussian (=) with DRAG (red dot-dashed), and GRAPE (purple dotted) with 1ns pixels. Other parameters are the same as in Fig. 1 and pulses at the optimal times in A are plotted in Fig. 2.

V Conclusions

The analytical and numerical results presented here demonstrate that it is possible to eliminate leakage in physical qubits by using a second quadrature control approximately equal to the derivative of the first. The predicted single-qubit error rates are extremely low, even when we include decoherence. While in this paper this procedure is only applied to leakage, it will also find application to other systems, for example turning off unwanted coupling in multi-qubit and qubit-oscillator systems for single and two qubit operations.

We acknowledge W. A. Coish, Jens Koch, B. Osberg, C. A. Ryan, and M. Laforest for valuable discussions. JMG was supported by CIFAR, MITACS, and ORDCF. FM and FKW were supported by NSERC through the discovery grants and QuantumWorks.


  1. Note this approach is equivalent to a perturbation using the time-averaged Hamiltonian theory, where . An adiabatic expansion of in terms of is used instead of a Magnus expansion of the unitary operator Warren:1984a ().
  2. R. J. Schoelkopf, private communication (2008)
  3. At this really short gate time the RWA is no longer valid for the Phase Martinis:2002a (); Martinis:2005a (); Lucero:2008a () and the Transmon Koch:2007a (); Schreier:2008a (); Houck:2008a (); Chow:2008a () qubits.


  1. J. Clarke and F. K. Wilhelm, Nature 453, 1031 (2008).
  2. Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature 398, 786 (1999).
  3. D. Vion et al., Science 296, 886 (2002).
  4. T. Yamamoto et al., Nature 425, 941 (2003).
  5. J. H. Plantenberg et al., Nature 447, 836 (2007).
  6. M. A. Sillanpaa, J. I. Park, and R. W. Simmonds, Nature 449, 438 (2007).
  7. J. Majer et al., Nature 449, 443 (2007).
  8. J. M. Martinis et al., Phys. Rev. Lett. 89, 117901 (2002).
  9. J. M. Martinis et al., Phys. Rev. Lett. 95, 210503 (2005).
  10. E. Lucero et al., Phys. Rev. Lett. 100, 247001 (2008).
  11. J. M.Chow et al., Phys. Rev. Lett. 102, 090502 (2009)
  12. J. Koch et al., Phys. Rev. A 76, 042319 (2007).
  13. J. A. Schreier et al., Phys. Rev. B 77, 180502(R) (2008).
  14. A. A. Houck et al., Phys. Rev. Lett. 101, 080502 (2008).
  15. P. Bertet et al., Phys. Rev. Lett. 95, 257002 (2005).
  16. S. Maneshi et al., Phys. Rev. A 77, 022303 (2008).
  17. J. Schliemann, D. Loss, and A. H. MacDonald, Phys. Rev. B 63 085311 (2001).
  18. E. Knill, Nature 434, 39 (2005).
  19. M. Steffen, J. M. Martinis, and I. L. Chuang, Phys. Rev. B 68, 224518 (2003).
  20. C. Bauer et al., J. Magn. Reson. 58, 442 (1984).
  21. P. Rebentrost and F. K.Wilhelm, Phys. Rev. B 79, 060507 (2009).
  22. S. Safaei et al., Phys. Rev. B 79, 064524 (2009).
  23. H. Jirari and W. Pötz, Phys. Rev. A 72, 013409 (2005).
  24. W. S. Warren, J. Chem. Phys 81, 5437 (1984).
  25. M. D. Bowdrey et al., Phys. Lett. A 294, 258 (2002).
  26. S. L. Patt, J. Magn. Reson. 96, 94 (1992).
  27. N. Khaneja et al., J. Magn. Reson. 172, 296 (2005).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description