Abstract
Now that it is possible to achieve measurement and control fidelities for individual quantum bits (qubits) above the threshold for fault tolerance, attention is moving towards the difficult task of scaling up the number of physical qubits to the large numbers that are needed for fault-tolerant quantum computing1,2. In this context, quantum-dot-based spin qubits could have substantial advantages over other types of qubit owing to their potential for all-electrical operation and ability to be integrated at high density onto an industrial platform3,4,5. Initialization, readout and single- and two-qubit gates have been demonstrated in various quantum-dot-based qubit representations6,7,8,9. However, as seen with small-scale demonstrations of quantum computers using other types of qubit10,11,12,13, combining these elements leads to challenges related to qubit crosstalk, state leakage, calibration and control hardware. Here we overcome these challenges by using carefully designed control techniques to demonstrate a programmable two-qubit quantum processor in a silicon device that can perform the DeutschâJosza algorithm and the Grover search algorithmâcanonical examples of quantum algorithms that outperform their classical analogues. We characterize the entanglement in our processor by using quantum-state tomography of Bell states, measuring state fidelities of 85â89 per cent and concurrences of 73â82 per cent. These results pave the way for larger-scale quantum computers that use spins confined to quantum dots.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 /Â 30Â days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
$199.00 per year
only $3.90 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
References
Barends, R. et al. Superconducting quantum circuits at the surface code threshold for fault tolerance. Nature 508, 500â503 (2014)
Debnath, S. et al. Demonstration of a small programmable quantum computer with atomic qubits. Nature 536, 63â66 (2016)
Loss, D. & DiVincenzo, D. P. Quantum computation with quantum dots. Phys. Rev. A 57, 120â126 (1998)
Maurand, R. et al. A CMOS silicon spin qubit. Nat. Commun. 7, 13575 (2016)
Vandersypen, L. M. K. et al. Interfacing spin qubits in quantum dots and donorsâhot, dense, and coherent. npj Quantum Inf. 3, 34 (2017)
Shulman, M. D. et al. Demonstration of entanglement of electrostatically coupled singlet-triplet qubits. Science 336, 202â205 (2012)
Kim, D. et al. Quantum control and process tomography of a semiconductor quantum dot hybrid qubit. Nature 511, 70â74 (2014)
Veldhorst, M. et al. A two-qubit logic gate in silicon. Nature 526, 410â414 (2015)
Medford, J. et al. Self-consistent measurement and state tomography of an exchange-only spin qubit. Nat. Nanotechnol. 8, 654â659 (2013)
Vandersypen, L. M. K. et al. Experimental realization of Shorâs quantum factoring algorithm using nuclear magnetic resonance. Nature 414, 883â887 (2001)
DiCarlo, L. et al. Demonstration of two-qubit algorithms with a superconducting quantum processor. Nature 460, 240â244 (2009)
Gulde, S. et al. Implementation of the DeutschâJozsa algorithm on an ion-trap quantum computer. Nature 421, 48â50 (2003)
van der Sar, T. et al. Decoherence-protected quantum gates for a hybrid solid-state spin register. Nature 484, 82â86 (2012)
Bluhm, H. et al. Dephasing time of GaAs electron-spin qubits coupled to a nuclear bath exceeding 200 μs. Nat. Phys. 7, 109â113 (2011)
Petta, J. R. et al. Coherent manipulation of coupled electron spins in semiconductor quantum dots. Science 309, 2180â2184 (2005)
Zwanenburg, F. A. et al. Silicon quantum electronics. Rev. Mod. Phys. 85, 961â1019 (2013)
Tyryshkin, A. M. et al. Electron spin coherence exceeding seconds in high-purity silicon. Nat. Mater. 11, 143â147 (2012)
Veldhorst, M. et al. An addressable quantum dot qubit with fault-tolerant control-fidelity. Nat. Nanotechnol. 9, 981â985 (2014)
Muhonen, J. T. et al. Storing quantum information for 30 seconds in a nanoelectronic device. Nat. Nanotechnol. 9, 986â991 (2014)
Kawakami, E. et al. Gate fidelity and coherence of an electron spin in an Si/SiGe quantum dot with micromagnet. Proc. Natl Acad. Sci. USA 113, 11738â11743 (2016)
Yoneda, J. et al. A quantum-dot spin qubit with coherence limited by charge noise and fidelity higher than 99.9%. Nat. Nanotechnol. 13, 102â106 (2018)
Kawakami, E. et al. Electrical control of a long-lived spin qubit in a Si/SiGe quantum dot. Nat. Nanotechnol. 9, 666â670 (2014)
Elzerman, J. M. et al. Single-shot read-out of an individual electron spin in a quantum dot. Nature 430, 431â435 (2004)
Srinivasa, V., Nowack, K. C., Shafiei, M., Vandersypen, L. M. K. & Taylor, J. M. Simultaneous spin-charge relaxation in double quantum dots. Phys. Rev. Lett. 110, 196803 (2013)
Pioro-Ladrière, M. et al. Electrically driven single-electron spin resonance in a slanting Zeeman field. Nat. Phys. 4, 776â779 (2008)
Knill, E. et al. Randomized benchmarking of quantum gates. Phys. Rev. A 77, 012307 (2008)
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)
Meunier, T., Calado, V. E. & Vandersypen, L. M. K. Efficient controlled-phase gate for single-spin qubits in quantum dots. Phys. Rev. B 83, 121403 (2011)
Zajac, D. M. et al. Resonantly driven CNOT gate for electron spins. Science 359, 439â442 (2018)
Deutsch, D. & Jozsa, R. Rapid solution of problems by quantum computation. Proc. R. Soc. Lond. A 439, 553â558 (1992)
Grover, L. K. Quantum mechanics helps in searching for a needle in a haystack. Phys. Rev. Lett. 79, 325â328 (1997)
Reed, M. D. et al. Reduced sensitivity to charge noise in semiconductor spin qubits via symmetric operation. Phys. Rev. Lett. 116, 110402 (2016)
Martins, F. et al. Noise suppression using symmetric exchange gates in spin qubits. Phys. Rev. Lett. 116, 116801 (2016)
Yang, C. H. et al. Spin-valley lifetimes in a silicon quantum dot with tunable valley splitting. Nat. Commun. 4, 2069 (2013)
James, D. F. V., Kwiat, P. G., Munro, W. J. & White, A. G. Measurement of qubits. Phys. Rev. A 64, 052312 (2001)
Das Sarma, S., Wang, X. & Yang, S. Hubbard model description of silicon spin qubits: charge stability diagram and tunnel coupling in Si double quantum dots. Phys. Rev. B 83, 235314 (2011)
Dial, O. E. et al. Charge noise spectroscopy using coherent exchange oscillations in a singlet-triplet qubit. Phys. Rev. Lett. 110, 146804 (2013)
Ithier, G. et al. Decoherence in a superconducting quantum bit circuit. Phys. Rev. B 72, 134519 (2005)
Acknowledgements
This research was sponsored by the Army Research Office (ARO) under grant numbers W911NF-17-1-0274 and W911NF-12-1-0607. The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the official policies, either expressed or implied, of the ARO or the US Government. The US Government is authorized to reproduce and distribute reprints for government purposes notwithstanding any copyright notation herein. Development and maintenance of the growth facilities used for fabricating samples is supported by DOE (DE-FG02-03ER46028). We acknowledge the use of facilities supported by NSF through the University of Wisconsin-Madison MRSEC (DMR-1121288). E.K. was supported by a fellowship from the Nakajima Foundation. We acknowledge financial support from the Marie SkÅodowska-Curie actionsâNanoscale solid-state spin systems in emerging quantum technologiesâSpin-NANO, grant agreement number 676108. We acknowledge discussions with S. Dobrovitski, C. Dickel, A. Rol, J. P. Dehollain, Z. Ramlakhan and members of the Vandersypen group, and technical assistance from R. Schouten, R. Vermeulen, M. Tiggelman, M. Ammerlaan, J. Haanstra, R. Roeleveld and O. Benningshof.
Author information
Authors and Affiliations
Contributions
T.F.W. performed the experiment with help from E.K. and P.S., T.F.W. and S.G.J.P. analysed the data, S.G.J.P. performed the simulations of the algorithms, T.F.W., S.G.J.P., E.K., P.S., M.V., M.F., S.N.C., M.A.E. and L.M.K.V. contributed to the interpretation of the data and commented on the manuscript, D.R.W. fabricated the device, D.E.S. and M.G.L. grew the Si/SiGe heterostructure, T.F.W. wrote the manuscript (S.G.J.P. wrote parts of Methods) and L.M.K.V. conceived and supervised the project.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Reviewer Information Nature thanks H. Bluhm and the other anonymous reviewer(s) for their contribution to the peer review of this work.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Figure 1 Schematic of the measurement set-up.
The sample was bonded to a printed circuit board (PCB) mounted onto the mixing chamber of a dilution refrigerator. All measurements were performed at the base temperature of the fridge, Tbaseâââ20âmK. DC voltages were applied to all of the gate electrodes using room-temperature digital-to-analogue converters via filtered lines (not shown). Voltage pulses were applied to plunger gates P1 and P2 using a Tektronix 5014C arbitrary waveform generator (AWG) with 1-GHz clock rate. The signals from the AWGs passed through a room-temperature low-pass filter and attenuators at different stages of the fridge and were added to the DC signals via bias tees mounted on the PCB. Two Keysight E8267D vector microwave sources, MW1 and MW2, were used to apply microwaves (18â20âGHz) to perform EDSR on Q1 and Q2, respectively. The signals passed through room-temperature DC blocks and custom-built 15-GHz high-pass filters and attenuators at different stages of the fridge and were added to the DC signals via bias tees mounted on the PCB. The output of the microwave source (phase, frequency, amplitude and duration) was controlled with I/Q vector modulation. The I/Q signals were generated with another Tektronix 5041C, which was the master device for the entire set-up and provided trigger signals for the other devices. In addition to the vector modulation, we used pulse modulation to give an on/off microwave power output ratio of 120âdB. Although I/Q modulation can be used to output multiple frequencies, the bandwidth of the AWG was not sufficient to control both qubits with one microwave source owing to their large separation in frequency (1.3âGHz). The sensor current I was converted to a voltage signal with a custom-built preamplifier, and an isolation amplifier was used to separate the signal ground, with the measurement equipment grounded to reduce interference. Following this, a 20-kHz Bessel low-pass filter was applied to the signal using a SIM965 analogue filter. A field-programmable gate array (FPGA) analysed the voltage signal during the readout and assigned the trace to be spin-up if the voltage fell below a certain threshold. The voltage signal could also be measured with a digitizer card in the computer. The shapes of the pulses generated by the AWGs and microwave sources during qubit manipulation, along with the typical timescales, are shown in the lower right. Square pulses were used to perform the CZ gate and as the input for the I/Q modulation to generate microwave pulses. The pulse modulation was turned on 40âns before turning on the I/Q signal, owing to the time needed for the modulation to switch on.
Extended Data Figure 2 Measurement protocol for two electron spins.
a, Stability diagram of the double quantum dot, showing the positions in gate space used to perform single-qubit gates (red circle) and two-qubit gates (yellow circle). The map shows the numerically obtained differential current dI/dVP1 through the quantum-dot charge sensor as a function of the voltages on P1 and P2. The white dashed line is the (1, 1)â(0,2) inter-dot transition line. The white arrow indicates the detuning axis ε used in the experiments. Although the detuning pulse for the two-qubit gate crosses the charge-addition lines of D1 and D2, the quantum dots remain in the (1, 1) charge state because the pulse time is much shorter than the electron tunnel times to the reservoirs. b, Plot of the voltage pulses applied to plunger gates P1 and P2 and the response of the quantum-dot charge sensor over one measurement cycle. D2 is unloaded by pulsing into the (1, 0) charge region for 1.5âms (purple circle). The electron on D1 is initialized to spin-down by pulsing to a spin-relaxation hotspot at the (1, 0) and (0, 1) charge degeneracies (orange circle) for 50âμs (see Extended Data Fig. 5). D2 is loaded with a spin-down electron by pulsing to the readout position for 4âms (blue circle). During manipulation, the voltages on the plunger gates are pulsed to the red circle for single-qubit gates and to the yellow circle for two-qubit gates, where the exchange is about 6âMHz. After manipulation, the spin of the electron on D2 is measured by pulsing to the readout position (blue circle) for 0.7âms, where the Fermi level of the reservoir is between the spin-up and spin-down electrochemical potentials of D2. If the electron is spin-up then it can tunnel out, after which a spin-down electron tunnels back in. These two tunnel events are detected by the quantum-dot sensor as a single blip in the current signal. An additional 1.3âms is spent at the readout position so that D2 is initialized to spin-down with high fidelity. Following this, Q1 is measured by first performing a CROT gate at the yellow circle so that , where CROT12 indicates a CROT gate with Q1 as the control and Q2 as the target. A projective measurement of Q1 is then performed by measuring Q2 at the readout position for 0.7âms (blue circle). Finally, we add a compensation pulse to VP1 and VP2 so that over the measurement cycle VDCâ=â0 to mitigate charging effects in the bias tees. c, Close-up of the stability diagram in a showing the positions in gate space used for initialization and readout.
Extended Data Figure 3 Single-qubit properties and two-axis control.
The purple (top) and orange (bottom) data points correspond to measurements performed on Q1 and Q2, respectively, in the (1, 1) regime (red circle in Extended Data Fig. 2). a, Spin-up fraction as a function of the microwave frequency of an applied Ï pulse, showing a resonant frequency of 18.387âGHz (19.670âGHz) for Q1 (Q2). b, The spin-relaxation time is measured by preparing the qubit to spin-up and varying the wait time before readout. From the exponential decay in the spin-up probability we measure T1â>â50âms (T1â=â3.7â±â0.5âms) for Q1 (Q2). c, Spin-up probability as a function of microwave duration, showing Rabi oscillations of 2.0âMHz for Q1 and Q2. d, The dephasing time is measured by applying a Ramsey pulse sequence and varying the free evolution time Ï. Oscillations were added artificially to improve the fit of the decay by making the phase of the last microwave pulse dependent on the free evolution time, Ïâ=âsin(ÏÏ), where Ïâ=â4âMHz. By fitting the data with a Gaussian decay, , we extract () for Q1 (Q2). In the measurement for Q1, the first Ï/2 microwave pulse is a Y gate. The Ramsey measurement was performed over about 20âmin with the frequency calibrated every approximately 1âmin. e, The coherence time of Q1 (Q2) can be extended to T2Hahnâ=â19â±â3âμs (T2Hahnâ=â7â±â1âμs) by a Hahn echo sequence. The coherence time is extracted from an exponential fit to the spin-up probability as a function of the free evolution time in the Hahn echo sequence. f, Full two-axis control is demonstrated by applying two Ï/2 pulses and varying the phase of the second one. All error bars are 1Ï from the mean.
Extended Data Figure 4 Randomized benchmarking of single-qubit gates.
Randomized benchmarking of the single-qubit gates for each qubit is performed by applying a randomized sequence of a varying number of Clifford gates m to either the |1ã or |0ã state and measuring the final spin-up probability or P|1ã, respectively. All gates in the Clifford group are decomposed into gates from the set {I, ±X, ±X2, ±Y, ±Y2}. The purple (orange) data points show the difference in the spin-up probabilities for Q1 (Q2) as a function of sequence length. For each sequence length m we average over 32 different randomized sequences. From an exponential fit (solid lines) of the data, , we estimate average Clifford-gate fidelities FCâ=â1 â (1 â p)/2 of 98.8% and 98.0% for Q1 and Q2, respectively. The last three data points from both datasets were omitted from the fits because they deviate from a single exponential20. All error bars are 1Ï from the mean.
Extended Data Figure 5 Spin-relaxation hotspots used for high-fidelity initialization.
a, Close-up stability diagram of the (1, 0)â(0, 1) charge transition. The white arrow defines the detuning axis between D1 and D2 controlled with P1. b, Schematic of the energy-level diagram as a function of detuning for one electron spin in a double quantum dot. c, Spin-relaxation hotspots are measured by preparing the electron on D1 to spin-up using EDSR, applying a voltage pulse along the detuning axis (white arrow in a) for a wait time of 200âns and performing readout of the electron spin. We observe three dips in the spin-up probability, which correspond to spin-relaxation hot spots. The first and third hotspot are due to anticrossings between the (0, â) and (â, 0) states and the (â, 0) and (0, â) states24. The second hotspot occurs at zero detuning. The voltage separation between the first and third hot spot corresponds to the sum of the Zeeman energy of D1 and D2 divided by the gate lever arm α along the detuning axis. Knowing precisely the Zeeman energies from EDSR spectroscopy, we can accurately extract the gate lever arm to be αâ=â0.09e. d, The spin-relaxation time at zero detuning (orange circle in a) is found to be T1â=â220âns by measuring the exponential decay of the spin-up probability as a function of wait time Ï at zero detuning. All error bars are 1Ï from the mean.
Extended Data Figure 6 Two-qubit CROT gate.
a, Microwave spectroscopy of Q2 close to zero detuning between the (1, 1) and (0, 2) states (yellow dot in Extended Data Fig. 2a) with the exchange coupling on. The blue and red curves show the resonance of Q2 after preparing Q1 into spin-down and spin-up, respectively. The resonance frequency of Q2 shifts by the exchange coupling, and by applying a Ï pulse at one of these frequencies we can perform a CROT gate, which is equivalent to a CNOT gate up to a rotation. As discussed in the main text, this CROT gate is used to perform the projective measurement of Q1. All error bars are 1Ï from the mean.
Extended Data Figure 7 Measurement of Joff using a decoupling sequence.
The exchange coupling Joff during single-qubit gates is measured using a two-qubit Hahn echo sequence, which cancels out any unconditional rotations during the free evolution time Ï. Fitting the spin-up probability as a function of free evolution time Ï using the functional form sin(2ÏJoffÏ), we extract Joffâ=â270âkHz. All error bars are 1Ï from the mean.
Extended Data Figure 8 Microwave spectroscopy of Q1 and Q2.
a, b, Spectroscopy of Q1 (a) and Q2 (b) versus detuning energy ε after initializing the other qubit to . Towards εâ=â0 there are two resonances each for Q1 and Q2, which are separated by the exchange energy J(ε)/h. As discussed, the Zeeman energy EZ(ε) of Q1 and Q2 also depends on detuning because changes to the applied voltages shift the position of the electron in the magnetic-field gradient. The four resonance frequencies are fitted (green, blue, red and yellow lines) with fjkâ=âEZj(ε)â+â(â1)k+1J(ε), where j denotes the qubit and k denotes the state of the other qubit. The data are fitted well using , , and EZ2(ε) â ε. The fitted Zeeman energies of Q1 and Q2 are shown by the black lines. We observe that the Zeeman energy of Q1 has an exponential dependence towards the (0, 2) charge regime (εâ=â0). This observation can be explained by the electron delocalizing from D1 towards D2, which has a much higher Zeeman energy. c, Schematic of the colour-coded transitions that correspond to the resonances in a and b.
Extended Data Figure 9 Decay of the decoupled CZ oscillations.
The normalized spin-up probability of Q1 as a function of the total duration time 2Ï of the two CZ gates in the DCZ sequence. The data are fitted using a sinusoid, P|1ãâ=â0.5sin(2ÏJÏ)â+â0.5, with either a Gaussian (black line; ) or exponential (red line; ) decay. From these fits we find a decay time of T2â=â1.6âμs. All error bars are 1Ï from the mean.
Extended Data Figure 10 Simulation of the DeutschâJosza and Grover search algorithms using the DCZ gate.
a, b, Two-spin probabilities as a function of the sequence time during the DeutschâJosza algorithm (a) and the Grover search algorithm (b) for each function, using the decoupled version of the two-qubit CZ gate (the DCZ gate). The solid lines show the outcome of the simulations, which include decoherence due to quasi-static charge noise and nuclear-spin noise. All error bars are 1Ï from the mean.
Supplementary information
Supplementary Information
This file contains supplementary notes S1-S3 and supplementary figures S1-S5. (PDF 1317 kb)
Rights and permissions
About this article
Cite this article
Watson, T., Philips, S., Kawakami, E. et al. A programmable two-qubit quantum processor in silicon. Nature 555, 633â637 (2018). https://doi.org/10.1038/nature25766
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nature25766
This article is cited by
-
The SpinBus architecture for scaling spin qubits with electron shuttling
Nature Communications (2024)
-
Spin-EPR-pair separation by conveyor-mode single electron shuttling in Si/SiGe
Nature Communications (2024)
-
Anisotropic exchange interaction of two hole-spin qubits
Nature Physics (2024)
-
Data needs and challenges for quantum dot devices automation
npj Quantum Information (2024)
-
A long lifetime floating on neon
Nature Physics (2024)