This invention relates generally to quantum computing and, more specifically, relates to cross-resonance (CR) gates in quantum computing.
This section is intended to provide a background or context to the invention disclosed below. The description herein may include concepts that could be pursued, but are not necessarily ones that have been previously conceived, implemented or described. Therefore, unless otherwise explicitly indicated herein, what is described in this section is not prior art to the description in this application and is not admitted to be prior art by inclusion in this section. Abbreviations and acronyms used in this document and/or the drawings are defined below, prior to the claims.
The cross-resonance (CR) gate is an entangling gate for superconducting qubits that uses only microwave control (see C. Rigetti and M. Devoret, Phys. Rev. B 81, 134507 (2010); and J. M. Chow, et al., Phys. Rev. Lett. 107, 080502 (2011)) and has been the standard for multi-qubit experiments in superconducting architectures using fixed-frequency transmon qubits (see J. M. Chow, et al., Nature Communications 5, 4015 (2013), arXiv:1311.6330 [quant-ph]; and A. Corcoles, et al., Nature Communications 6 (2015), 10.1038/ncomms7979). Superconducting qubits arranged with shared quantum buses (see J. Majer, et al., Nature 449, 443 (2007)) allow qubit networks to be designed with any desired connectivity. This flexibility of design also translates into a flexibility of control and many choices in entangling gate implementations. The CR gate is one choice of two-qubit gates that uses only microwave control, as opposed to using magnetic flux drives to tune two qubits into a specific resonance condition to entangle, as in the controlled-phase gate (see R. Barends, et al., Nature 508, 500 (2014); and L. DiCarlo, et al., Nature 460, 240 (2009)), or to tune a coupler directly. For the latter, see the following: A. O. Niskanen, et al., Science 316, 723 (2007); R. C. Bialczak, at al., Phys. Rev. Lett. 106, 060501 (2011); Y. Chen, et al., Phys. Rev. Lett. 113, 220502 (2014); and D. C. McKay et al., Phys. Rev. Applied 6, 064007 (2016).
The CR gate requires a small static coupling of the qubit pair that slightly hybridizes the combined system and requires one additional microwave drive. The relatively low overhead of the CR scheme (the additional control line is combined with a single-qubit drive at room temperature) makes it an attractive gate for use in quantum computing architectures based on planar superconducting qubits. Additionally, the CR gate is well-suited to transmon qubits (see J. Koch, et al., Phys. Rev. A 76, 042319 (2007)), which have become the superconducting qubit of choice due to promising long coherence and lifetimes (see A. A. Houck, et al., Phys. Rev. Lett. 101, 080502 (2008); and R. Barends, et al., Phys. Rev. Lett. 111, 080502 (2013)), limited charge noise (see J. M. Gambetta, et al., Phys. Rev. Lett. 109, 240504 (2012)), and high single-qubit gate fidelities (see S. Sheldon, et al., Phys. Rev. A 93, 012301 (2016)). The microwave-only control allows the use of fixed-frequency transmons, further reducing the sources of possible noise (see D. J. Van Harlingen, et al., Phys. Rev. B 70, 064517 (2004)). Given all of these qualities, the CR gate has been a useful tool for application in multi-qubit experiments, including demonstrations of parity measurements required for the surface code (see J. M. Chow, et al., Nature Communications 5, 4015 (2014)).
Despite the appeal of the CR gate, its implementation has been hindered by slow gate times. The CR gate relies on an always-on qubit-qubit coupling, but large couplings can lead to crosstalk between qubits. This leads to a trade-off between fast, high-fidelity two-qubit gates and high-fidelity simultaneous single-qubit gates. As a result, typical CR gates between transmon devices have resulted in gate times >300˜400 ns, with measured fidelities of 94-96% (see A. Corcoles, et al., Nature Communications 6 (2015), 10.1038/ncomms7979; and A. Corcoles, et al., Phys. Rev. A 87, 030301 (2013)).
This section is intended to include examples and is not intended to be limiting.
A method includes generating, in a system comprising a cross-resonance gate having a superconducting control qubit and having a superconducting target qubit coupled through a bus resonator, echo pulses at a first frequency and directing the echo pulses to the control qubit, wherein the first frequency is on resonance with the control qubit. The method also includes generating cross-resonance pulses at a second frequency on resonance with the target qubit and applying the cross-resonance pulses to the control qubit, wherein the generating and applying the cross-resonance pulses induce rotations on the target qubit through an interaction that is mediated by the bus resonator. The method further includes generating cancellation pulses at the second frequency and applying the cancellation pulses to the target qubit.
An apparatus comprises a system comprising a cross-resonance gate having a superconducting control qubit and having a superconducting target qubit coupled through a bus resonator. The apparatus further comprises system controller circuitry comprising one or more processors and one or memories having computer readable code. The one or more processors, in response to retrieving and executing the computer readable code, cause the system to perform operations comprising: generating, in the system, echo pulses at a first frequency and directing the echo pulses to the control qubit, wherein the first frequency is on resonance with the control qubit; generating cross-resonance pulses at a second frequency on resonance with the target qubit and applying the cross-resonance pulses to the control qubit, wherein the generating and applying the cross-resonance pulses induce rotations on the target qubit through an interaction that is mediated by the bus resonator; and generating cancellation pulses at the second frequency and applying the cancellation pulses to the target qubit.
Another method comprises measuring, in a system comprising a cross-resonance gate having a superconducting control qubit and a superconducting target qubit coupled through a bus resonator, a first set of Hamiltonian tomographies, where cross-resonance pulses at a frequency on resonance with the target qubit are applied to the control qubit, and where Hamiltonian tomographies are measured at specific points as the cross-resonance pulses are swept from a beginning phase to an ending phase at a specific amplitude. The method also comprises determining, using output from the first set of first Hamilton tomographies, a phase to be used for cancellation pulses. The method further comprises measuring a second set of Hamiltonian tomographies by applying cross-resonance pulses at the frequency to the control qubit and applying cancellation pulses at the frequency and at the determined phase to the target qubit, where the cross-resonance pulses have a specific amplitude and a specific phase, where Hamiltonian tomographies are measured at specific points as an amplitude of the cancellation pulses are swept from a beginning amplitude to an ending amplitude. The method additionally comprises determining, using output from the second set of second Hamilton tomographies, an amplitude to be used for the cancellation pulses, and outputting the determined phase and amplitude for the cancellation pulses, to be used to configure the system for subsequent operations of the cross-resonance gate.
In the attached Drawing Figures:
The word “exemplary” is used herein to mean “serving as an example, instance, or illustration.” Any embodiment described herein as “exemplary” is not necessarily to be construed as preferred or advantageous over other embodiments. All of the embodiments described in this Detailed Description are exemplary embodiments provided to enable persons skilled in the art to make or use the invention and not to limit the scope of the invention which is defined by the claims.
Here we describe improvements to the CR gate through a careful Hamiltonian analysis and novel tune-up procedure that reduce the gate time by a factor of two with corresponding fidelities over 99% (as described below, e.g., in reference to
The system tested and used comprised two fixed-frequency transmon qubits coupled by a bus resonator as shown in
In
The qubit frequencies are 4.914 GHz (target qubit 120) and 5.114 GHz (control qubit 110) with an harmonicities of −330 MHz for both, and the bus frequency is 6.31 GHz. The coupling between the two qubits is estimated to be J/2π=3.8 MHz for these parameters. The single-qubit gate fidelities measured with simultaneous randomized benchmarking (see J. M. Gambetta, et al., Phys. Rev. Lett. 109, 240504 (2012)) are 0.9991±0.0002 for the target and 0.9992±0.0002 for control. We characterize the two-qubit gate fidelities in the system by using interleaved randomized benchmarking (RB) (see the following: E. Magesan, et al., Phys. Rev. Lett. 109, 080505 (2012); and J. P. Gaebler, et al., Phys. Rev. Lett. 108, 260503 (2012)). For this measurement, we first find the average fidelity per Clifford in the two-qubit system using standard randomized benchmarking (see E. Magesan, et al., Phys. Rev. A 85, 042311 (2012)), and then repeat the measurement interleaving the CR gate between random Cliffords in the sequence. As is known, fidelity is a measure of how close a gate is to a desired gate, and the higher the fidelity, the closer the implemented gate is to the desired gate. The fidelities are extracted from exponential fits to the average over 35 random sequences each with a total length of 100 Clifford gates. By applying the theoretical Hamiltonian understanding of the teachings below, we were able to benchmark results showing a CR gate fidelity, f=0.991±0.002 for a 160 ns gate. This gate time includes 20 ns added by the single qubit echo, which is buffered by two 10 ns delays. The CR pulses are rounded square pulses with Gaussian rise times of 15 ns.
The critical experimental technique for the improvement of the CR gate is an active cancellation pulse on the target qubit drive to eliminate unwanted interactions of the CR drive Hamiltonian. A circuit to implement the RF signals to the control qubit 110 and the target qubit 120 is illustrated in
The qubit control 180 can configure this circuit as follows. In particular, one may program the system 100 to cause the qubit control 180-1 (the computer readable code) to configure the circuit in
Turning to
It is noted that the CR pulses 195 and cancellation pulses 197 have the same analytical form, the same duration, and the same rise/fall time. However, the amplitude and phase of the two pulses will be different (e.g., the point of the calibration described below is to determine the amplitude and phase of the cancellation pulse). The pulse amplitude may be thought of as indicating an envelope around the oscillating carrier tone at the qubit frequencies. Another way of looking at this is that positive amplitude indicates rotation about the +x-axis and negative amplitude indicates rotation about the −x-axis (although it's not necessary the x-axis, the exact axis depends on the phase of the pulse). Equivalently, the negative amplitude pulse is exactly π out of phase with the positive amplitude pulse.
Turning to
The addition of the cancellation drive on the target is the consequence of measurements and simulation of the CR Hamiltonian indicating single-qubit errors that are not fully refocused by the standard CR gate (e.g. the echoed sequence in
where X is the eigenvector matrix of H and is assumed to be non-singular, and XBD is the block-diagonal matrix of X.
In our implementation, the different blocks correspond to the different states of the control qubit and off-resonant (higher energy) qubit subspaces. The model predicts ZX and IX components of similar magnitude, negligible IZ and ZZ contributions, and a large ZI term arising from a Stark shift of the control qubit from off-resonant driving. The complete CR Hamiltonian has the following structure:
where A and B are distinct unitary rotations on the target qubit. Note that in this context the operators A and B are used to illustrate the structure of the CR Hamiltonian, and a technique for determining these is described below.
Much of the rest of this disclosure is described in reference to
Motivated by this understanding of the CR Hamiltonian, we have developed a protocol for experimentally measuring the CR Hamiltonian that allows us to determine the real error terms in the gate. See block 410 of
1. Set the phase of the CR pulse to zero (and the amplitude of the CR pulse to a specific amplitude).
2. Run Hamiltonian tomography and extract the Hamiltonian for this phase.
3. Increment the phase of the CR pulse to 2π/N.
4. Loop back to #2 N times, until the Hamiltonian has been measured for N CR phases between 0 and 2π.
5. From the Hamiltonian versus phase data, extract the correct phase for the CR pulse and the correct phase for the cancellation pulse (as described in more detail below).
A similar procedure is used for sweeping the amplitude.
The Hamilton tomography 405 is shown in more detail in
When this quantity goes to zero, the two-qubits are maximally entangled (unless the dynamics are completely mixed, but because ∥{right arrow over (R)}∥ oscillates between zero and one we believe this is not the case). We use ∥{right arrow over (R)}∥ to estimate the gate length required to perform the entangling gate. We fit the set of Rabi oscillations that correspond to either state (|0> or |1>) of the control qubit separately with a Bloch equation model function:
{right arrow over (r)}(t)=eGt{right arrow over (r)}(0), (4)
with the matrix G defined as the following:
Here Δ is the control drive detuning, and Ωx,y is the Rabi drive amplitude along {x, y}. Additionally, {right arrow over (r)}(t) is the vector composed of the measured expectation values as a function of the length of the applied Rabi drive, (<X(t)>,<Y(t)>,<Z(t)>), and {right arrow over (r)}(0) is {right arrow over (r)}(t) at t=0. We find two generators, G{0,1}, corresponding to the control qubit in either |0> or |1>, respectively, characterized by the terms Ωx{0,1}, Ωy{0,1}, and Δ{0,1}. See block 480, where #1 (block 465) and #2 (block 470) are fit separately to the matrix G. This gives G0 and G1. From these terms, we derive the CR drive Hamiltonian in terms of the six possible parameters (see block 490): IX, IY, IZ, ZX, ZY, ZZ. For example, IX=(Ωx0+Ωx1)/2 and ZX=(Ωx0−Ωx1)/2. By measuring the Bloch vectors as a function of time dependent on the state of the control, we can determine the operators A and B of Eq. 2.
We can vary the CR drive amplitude, perform the CR-Rabi experiments, and extract the CR Hamiltonian using this Bloch equation fitting method. This provides a measurement of the CR Hamiltonian as a function of the two-qubit drive amplitude, as shown in
There is an unexpected feature in the experiment as there is also an IY term present when the CR phase is set as above. We attribute this phase difference between conditional and single-qubit terms to classical crosstalk. Although such crosstalk has little effect on simultaneous RB since it is off-resonance, since the CR drive is applied at the frequency of the target qubit, this crosstalk has a significant impact on the two-qubit gate calibration. The standard CR gate is performed using an echo to refocus the parameters IX, ZZ, and ZI, depicted in
The gate calibration procedure is based on these Hamiltonian measurements. Ultimately, the goal is to tune up a ZX90, which is a generator of a controlled-NOT (CNOT) with single-qubit Clifford rotations. The first step is to find the Hamiltonian parameters as a function of the phase of the CR drive (see
With the CR phase fixed at ϕ0, the two-qubit drive produces nonzero IX and IY components. We find the phase, ϕ1, at which the single-qubit component, IY, is zero. The correct phase for the cancellation pulse is ϕ=ϕ0−ϕ1, at which phase the single-qubit drive on the target matches tan−1(IY/IX) for the two-qubit drive. See block 430 of
The second step is to set the cancellation pulse to the correct amplitude for canceling cos(ϕ)IX+sin(ϕ)IY. See block 430 of
In block 460 of
The CR-Rabi oscillations on the target qubit with the fully calibrated cancellation pulse (
Plotting the same data on the Bloch sphere, as in
While the measured fidelity is higher than previously reported fidelities for CR two-qubit gates, it is still not yet limited by coherence. The limit placed on the fidelity by T1 (38±2 μs/41±2 μs for the control/target) and T2 (50±4 μs/61±6 μs) for the two-qubits is 0.996. There is evidence from studies of single-qubit gates that drive-activated dephasing may be a limiting factor for the error rates of gates at similar drive powers (see the following: S. Sheldon, et al., Phys. Rev. A 93, 012301 (2016); and H. Ball et al., ArXiv e-prints (2016), arXiv:1602.04551 [quant-ph]). Even so, it appears from the CR Rabi data that coherent errors have not been fully eliminated in the CR gate.
While it is clear that some coherent errors remain on the cross-resonance gate, the inclusion of an active cancellation tone has produced a dramatic improvement in the two-qubit gate fidelity. The calibration procedure that we have developed has provided insight into the full CR Hamiltonian, and revealed that a single-qubit phase shift due to crosstalk (either classical or quantum) is a significant source of error in the echoed CR gate. Future work might focus on shortening the gate times further by tuning up a CR gate without an echo and developing faster and more robust calibration procedures.
These improvements to the microwave-driven CR gate are encouraging for a quantum computing architecture built out of fixed-frequency qubits. With greater knowledge of the CR Hamiltonian, we are confident that further improvements can be made to reduce the two-qubit gate time and increase fidelity.
Note that the previously described method of Hamiltonian tomography is applicable to any system with a Hamiltonian with the same form as Eq. 2. In addition, due to the symmetry of , this method scales efficiently for an n-qubit system, since there are n(n−1)/2 different control/target qubit pairs and each pair requires six Rabi measurements described above.
As an example, we can extend the Hamiltonian tomography method to measure the effects of crosstalk from neighboring qubits in multi-qubit system. For instance, in a system with three qubits, there is a control qubit, a target qubit, and a spectator qubit. The spectator qubit should do nothing during the two qubit gate that is applied between the control and target qubits. Since all the qubits are coupled, however, the spectator qubit can cause additional Z errors on the target qubit. We can measure these errors by performing Hamiltonian tomography with an additional step of measuring the Hamiltonian when the spectator qubit is in both |0> and |1> states. In this case, there are 12 total Rabi experiments: three projections of the target qubit for the cases of (1) control in |0> and spectator in |0>, (2) control in |1> and spectator in |0>, (3) control in |0> and spectator in |1>, and (4) control in |1> and spectator in |1>.
The present invention may be a system, a method, and/or a computer program product at any possible technical detail level of integration. The computer program product may include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.
The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A non-exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or raised structures in a groove having instructions recorded thereon, and any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.
Computer readable program instructions described herein can be downloaded to respective computing/processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and/or a wireless network. The network may comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. A network adapter card or network interface in each computing/processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing/processing device.
Computer readable program instructions for carrying out operations of the present invention may be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, configuration data for integrated circuitry, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++, or the like, and procedural programming languages, such as the “C” programming language or similar programming languages. The computer readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate arrays (FPGA), or programmable logic arrays (PLA) may execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present invention.
Aspects of the present invention are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer readable program instructions.
These computer readable program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function/act specified in the flowchart and/or block diagram block or blocks.
The computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computer, other programmable apparatus, or other device implement the functions/acts specified in the flowchart and/or block diagram block or blocks.
The flowchart and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical function(s). In some alternative implementations, the functions noted in the blocks may occur out of the order noted in the Figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems that perform the specified functions or acts or carry out combinations of special purpose hardware and computer instructions.
The following abbreviations that may be found in the specification and/or the drawing figures are defined as follows:
AWG arbitrary waveform generator
a tan arctangent (tan−1(•))
a.u. arbitrary units
BD block-diagonal
CNOT controlled-NOT
CR cross-resonance
CRC computer-readable code
DRAG derivative removal via adiabatic gate
GHz gigahertz
LO local oscillator
MHz megahertz
ns nanoseconds
RB randomized benchmarking
RF radio frequency
ssb single sideband
This application claims the benefit of U.S. Provisional Application No. 62/308,137, filed on Mar. 14, 2016, the disclosure of which is hereby incorporated by reference in its entirety.
This invention was made with Government support under contract W911NF-14-1-0124 awarded by the Army Research Office. The Government has certain rights in this invention.
Number | Date | Country | |
---|---|---|---|
62308137 | Mar 2016 | US |