Implementing a strand of a scalable

Full text

Turn on search term navigation

ARTICLE

Received 11 Feb 2014 | Accepted 30 Apr 2014 | Published 24 Jun 2014

DOI: 10.1038/ncomms5015

Implementing a strand of a scalable fault-tolerant quantum computing fabric

Jerry M. Chow1, Jay M. Gambetta1, Easwar Magesan1, David W. Abraham1, Andrew W. Cross1, B.R. Johnson2, Nicholas A. Masluk1, Colm A. Ryan2, John A. Smolin1, Srikanth J. Srinivasan1 & M. Steffen1

With favourable error thresholds and requiring only nearest-neighbour interactions on a lattice, the surface code is an error-correcting code that has garnered considerable attention. At the heart of this code is the ability to perform a low-weight parity measurement of local code qubits. Here we demonstrate high-delity parity detection of two code qubits via measurement of a third syndrome qubit. With high-delity gates, we generate entanglement distributed across three superconducting qubits in a lattice where each code qubit is coupled to two bus resonators. Via high-delity measurement of the syndrome qubit, we deterministically entangle the code qubits in either an even or odd parity Bell state, conditioned on the syndrome qubit state. Finally, to fully characterize this parity readout, we develop a measurement tomography protocol. The lattice presented naturally extends to larger networks of qubits, outlining a path towards fault-tolerant quantum computing.

1 IBM T.J. Watson Research Center, Yorktown Heights, New York 10598, USA. 2 Raytheon BBN Technologies, Cambridge, Massachusetts 02138, USA. Correspondence and requests for materials should be addressed to J.M.C. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5015

Quantum error correction is an essential step towards realizing scalable quantum computers. Theoretically, it is possible to achieve arbitrarily long protection of quantum

information from corruption due to decoherence or imperfect controls, so long as the error rate is below a threshold value1,2. However, most of the initial fault-tolerant quantum computing proposals were purely theoretical studies that would be impractical to implement in a physical system. Knills C4/C6 architecture3 showed that it was possible to have a pseudo-threshold as high as 3% but with very long-range interactions. The two-dimensional surface code (SC) is a topological error-correcting code4,5 that achieves high error thresholds B1% in a two-dimensional lattice of qubits supporting nearest-neighbour interactions6,7.

With ongoing improvements to coherence times811, and with gate12,13 and readout delities1416 at or approaching relevant threshold values, superconducting qubits are prime candidates for scaling up towards a fault-tolerant architecture. Our scheme for building a network of superconducting qubits employs microwave resonators as the links, as it has been shown that resonators can be used as quantum buses to mediate interactions between qubits17, and that multiple resonators can be coupled to a single qubit18. In the future, larger superconducting qubit systems and networks will also leverage circuit integration techniques that come along with a solid-state architecture.

The original SC layout comprises qubits arranged in a square lattice5 with the qubits in the lattice coming in two distinct avours, either code qubits, which carry logical information, or syndrome qubits, which are used to measure stabilizer operators of surrounding code qubits. Performing a round of error correction in the SC consists of applying controlled-NOT (CNOT) gates to map the parity of the surrounding code qubits into the states of the syndrome qubits, which are then measured to determine the plaquette parities. Using superconducting resonators as the links of the lattice, it is possible to construct a SC architecture using superconducting qubits at each vertex along with the ability to couple a single qubit to four resonators. As it is also important to be able to read out and address the qubits individually, this may require an additional fth resonator per qubit. A layout proposed by DiVincenzo19 allows for the number of resonators per qubit to be reduced to two (three with readout). However, each code and syndrome qubit in this scheme is replaced by four qubits, which results in a fourfold increase in the number of two-qubit gates.

Here we introduce an optimal layout for implementing the SC on a two-dimensional lattice that does not carry additional overhead as in the original DiVincenzo proposal. In this work, we demonstrate a subsection of this optimized lattice using three transmon qubits, in which two independent outer code qubits are joined to a central syndrome qubit via two independent bus resonators. By combining all-microwave high-delity single- and two-qubit nearest-neighbour entangling gates, with high-delity measurement of the syndrome qubit, we deterministically measure the parity of the code qubits. This protocol is tested in the case of a superposition of the code qubits, and we are able to detect the outer code qubits in either an even or odd parity Bell state, conditioned on the syndrome state. As a full characterization of this parity detection, we implement a measurement tomography protocol to obtain a delity metric (90% for odd and 91% for even). Our results reveal a straightforward path for expanding superconducting circuits towards larger networks for the SC and eventually a primitive logical qubit implementation.

ResultsModied lattice for SC with superconducting qubits. We introduce a direct one-to-one mapping of the SC to a

two-dimensional lattice (Fig. 1) that does not carry the additional overhead of qubits and gates as in the original proposal by DiVincenzo19. Rather than requiring 12 qubits to demonstrate a half-plaquette strand operation on the surface fabric, this arrangement allows such an achievement with just 3 qubits (details about the mapping are given in the Methods). Syndrome qubits can be used as either X-parity checks (red circles) or Z-parity checks (green circles) of four surrounding code qubits (blue circles). The ve-qubit block consisting of X (Z)-syndrome and four surrounding code qubits denes a unit cell X (Z)-plaquette.

An advantage of our lattice architecture is that it is commensurate with experiments already demonstrated, where a single qubit can be coupled to two separate resonators18, a single bus has been used to couple up to three qubits20, and an independent readout can be used to measure a single qubit that has been entangled with a separate qubit via a bus resonator21. The experiments presented in this article are performed effectively on a half-plaquette subsection, consisting of three qubits (all with individual readout resonators) and two bus resonators, depicted by the pink square of Fig. 1, where we demonstrate all necessary gate operations and measurements that constitute the SC protocol.

Physical subsection and control characterization. The half-plaquette device (Fig. 2a) contains three single-junction transmon qubits connected by two coplanar waveguide (CPW) resonators serving as the buses, and each qubit is coupled to its own separate CPW resonator for independent readout and control. Here we label the code qubits Q1 and Q3, and the middle qubit Q2 serves

Z-syndrome qubit

X-syndrome qubit

Code qubit

R12

R23

Figure 1 | Implemented two-dimensional SC layout. Lattice of nearest-neighbour qubits for realizing the SC. The lattice contains three avours of qubits, code (blue circles), X-syndrome (red circles) and Z-syndrome (green circles) qubits. The lattice is composed of two types of plaquettes, an X-plaquette (square enclosing a red X-syndrome qubit) and a Z-plaquette (square enclosing a green Z-syndrome qubit). Critical to the SC is performing CNOT operations between code qubits and their neighbouring syndrome qubits, followed by measuring the Z- and X-syndrome qubits to determine the parity of the code qubits. Bus resonators that link qubits are denoted by grey bars and each qubit need only connect to two bus resonators instead of the required four in the strict square lattice approach. The smallest logical qubit is shown in the shaded blue area, made up of four plaquettes. The half-plaquette device studied here is enclosed in the pink square and comprises three qubits (Q1, Q2, Q3) and two bus resonators R12 and R23. Note that Q3 is also a code qubit, but coloured in teal for clarity to be distinguished from Q1.

2 NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5015 ARTICLE

JPA

[afii9853]p

[afii9853]M3

|[afii9852]

HEMT

R23

[afii9853]3

[afii9853]2

[afii9853]M2

|[afii9852]

[afii9853]2

HEMT

X()

n90

Z+90

X90

HEMT

100 m

[afii9853]2 [afii9853]1 [afii9853]M1

R12 Q1

X()

n90

X+90

X90

Figure 2 | Half-plaquette device schematic and parity check quantum circuits. (a) The optical image of the half-plaquette device shows in false colour all the different components of the device: three qubits, Q1 (blue), Q2 (green) and Q3 (teal), each with individual readout resonators, and two bus resonators (maroon) R12 and R23. Each transmon qubit (zoom view inset) is independently addressed via its corresponding readout resonator, with single-qubit gates applied on resonance with each qubit at oi, iA[1, 2, 3] and readout performed at the measurement frequencies oMi. Whereas

Q1 and Q3 readout signals are only amplied through HEMTs, the Q2 readout is reected off a JPA stage rst before going on to a HEMT. Two-qubit gates are performed in the cross-resonance scheme, applying o2 on both control qubits, Q1 and Q3. (b) The PCP for qubit Q1 and Q3 where the Z-parity operator PZZ is applied, giving a single classical bit of information b2 (double lines indicate classical channel). (c) The quantum circuit which implements the Z-parity check consists of a pair of CNOT gates from the code qubits (Q1 and Q3) to the syndrome (Q2) followed by a measurement M2, which gives the classical bit b2. (d) The CNOT can be decomposed into the ZX90 gate and single-qubit rotations. Using the cross-resonance microwave interaction, we have at our disposal the gate combination boxed in dashed red, composed of a ZX90 followed by a NOT (or X) gate on the control qubit. The n in the depiction of the ZX90 gate can be either 0 or 1, indicating the state-dependent rotation.

as the syndrome. A simplied control and readout diagram is also shown in Fig. 2a, with all microwave sources for single- and two-qubit gates and independent readout indicated. The syndrome Q2 is read out with the assistance of a Josephson parametric amplier (JPA)16,22 for high-delity single-shot state discrimination, which is a critical component for demonstrating the SC parity check protocol. All device parameters, relevant coherence times and a complete schematic are given in the Methods. All single-qubit controls are performed in 40 ns and characterized to 499.7% average gate delity via simultaneous randomized benchmarking23 (see Methods).

At the crux of the SC protocol is the detection of the parity of the code qubits. Figure 2b shows the parity check protocol (PCP) in a circuit depiction on an arbitrary state c

j i of Q1 and Q3, with

the parity being indicated through the classical detection of an indication bit, b2. The PCP is realized in a system of three qubits via the quantum circuit shown in Fig. 2c, where the parity of the Q1 and Q3 state c

j i is mapped on to the syndrome Q2 via two

CNOT gates, between Q1 and Q2, Q3 and Q2, and then subsequently the classical indication bit b2 is obtained via a quantum measurement of Q2.

High-delity CNOT entangling gates are critical for the PCP. To realize these CNOTs in our device, we implement entangling ZX90 gates using the cross-resonance interaction24,25. By driving

Q1 (Q3) at the Q2 transition frequency o2, we use a simple decoupling sequence12 to implement the two-qubit Clifford generator ZX90 gate between Q2 and Q1 (Q3). The ZX90 gate is equivalent to a CNOT up to single-qubit rotations, as illustrated in Fig. 2d, and thus can be used interchangeably in the PCP.

Both pairs of ZX90 entangling gates are characterized with two-qubit Clifford randomized benchmarking. With a total gate time of 350 ns, the ZX90 gate between Q2 and Q1 (Q3) is experimentally shown to have a gate delity of 0.9620.002

(0.9570.001). The delities are consistent with the gate and coherence times in the system. Further details about the gate tune-up, calibration and timescales can be found in the Methods.

To determine the collective state of all qubits in the system, we can perform independent readouts of each qubit through the individually coupled resonators. The syndrome qubit, Q2, is read out dispersively with the JPA, pumped 4 MHz from the

readout frequency (oM2 2p 6.584 GHz) and an optimized

single-shot assignment delity (with no preparation corrections, and dened in Methods) of 91% is achieved, although uctuations on the order of 23% are observed. The code qubits are read out using the high-powered Josephson nonlinearity of the readout cavities26. Further detail and parameters about the readouts are given in the Methods.

Distributed entanglement and parity measurement. To observe the action of the gate portion of the PCP, in which CNOT (in our case ZX90) gates are performed between the syndrome and the code qubits, it is insightful to perform tomographic reconstruction of the complete three-qubit system. State tomography in our system is achieved by correlating individual single shots of all three individual readouts27, Mi, iA[1,3]. Figure 3 shows reconstructed three-qubit Pauli state vectors for the entanglement processes necessary for the PCP. In Fig. 3a, a ZX90 entangling operation between the code qubit Q3 and the syndrome Q2 is implemented giving a state delity F

state

0.9490.002 (SDP)0.9540.002 (raw), where SDP refers to a semi-denite program reconstruction of the state and raw reects unconstrained inversion28. The difference between the SDP and raw estimates exceeds the statistical error estimated from the readout signal-to-noise and suggests that the main source of error in our experiment is systematic. We estimate that our systematic

NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5015

X+90

I,X X,Y

Q1 |0

0.5

X+90

I,X X,Y

0.5

X()

I,X X,Y

X+90

I,X X,Y

0.5

X+90

I,X X,Y

0.5

X()

I,X X,Y

X+90

I,X X,Y

0.5

X+90

I,X X,Y

0.5

X()

I,X X,Y

State tomography

I I I

I I X

I I Y

I I Z

I X

I Y

I Z

I I

I X X

I Y X

I Z X

I X Y

I Y Y

I Z Y

I X

I Y

I Z Z

I X

I Y

I Z

X X

Z X

Y X

Y Y

Z Y

Z Z

X X X

Y X X

Z X X

X Y X

Y Y X

Z Y X

Z X

Z Z X

X X Y

Y X Y

Z X Y

X Y Y

Y Y Y

Z Y Y

Z Y

Z Z Y

X X

Y X

Z X

Y Y

Z Y

Z Z

Z Z Z

Pauli operator

Figure 3 | Three-qubit state tomography via correlated single-shot readout. Building towards the action of the PCP on a superposition state of the code qubits (c), we also show the operation of entangling the syndrome Q2 with each of the data qubits Q1 (a) and Q3 (b). Reconstructed three-qubit state represented in Pauli state vector form (black lled bars, experiment, white bars in background, ideal) after entangling syndrome Q2 with, a, code qubit Q3 via a ZX90 two-qubit gate (entangled two-qubit state with delity F

stateB0.95), (b) code qubit Q1 via a ZX90 two-qubit gate

(entangled two-qubit state with delity F

stateB0.94). In the shown Pauli vector plots, the blue-, pink- and purple-shaded regions signify single-, two- and three-qubit Pauli operators, respectively.

stateB0.95), (c) both code qubits Q1 and Q3 via ZX90 between Q1 and Q2, and Q3 and Q2 simultaneously (which comprise performing the gate portion of the PCP) giving the maximally entangled three-qubit GHZ state (F

errors are of order B23%. A rigorous error bound on delity in the presence of systematic noise is still the subject of future study29. In Fig. 3b, a ZX90 entangling operation between Q1 and

Q2 is implemented (F

state

0.9510.002 (SDP) 0.9530.003 (raw)) and nally, in Fig. 3c both ZX90 gates are applied simultaneously, generating a maximally entangled GHZ state of all three qubits (F

state

0.9350.002 (SDP) 0.9420.003 (raw)).

We thus show the ability to distribute entanglement across the full network, rst between nearest-neighbour qubits Q3 and Q2 or Q1 and Q2, and then across the entire system with the GHZ state, spanning both bus resonators.

With the demonstrated high-delity ZX90 gate primitives for the PCP, the next step is to observe the single-shot readout of the syndrome Q2 to signal the parity of the code qubits Q1 and Q3. We begin by obtaining a syndrome readout calibration histogram (a typical one is shown in Fig. 4a). Then, we prepare the simple computational basis states, 00

j i, 01

j i, 10

j i and 11

j i, as inputs for

the code qubits and observe the proper parity assignment via the PCP, shown via the M2 histograms in Fig. 4b. By thresholding the measurement outcomes of M2 based on the readout calibration traces, we reconstruct the state of Q1 and Q3 conditioned on M2 using standard quantum state tomography techniques (see Methods). In the case of the four computational states, we obtain delity of F

SDP

error-free situation, these delities should not exceed the assignment delity; however, due to systematic errors it is possible for these to deviate. The observed deviation between the assignment and state delity is within our estimate of systematic error. Nonetheless, the ability to generate these entangled states of non-nearest-neighbour qubits Q1 and Q3 in two separate bus resonators is a crucial element for scaling up towards larger quantum networks.

Finally, to characterize the complete ideal projective nature of the PCP, we perform measurement tomography. This is accomplished via quantum process tomography of the code qubits, for which further details are given in the Methods. Conditioned on the measurement of the syndrome M2, we obtain the two maps for the odd and even parity projection operators shown in Fig. 4f,h (ideal maps are shown in Fig. 4i,k). To quantify the performance of the PCP, we introduce a measurement delity metric that takes into account the full quantum dynamics of the measurement, including projection and back-action (see Methods for more details). We obtain a measurement delity of Foddmeas

0:904 0:002 and Fevenmeas 0:912 0:003. The loss in

measurement delity corresponds mostly to the 91% assignment delity of M2, best illustrated by the unconditional map (shown in Fig. 4g) having a measurement delity F

meas

0.975,

0.989, 0.999, 0.905, respectively.

A more complete stress test of the PCP is to observe its

function on the maximal superposition state of the code qubits Q1 and Q3. The gate protocol now mimics that of the GHZ state generation from Fig. 3c, and over repeated state preparations and measurements of the syndrome, M2, we obtain a bi-modal histogram, indicating that instances of both parities exist (Fig. 4c). We observe the probabilistic entanglement of either the odd or even Bell states codd

j i 01

j i 10

j i

0.968.

It is important to note that all the gates used are calibrated and run to achieve these results without any Hamiltonian corrections for either single-qubit or two-qubit errors. An X-parity check is a simple extension through the appropriate application of single-qubit Hadamard pulses on the code qubits.

DiscussionThe experiment described implements a subsection of the SC fault-tolerant architecture. By combining high-coherence transmon qubits, high-delity nearest-neighbour two-qubit gates and high-delity quantum non-demolition single-shot readout, we use a syndrome qubit to determine the parity of its neighbouring

0.984, 0.987, 0.989, 0.909 and F raw

p or ceven

j i

p conditioned on M2. For these conditioned entangled states, we nd state delities of F

j i 11

j i

odd

0.891 (raw)

0.891 (SDP) and F

even

0.970 (raw) 0.948 (SDP). In an

4 NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5015 ARTICLE

0.5

0 0.5 1

Experiment Ideal

0.4

|Q2 = |0

|Q2 = |1

Fit

Threshold

|Q1Q3 =

|01 |10

|11 |00

0.2

M2voltage (mV)

0.2

0.4

Z Y

X Y

Y X Y

Output Pauli operator

...

X I Z

Y Y

I X X

Pauli operator

0.02

0.04

0.02

0.04

0.02

0.04

...

Count fraction

Input Pauli operator

Figure 4 | Parity check experiment and non-nearest-neighbour pair entanglement generation. (a) Readout histograms of M2 having prepared Q2 in the state 0

j i and 1

j i averaged over all basis states for Q1 and Q3. (b) Readout histograms of M2 after applying the PCP when Q1 and Q3 are prepared

in the four different basis states, 00

j i (black), 01

j i (teal), 10

j i (blue), 11

j i (purple). (c) Readout histogram of measurement M2 after applying the PCP

when Q1 and Q3 are prepared in an equal superposition state. (d,e) Pauli state vectors of Q1 and Q3 conditioned on the single-shot measurement of Q2. In the case of Q2 in 0

j i 1

j i

, state tomography conrms the odd (even) parity Bell state, codd

j i 01

j i 10

j i

p ceven

j i 00

j i 11

j i

with

delity 0.89 (0.95). (fk) Pauli transfer matrices of the Z-parity check measurement operation. Top (bottom) corresponds to the odd (even) projection and has a measurement delity of 0.90 (0.91). The central measurement operation corresponds to the unconditional map and has a delity of 0.968 with a map that is completely dephasing in the even and odd parity basis.

qubits. With this device, we demonstrate the versatility of superconducting qubits in the extended quantum bus architecture for application towards a larger fault-tolerant quantum computing device. Looking ahead, direct extensions to simple error detection and error correction demonstrations will be feasible with existing integration techniques and coherence levels. Overcoming integrated circuit engineering hurdles while preserving long coherence times should pave the way for larger surfaces for quantum error correction.

During the completion of this manuscript, we became aware of similar work by Saira et al.33

Methods

Mapping to SC. The SC requires CNOT gates to be performed between the code and syndrome qubits. To optimize the delity of these gates, it is important to use a connective element (connector) to mediate a strong interaction between these qubits. For superconducting transmon qubits, the connections are typically achieved by capacitively coupling two (or more) qubits to a resonator. Since the SC is a two-dimensional geometry, connectors should not cross.

The most straightforward layout of the SC is to have the qubits arranged ona square lattice, each connected to its nearest neighbours (see Fig. 5a). This arrangement has considerable overhead, in that each qubit is attached to four connectors. For a realization in a transmon-based architecture it will require each qubit having at least four capacitive couplers (not including readout and control lines). As a result, a signicant amount of the qubit charging energy will come from the presence of the couplers.

Previously, DiVincenzo19 has shown that it is possible to have a lattice that realizes the SC with each qubit attached only to two connectors. The DiVincenzo lattice is shown in Fig. 6. In this lattice, each small square contains four qubits and is used to represent a single code or syndrome qubit. It is apparent that this arrangement has connectivity between the small squares that is equivalent to the standard square lattice. However, the mapping is inefcient in that it requires four qubits to represent a single qubit in the original SC and furthermore it requires as many as four CNOT operations spanning four connections (resonators) for the furthest reaches of adjacent squares to communicate.

Figure 5b shows three additional arrangements of connections in the square lattice. Connectors are located in every other square of the checkerboard, shown as pink, and each pink square represents either four, two or one connector. In the rst case, the pink square is broken up into four connectors, each attached to two qubits and each qubit to four connectors, recovering Fig. 5a. We call this the [2,4] arrangement. The second arrangement is to divide the pink regions into two, with each connector joined to three qubits and each qubit attached to three couplers ([3,3] arrangement).

Finally, the connector can be the entire pink square connecting four qubits with each qubit connecting only to two connectors ([4,2] arrangement). This nal arrangement can easily been seen to achieve the desired simultaneous goals of having each qubit attached to just two connections while retaining the connectivity of the square lattice. As shown in Fig. 1, the connector can be realized by using a resonator that connects to four qubits. The alternating horizontal and vertical lines inside the pink squares correspond precisely to the bus resonators.

Two nal points are worth discussing. First, note that the [4,2] arrangement is actually more connected than the square lattice, that is, diagonal neighbouring qubits can be coupled. This extra connectivity may be of some advantage, although it comes at the cost of reducing the isolation between these qubits. This is an unavoidable consequence of having more than two qubits per connector. Lastly, DiVincenzos lattice is a simple distortion, a skewing, of the [4,2] lattice. Our key observation can be thought of as seeing that his lattice can be directly mapped to the square lattice without having four qubits being used to represent just one qubit from the SC.

Device fabrication. The device is fabricated on a 720-mm thick silicon substrate. All superconducting CPW resonators are dened via optical lithography and subtractive reactive ion etching of a sputtered niobium lm (200 nm thick). The three single-junction transmon qubits are patterned using electron-beam lithography, followed by double-angle deposition of aluminium, with layer thicknesses of 35 and 85 nm. Lift-off process is used to form the nal junction structure.

Device parameters. The three transmon qubits (iA[1,3]) have transition frequencies {oi}/2p {5.0388, 5.0080, 5.2286} GHz, with readout resonators at

{oRi}/2p {6.698, 6.585, 6.695} GHz, relaxation times {T1(i)} {24, 29, 20} ms,

Techo2in o

32; 25; 18

f g ms. The bus resonators are unmeasured but R12 (R23) is designed to resonate at 8 (8.5) GHz. The dispersive cavity shifts of the readout

NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5015

Z-syndrome qubit

X-syndrome qubit

Code qubit

Figure 6 | Original skew-square lattice of DiVincenzo. Each code and syndrome qubit is represented by the small coloured squares, each containing four physical qubits. The black lines represent resonators (couplers).

[2.4] [3.3] [4.2]

Figure 5 | Square lattice with connectors. (a) The conventional square lattice. Each qubit is connected to its four nearest neighbours. The lattice contains three avours of qubits, code (blue circles), X-syndrome (red circles) and Z-syndrome (green circles) qubits. (b) Three kinds of connectivity: the pink squares shows a connector that is used to allow the application of gates required for the SC. This square can be made from four connectors, two connectors or a single connector.

resonators are measured to be {wi}/p { 2.0, 2.0, 2.3} MHz and the readout

resonators have line widths {kRi}/2p {443, 976, 793} kHz. All qubits have

measured anharmonicities of B 340 MHz. From the above, we calculate

coupling strengths {|gi|}/2p {70, 67, 67} MHz to the readout resonators, which is

consistent with electromagnetic simulations.

Experimental set-up. The half-plaquette device is cooled to 15 mK in an Oxford Triton dilution refrigerator. A full schematic of the wiring and experimental control hardware is depicted in Fig. 7. Each qubit has its own dedicated readout line with an associated set of isolators and Caltech HEMT (noise temp B6 K)

ampliers. Q2 is unique in that its readout signal is reected off of a UC Berkeley JPA before going on to the isolator and high-electron mobility transistor (HEMT) chain. The device is housed in a light-tight Ammuneal cryoperm-shield that is coated throughout with a layer of lossy eccosorb (Emerson & Cuming CR-124). Besides explicit cryogenic attenuators at the different stages of the cryostat, all qubits are also attenuated at the lowest temperature stage with in-house eccosorb coaxial lters.

Outside the cryostat, all microwave qubit control signals are generated via vector modulation combining off-the-shelf electronics. The microwave readout signals are pulse modulated using Arbitrary Pulse Sequencers built by Raytheon BBN Technologies. The readout signals are processed via two Alazartech ATS9870 acquisition cards. All single-shot readout traces are processed with various lters, down-sampling and an optimal quadrature rotation lter. The details of this full process, which is currently performed explicitly in the acquisition computer, is described in detail in parallel work27. As the devices scale up in the near term towards plaquette and logical qubit demonstrations, the multiple readout signals will need to be processed directly on Field-Programmable Gate Arrays and then correlated in post processing.

Calibration sequences. Complete tune-up of all microwave gates is accomplished using sets of automated repeated sequences. For single-qubit gates, the repeated calibration sequences are described in a previous publication28.

The cross-resonance pulse amplitude is calibrated in close analogy to single-qubit amplitude calibrations. An odd number 2N 1 of ZX90 pulses are applied

and the amplitude is adjusted so that for each N the expected signal is halfway between 0 and 1. Any amplitude miscalibrations lead to departures from this expected signal and are amplied for increasing N.

In addition to amplitude, we must also calibrate the phase of the ZX90 pulse between Q3 and Q2 (as well as Q1 and Q2). In our experiment, we use a separate microwave generator to supply the cross-resonance pulse on Q3 at the frequency of Q2. The phase of this microwave signal must be calibrated to match that of the microwave generator supplying the single-qubit pulses on Q2. This is done by applying the pulse sequence IY90(ZU180IX)NIX90. The U denotes the rotation axis

dened by the second generator and the goal is to calibrate for an X-rotation. In the case of an X-rotation, we expect the signal to be halfway between 0 and 1 for each N and miscalibration of the phase leads to deviations that are amplied with increasing N. These methods provide a routine for automated calibration with high precision. In the experiments, all cross-resonance pulses were calibrated on a regular basis because of phase drift between the two microwave generators.

Randomized benchmarking. All single-qubit gates are 40-ns Gaussian-shaped microwave pulses (Gaussian width s 10 ns) resonant with the transition fre

quencies of the qubits, with scaled derivative-of-Gaussian shapes applied on the quadrature channel to minimize leakage effects30. The gates are all autonomously calibrated with a set of repeated pulse experiments, correcting for: amplitude of X90

and X gates, amplitude imbalance between X- and Y-rotations, mixer skew and derivative-of-Gaussian-shape parameter. Single-qubit gates are all independently characterized via Clifford31 randomized benchmarking (RB) and summarized in Table 1. To characterize the addressability error of the system, we perform simultaneous23 RB, applying different sets of randomized single-qubit Clifford gates to all three qubits at the same time. These results are also summarized in Table 1 and essentially indicate that addressability errors are at the 0.1% error level.

The two-qubit ZX90 gates for both pairs of qubits are composed as a two-qubit refocusing sequence previously described32 and include shaped Gaussian turn-on (3s, s 24 ns), a at section and then a Gaussian turn-off, for a total gate time of

350 ns. The ZX90 gates are tuned-up also using repeated pulse experiments (described in previous section). It is also important to note that the pair of twoqubit gates can be applied simultaneously, as they commute with one another. To characterize the gates, we generate two-qubit Clifford operations12 and perform RB. The results for the two cases are shown in Fig. 8, where we show the average delity decay over 35 different randomized two-qubit Clifford sequences. Analysing the decay curves gives us error per two-qubit Clifford gate of0.0580.003 for Q1 and Q2, and 0.0650.002 for Q3 and Q2. We nd that the reduced chi-square for these ts are 0.583 and 0.385, respectively. This demonstrates that the model is a faithful representation of the data. As each two-qubit Clifford gate is composed of 1.5 ZX90 generators, we estimate the two-qubit ZX90 gate errors to be 3.8 and 4.3%.

Readout characterization. For this experiment, each qubit has its own measurement resonator. On Q1 and Q3 high-power readout was used and for Q2 a dispersive linear readout with a JPA was used. The readout was performed by using an integrating kernel that takes into account the response of the cavity (see ref. 27 for more details). This is important when most of the information is in the initial transients of the signal. The integration time for the experiment was 4 ms for the high-power readout and 2 ms for the dispersive readout with the JPA.

Shown in Fig. 9 are typical histograms for the three readout channels averaged over all computational basis for the qubits not measured. Here we see that the assignment delity, dened by

1 P 0j1

=2 P 1j0

=2 1

6 NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5015 ARTICLE

Master Slave

Slave

Q I

Slave

TEK5014

Ch1 Ch2 Ch3 Ch4

BBNAPS

Ch1 Ch2 Ch3 Ch4

BBNAPS

Ch1 Ch2 Ch3 Ch4

TEK5014

Ch1 Ch2 Ch3 Ch4

TEK7052

Ch1 Ch2

MXG

N5128A

MXG

N5128A

MXG

N5128A

E8257D

1 2

Q I

Q I Q I

1 2

IF2

IF1

IF3

Dilution refrigerator

YOKO GS200

300 K

Ch1 Ch2

ALAZARTECH ATS9870

Ch1 Ch2

ALAZARTECH ATS9870

Qubit 1

Qubit 2

Pump

Qubit 3

HEMT

5GHz

1GHz

JPA

Q1 Q2 Q3

10 dB

70 K

20 dB

4 K

Optimal quadrature rotation filter

800 mK

30 dB

15 mK

Eccosorb filter

Figure 7 | Detailed schematic of experimental setup. Wiring scheme for all room temperature controls as well as internal conguration of Oxford Instruments Triton dilution refrigerator.

Table 1 | Summary of single-qubit randomized benchmarking.

Randomize Error ( 10 3)

Q1 3.060.05

Q2 2.300.05

Q3 2.770.05 Q1, Q2, Q3 3.80.1 3.320.09 2.890.05

0.4 0 5 10 15 20 25

0.5

ZX90 for Q1 and Q2

ZX90 for Q3 and Q2

0.9

0.8

P 0

0.7

0.6

for the three channels is 0.84, 0.91 and 0.89, respectively. These are typical values and we see about a 23% uctuation over the course of a typical experiment. By tting a double-Gaussian model to the data, we nd that the ratio of the undesired state to the desired state for Q1 prepared in the ground (excited) is 9.9% (22%), for Q2 5.7% (8.3%) and for Q3 13.7% (6.0%). We believe most of the error is due to the high-power nonlinear readout of Q1 and Q3 and is not due to the qubits being initialized in the wrong state. With no power applied to the Q1 and Q3 resonator, the assignment delity is 0.95 and the ratios of the two Gaussians are 5.6% when Q2 is prepared in the excited state and negligible when Q2 is prepared in the ground state.

State tomography. For state tomography we used the correlation method as described in ref. 27. The single shots for each measurement resonator are correlated and from a set of complete post rotations we can use either linear inversion or an SDP (with constraints rr and tr(r) 1) to reconstruct the state. The complete set

of rotations used are {I, X, X90, X 90, Y90, Y 90}#n.

Typically, 20,000 shots for each post rotation are used and we nd that the statistical error in the measured voltages has signal-to-noise ratio (SNR) B1 104,

2 104, 2 104 for the three measurement channels M1, M2, M3, respectively. The

second-order correlators range in SNR from B2 to 5 103 and the third order has

SNR B1,000. Using these and a bootstrapping method28 we can estimate the state

delity and the statistical error. The state delity is given by

Fstate

Number of cliffords

Figure 8 | Two-qubit randomized benchmarking. Average P0, population of Q2 ground state, versus number of two-qubit Cliffords generated via ZX90 gates between Q1 (Q3) and Q2 is shown as red (blue) circles.

Experiments are performed randomizing over 35 different sequences of Cliffords. Fits to the RB experiment for Q1 (Q3) and Q2 are shown as solid red (blue) lines, from which we extract an error per two-qubit Clifford of0.0580.003 (0.0650.002).

h i

; 2

rideal

p rnoisy

rideal

where rideal is the ideal state and rnoisy is the reconstructed state.

We nd that in all cases the uctuations in the state delity from statistics is much smaller than the difference between the linear reconstruction and the SDP. Furthermore, we nd typically the sum of all the negative eigenvalues in the three-qubit space to be o0.03.

NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5015

0.06

0.4 0.3 0.2 0.1 0 0.1 0.2 0.3 0.4

0.04

0.02

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45

0.06

Count fraction

0.04

0.02

0.06

0.04

0.02

0.2 0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

M2 voltage

Figure 9 | Histogram of three independent readouts. (a) Histograms of the high-power readout used for Q1. (b) Histograms of the dispersive linear readout with a JPA for Q2. (c) Histograms of the high-power readout used for Q3. In all cases blue is preparing ground and red is excited. Solid lines are double-Gaussian ts to the histograms.

Measurement tomography. An ideal Z-parity check can be described by the quantum operation

evenr even 1

j i 1

h j oddr odd 0

j i 0

h j; 3

where

even II ZZ

=2 odd II ZZ

=2 4 and the extra system is used to label the outcome of measurement of the syndrome qubit. In the noisy case this is represented by the operation

E r

Eeven r

j i 1

h j Eodd r

j i 0

h j 5 and the goal of measurement tomography is to determine the conditional maps

Eeven(r) and Eodd(r). These quantum operations are completely positive but not

trace preserving.

By binning the results of the measurement on the syndrome qubit, tomography on the two-qubit subspace is performed by preparing a complete set of different input states and measurement bases via pre and post rotations, and reconstructing the operations from the measurement results. The complete set of rotations that we use are the same as those used in state tomography. We use both a linear reconstruction and a minimization to make the maps physical. For more details on how quantum process tomography can be performed, see ref. 28.

We use the Pauli transfer matrix28 dened by

RL i; j

Tr PiE Pj

d 6 to represent the measurement operations where Pjs are the standard Pauli operators {I, X, Y, Z}#2.

To quantify the measurement, we dene the measurement delity by a generalization of the average delity. Since the measurement maps are not trace preserving, we need to use normalized outputs 0xc, and E0xc, where

A0x c

Ax c

j i c

h j

odd r : 9

Since this is a quantum operation, the standard delity between quantum operations can be used.

References

1. Aharonov, D. & Ben-Or, M. Fault-tolerant quantum computation with constant error. In STOC 97, 176188 (ACM, New York, NY, 1997).

2. Preskill, J. Reliable quantum computers. Proc. R. Soc. Lond. A 454, 385410 (1998).

3. Knill, E. Quantum computing with realistically noisy devices. Nature 434, 3944 (2005).

4. Bravyi, S. & Kitaev, A. Quantum codes on a lattice with boundary. Preprint at http://arXiv.org/quant-ph/9811052

Web End =http://arXiv.org/quant-ph/9811052 (1998).

5. Dennis, E., Kitaev, A., Landahl, A. & Preskill, J. Topological quantum memory.J. Math. Phys. 43, 44524505 (2002).6. Raussendorf, R. & Harrington, J. Fault-tolerant quantum computation with high threshold in two dimensions. Phys. Rev. Lett. 98, 190504 (2007).

7. Fowler, A. G., Mariantoni, M., Martinis, J. M. & Cleland, A. N. Surface codes: towards practical large-scale quantum computation. Phys. Rev. A. 86, 032324 (2012).

8. Paik, H. et al. Observation of high coherence in Josephson junction qubits measured in a three-dimensional circuit QED architecture. Phys. Rev. Lett. 107, 240501 (2011).

9. Rigetti, C. et al. Superconducting qubit in a waveguide cavity with a coherence time approaching 0.1ms. Phys. Rev. B 86, 100506 (2012).

10. Chang, J. B. et al. Improved superconducting qubit coherence using titanium nitride. Appl. Phys. Lett. 103, 012602 (2013).

11. Barends, R. et al. Coherent Josephson qubit suitable for scalable quantum integrated circuits. Phys. Rev. Lett. 111, 080502 (2013).

12. Corcoles, A. D. et al. Process verication of two-qubit quantum gates by randomized benchmarking. Phys. Rev. A 87, 030301 (2013).

13. Barends, R. et al. Superconducting quantum circuits at the surface code threshold for fault tolerance. Nature 508, 500503 (2014).

14. Bergeal, N. et al. Phase-preserving amplication near the quantum limit with a Josephson ring modulator. Nature 465, 6468 (2010).

15. Riste, D., van Leeuwen, J. G., Ku, H.-S., Lehnert, K. W. & DiCarlo, L. Initialization by measurement of a superconducting quantum bit circuit. Phys. Rev. Lett. 109, 050507 (2012).

16. Johnson, J. E. et al. Heralded state preparation in a superconducting qubit. Phys. Rev. Lett. 109, 050506 (2012).

17. Majer, J. et al. Coupling superconducting qubits via a cavity bus. Nature 449, 443447 (2007).

18. Johnson, B. R. et al. Quantum non-demolition detection of single microwave photons in a circuit. Nat. Phys. 6, 663667 (2010).

Ax cj i c

h j

k ktr, tr refers to the trace norm and x {even,

odd}. Doing this gives

Fxmeas Z

dcTr

0x c

q E0x

0x c

: 7

Since the nullspace of a projection operation has measure zero and the noisy realization typically will also have a nullspace of zero measure, this integral is well dened. To compute this, we draw 150,000 different random states from the FubiniStudy measure and compute the average.

One could also dene a process delity by computing the state delity between normalized Choi matrices of the ideal and noisy operations

Fxpro

h i

Tr rxideal

rxnoisy

h i

rxideal

rxnoisy

rxideal

; 8

however, for non-unitary processes there is no simple relationship between them.

The unconditional map can be dened by tracing equation (5) over the syndrome qubit giving

E r

Eeven r

8 NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms5015 ARTICLE

19. DiVincenzo, D. P. Fault-tolerant architectures for superconducting qubits. Phys. Scr. T137, 014020 (2009).

20. DiCarlo, L. et al. Preparation and measurement of three-qubit entanglement in a superconducting circuit. Nature 467, 574578 (2010).

21. Groen, J. P. et al. Partial-measurement backaction and nonclassical weak values in a superconducting circuit. Phys. Rev. Lett. 111, 090506 (2013).

22. Hatridge, M., Vijay, R., Slichter, D. H., Clarke, J. & Siddiqi, I. Dispersive magnetometry with a quantum limited SQUID parametric amplier. Phys. Rev. B 83, 134501 (2011).

23. Gambetta, J. M. et al. Characterization of addressability by simultaneous randomized benchmarking. Phys. Rev. Lett. 109, 240504 (2012).

24. Paraoanu, G. S. Microwave-induced coupling of superconducting qubits. Phys. Rev. B 74, 140504 (2006).

25. Rigetti, C. & Devoret, M. Fully microwave-tunable universal gates in superconducting qubits with linear couplings and xed transition frequencies. Phys. Rev. B 81, 134507 (2010).

26. Reed, M. D. et al. High-delity readout in circuit quantum electrodynamics using the Jaynes-Cummings nonlinearity. Phys. Rev. Lett. 105, 173601 (2010).

27. Ryan, C. A. et al. Tomography via Correlation of Noisy Measurement Records. Preprint at http://arxiv.org/abs/1310.6448

Web End =http://arxiv.org/abs/1310.6448 (2013).

28. Chow, J. M. et al. Universal quantum gate set approaching fault-tolerant thresholds with superconducting qubits. Phys. Rev. Lett. 109, 060501 (2012).

29. Merkel, S. T. et al. Self-consistent quantum process tomography. Phys. Rev. A 87, 062119 (2013).

30. Motzoi, F., Gambetta, J. M., Rebentrost, P. & Wilhelm, F. K. Simple pulses for elimination of leakage in weakly nonlinear qubits. Phys. Rev. Lett. 103, 110501 (2009).

31. Magesan, E., Gambetta, J. M. & Emerson, J. Scalable and robust randomized benchmarking of quantum processes. Phys. Rev. Lett. 106, 180504 (2011).32. Corcoles, A. D. et al. Protecting superconducting qubits from radiation. Appl. Phys. Lett. 99, 181906 (2011).

33. Saira, O.-P. et al. Entanglement genesis by ancilla-based parity measurement in 2D circuit QED. Phys. Rev. Lett. 112, 070502 (2014).

Acknowledgements

We thank M.B. Rothwell and G.A. Keefe for fabricating devices. We thank J.R. Rozen andJ. Rohrs for experimental contributions, and D. DiVincenzo and G. Smith for engaging theoretical discussions. We thank I. Siddiqi for providing the Josephson Parametric Amplier. We acknowledge Caltech for HEMT ampliers. We acknowledge support from IARPA under contract W911NF-10-1-0324. All statements of fact, opinion or conclusions contained herein are those of the authors and should not be construed as representing the ofcial views or policies of the US Government.

Author contributions

J.M.C. and J.M.G. designed the experiments. J.M.G. and J.A.S. came up with the proposed lattice. D.W.A. and J.M.G. performed simulations of devices. J.M.C., D.W.A. and S.J.S. characterized devices. J.M.C., J.M.G., A.W.C. and E.M. interpreted and analysed the experimental data. N.A.M. set-up the microwave control hardware. E.M. and J.M.G. developed the measurement tomography formalism. B.R.J. and C.A.R. developed readout hardware and analysis for the correlated readouts. All authors contributed to the composition of the manuscript.

Additional information

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Chow, J. M. et al. Implementing a strand of a scalable fault-tolerant quantum computing fabric. Nat. Commun. 5:4015doi: 10.1038/ncomms5015 (2014).

NATURE COMMUNICATIONS | 5:4015 | DOI: 10.1038/ncomms5015 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

Word count: 7853

Show less

Abstract

Translate

With favourable error thresholds and requiring only nearest-neighbour interactions on a lattice, the surface code is an error-correcting code that has garnered considerable attention. At the heart of this code is the ability to perform a low-weight parity measurement of local code qubits. Here we demonstrate high-fidelity parity detection of two code qubits via measurement of a third syndrome qubit. With high-fidelity gates, we generate entanglement distributed across three superconducting qubits in a lattice where each code qubit is coupled to two bus resonators. Via high-fidelity measurement of the syndrome qubit, we deterministically entangle the code qubits in either an even or odd parity Bell state, conditioned on the syndrome qubit state. Finally, to fully characterize this parity readout, we develop a measurement tomography protocol. The lattice presented naturally extends to larger networks of qubits, outlining a path towards fault-tolerant quantum computing.

Details

Title

Implementing a strand of a scalable fault-tolerant quantum computing fabric

Author

Chow, Jerry M; Gambetta, Jay M; Magesan, Easwar; Abraham, David W; Cross, Andrew W; Johnson, B R; Masluk, Nicholas A; Ryan, Colm A; Smolin, John A; Srinivasan, Srikanth J; Steffen, M

Pages

4015

Publication year

2014

Publication date

Jun 2014

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms5015

ProQuest document ID

1539438908

Implementing a strand of a scalable fault-tolerant quantum computing fabric

Jump to:

Full text

Abstract

Details

Suggested sources