∗
A Quantum Engineer’s Guide to Superconducting Qubits P. Krantz1,2,†, M. Kjaergaard1, F. Yan1, T.P. Orlando1, S. Gustavsson1, and W. D. Oliver1,3,‡ 1Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 2Wallenberg Centre for Quantum Technology (WACQT), Chalmers University of Technology, Gothenburg, SE-41296, Sweden and 3MIT Lincoln Laboratory, 244 Wood Street, Lexington, MA 02420, USA
(Dated: 9 July 2021) The aim of this review is to provide quantum engineers with an introductory guide to the central concepts and challenges in the rapidly accelerating field of superconducting quantum circuits. Over the past twenty years, the field has matured from a predominantly basic research endeavor to one that increasingly explores the engineering of larger-scale superconducting quantum systems. Here, we review several foundational elements – qubit design, noise properties, qubit control, and readout techniques – developed during this period, bridging fundamental concepts in circuit quantum electrodynamics (cQED) and contemporary, state- of-the-art applications in gate-model quantum computation.
CONTENTS 1. Connecting T1 to S(ω) 19 2. Connecting Tϕ to S(ω) 20 I. Introduction 2 3. Noise spectroscopy 21 A. Organization of this article 2 E. Engineering noise mitigation 21 1. Materials and fabrication improvements 22 II. Engineering quantum circuits 3 2. Design improvements 22 A. From quantum harmonic oscillator to the 3. Dynamical error suppression 22 transmon qubit 3 4. Cryogenic engineering 22 B. Qubit Hamiltonian engineering 6 1. Tunable qubit: split transmon 6 IV. Qubit control 22 2. Towards larger anharmonicity: flux qubit A. Boolean logic gates used in classical and fluxonium 7 computers 23 C. Interaction Hamiltonian engineering 9 B. Quantum logic gates used in quantum 1. Physical coupling: capacitive and computers 24 inductive 9 C. Comparing classical and quantum gates 26 2. Coupling axis: transverse and 1. Gate sets and gate synthesis 27 longitudinal 10 2. Addressing superconducting qubits 27 D. Single-qubit gates 27 III. Noise, decoherence, and error mitigation 12 1. Capacitive coupling for X,Y control 28 A. Types of noise 12 2. Virtual Z gate 30 1. Systematic noise 12 3. The DRAG scheme 30 2. Stochastic noise 12 E. The iSWAP two-qubit gate in tunable 3. Noise strength and qubit susceptibility 12 qubits 31 B. Modeling noise and decoherence 12 1. Deriving the iSWAP unitary 31 1. Bloch sphere representation 12 2. Applications of the iSWAP gate 33 arXiv:1904.06560v5 [quant-ph] 7 Jul 2021 2. Bloch-Redfield model of decoherence 13 F. The CPHASE two-qubit gate in tunable 3. Modification due to 1/f-type noise 16 qubits 33 4. Noise power spectral density (PSD) 17 1. Trajectory design for the CPHASE gate 35 C. Common examples of noise 18 2. The CPHASE gate for quantum error 1. Charge noise 18 correction 35 2. Magnetic flux noise 18 3. Quantum simulation and algorithm 3. Photon number fluctuations 19 demonstrations using CPHASE 36 4. Quasiparticles 19 G. Two-qubit gates using only microwaves 37 D. Operator form of qubit-environment 1. The operational principle of the CR gate 37 interaction 19 2. Improvements to the CR gate and quantum error correction experiments using CR 38 3. Quantum simulation and algorithm † ‡ ∗ [email protected], [email protected] demonstrations with the CR gate 38 2
4. Other microwave-only gates: bSWAP, and polarized photons30–33, where the quantum informa- MAP, and RIP 39 tion is encoded in natural microscopic quantum systems, H. Gate implementations with tunable superconducting qubits are macroscopic in size and litho- coupling 40 graphically defined. One remarkable feature of superconducting qubits is V. Qubit readout 41 that their energy-level spectra are governed by circuit el- A. Dispersive readout 41 ement parameters and thus are configurable; they can be B. Measuring the resonator amplitude and designed to exhibit “atom-like” energy spectra with de- phase 43 sired properties. Therefore, superconducting qubits are 1. Representation of the readout signal 43 also often referred to as artificial atoms, offering a rich 2. I-Q mixing 44 parameter space of possible qubit properties and opera- 3. Homodyne demodulation 44 tion regimes, with predictable performance in terms of transition frequencies, anharmonicity, and complexity. 4. Heterodyne demodulation 45 While there are many other excellent reviews on su- C. Weak and strong qubit measurements: perconducting qubits, see e.g. Refs. 34–43, this work Impact of noise 46 specifically aims to introduce new quantum engineers D. “Purcell filters” for faster readout 48 (academic and industrial alike) to the terminology and E. Improve signal-to-noise ratio: Parametric state-of-the-art practices used in the rapidly accelerat- amplification 50 ing field of superconducting quantum computing. The 1. Quantum-limited amplification processes 50 reader is assumed to be familiar with basic concepts that 2. Operation of Josephson parametric span classical physics, quantum mechanics, and electrical amplifiers 51 engineering. In particular, readers will find it useful to 3. The traveling wave parametric amplifier 53 have had previous exposure to classical mechanics, the Schr¨odingerequation, the Bloch sphere representation of VI. Summary and outlook 54 qubit states, second quantization, basic concepts of su- perconductivity, electromagnetism, introductory circuit Acknowledgments 55 analysis, classical boolean logic, linear dynamical sys- tems, analog and digital signal processing, and famil- iarity with microwave components such as transmission I. INTRODUCTION lines and mixers. These topics will be introduced as they arise, but having basic prior knowledge will be helpful. Quantum processors harness the intrinsic properties of quantum mechanical systems – such as quantum paral- lelism and quantum interference – to solve certain prob- A. Organization of this article lems where classical computers fall short1–6. Over the past two decades, rapid developments in the science and This review is organized in the following four sections; engineering of quantum systems have advanced the fron- first, in Sec. II, we explore the parameter space avail- tier in quantum computation, from the realm of scien- able when designing superconducting circuits. In particu- tific explorations on single isolated quantum systems to- lar, we look at the promising capacitively-shunted planar ward the creation and manipulation of multi-qubit pro- qubit modalities and how these can be engineered with cessors7,8. In particular, the requirements imposed by desired properties, such as transition frequency, anhar- larger quantum processors have shifted of mindset within monicity, and reduced susceptibility to various sources the community, from solely scientific discovery to the de- of noise. In this section, we also introduce several ways velopment of new, foundational engineering abstractions in which interactions between qubits can be engineered, associated with the design, control, and readout of multi- in order to implement two-qubit entangling operations, qubit quantum systems. The result is the emergence of a needed for a universal gate set. new discipline termed quantum engineering, which serves In Sec. III, we discuss systematic and stochastic noise, to bridge the basic sciences, mathematics, and computer the concepts of noise strength and qubit noise suscepti- science with fields generally associated with traditional bility, and the common sources of noise which lead to engineering. decoherence in superconducting circuits. We introduce One prominent platform for constructing a multi-qubit the Bloch-Redfield model of decoherence, characterized quantum processor involves superconducting qubits, in by longitudinal and transverse relaxation times T1 and which information is stored in the quantum degrees of T2, and discuss the implications of 1/f noise. We then freedom of nanofabricated, anharmonic oscillators con- define the noise power spectral density, which is com- structed from superconducting circuit elements. In con- monly used to characterize noise processes, and describe trast to other platforms, e.g. electron spins in sili- how it drives decoherence. Finally, we close the section con9–14 and quantum dots15–18, trapped ions19–23, ul- with a review of coherent control methods used to miti- tracold atoms24–27, nitrogen-vacancies in diamonds28,29, gate certain types of coherence, reversible noise. 3
In Sec. IV, we provide a review of how single- and The time-independent Hamiltonian Hˆ governs the time ˆ two-qubit operations are typically implemented in super- evolution of the system through the operator e−iHt/~. conducing circuits, by using a combination of local mag- Thus, just as with classical systems, determining the netic flux control and microwave drives. In particular, we Hamiltonian of a system – whether the classical Hamil- discuss the family of two-qubit gates arising from a ca- tonian H or its quantum counterpart Hˆ – is the first step pacitive coupling between qubits, and introduce several to deriving its dynamical behavior. In Sec. IV, we con- recent advances that have been demonstrated to achieve sider the case when the Hamiltonian is time-dependent high-fidelity gates, as well as applications in quantum in the context of qubit control. information processing that use these gates. The contin- To understand the dynamics of a superconducting ued development of high-fidelity two-qubit gates in super- qubit circuit, it is natural to start with the classical de- conducting qubits is a highly active research area. For scription of a linear LC resonant circuit [Fig. 1(a)]. In this reason, we include sufficient technical details that a this system, energy oscillates between electrical energy reader may use this review as a starting point to critically in the capacitor C and magnetic energy in the inductor assess the pros and cons of the various gates, as well as L. In the following, we will arbitrarily associate the elec- develop an appreciation for the types of gate-engineering trical energy with the “kinetic energy” and the magnetic already implemented in state-of-the-art superconducting energy with the “potential energy” of the oscillator. The quantum processors. instantaneous, time-dependent energy in each element is Finally, in Sec. V, we discuss the physics and engi- derived from its current and voltage, neering associated with the dispersive readout technique, Z t typically used to measure the individual qubit states in E(t) = V (t0)I(t0)dt0, (3) modern quantum processors. After a discussion of the −∞ theory behind dispersive coupling, we give an introduc- 0 0 tion to design of Purcell filters and the development of where V (t ) and I(t ) denote the voltage and current of quantum-limited parametric amplifiers. the capacitor or inductor. To derive the classical Hamiltonian, we follow the stan- dard approach used in classical mechanics: the Lagrange- II. ENGINEERING QUANTUM CIRCUITS Hamilton formulation. Here, we represent the circuit el- ements in terms of one of its generalized circuit coor- In this section, we will demonstrate how quantum sys- dinates, charge or flux. In the following, we pick flux, tems based on superconducting circuits can be engineered defined as the time integral of the voltage to achieve certain desired properties. Using the most Z t common qubit modalities, we discuss how properties such 0 0 as the qubit transition frequency, anharmonicity, and Φ(t) = V (t )dt . (4) −∞ noise susceptibility can be tailored by the choice of circuit topology and element parameter values. We also discuss In this example, the voltage at the node is also the branch how to engineer the interactions between different quan- voltage across the element. In this section, we will simply tum systems, in particular the cases of qubit-qubit and refer to these as node voltages and fluxes for convenience. qubit-resonator couplings. For a more detailed discussion of nodes and branches in this context, we refer the reader to Ref. 44. Note that in the following, we could have exchanged A. From quantum harmonic oscillator to the transmon our associations with kinetic energy (momentum coor- qubit dinate) and potential energy (position coordinate), and instead start with the charge variable Q(t), which is the A quantum mechanical system is governed by the time- time integral of the current I(t). dependent Schr¨odingerequation, By combining Eqs. (3) and (4), using the relations V = L dI/dt and I = C dV/dt, and applying the integration ∂ Hˆ |ψ(t)i = i~ |ψ(t)i, (1) by parts formula, we can write down energy terms for the ∂t capacitor and inductor in terms of the node flux, where |ψ(t)i is the state of the quantum system at time t, is the reduced Planck’s constant h/2π, and Hˆ is the 1 2 ~ TC = CΦ˙ , (5) Hamiltonian that describes the total energy of the sys- 2 tem. The “hat” is used to indicate that Hˆ is a quantum operator. As the Schr¨odinger equation is a first-order lin- 1 2 ear differential equation, the temporal dynamics of the UL = Φ . (6) L quantum system may be viewed as a straightforward ex- 2 ample of a linear dynamical system with formal solution, The Lagrangian is defined as the difference between the kinetic and potential energy terms and can thus be ˆ |ψ(t)i = e−iHt/~|ψ(0)i. (2) expressed in terms of Eqs. (5) and (6) 4
across the inductor. These two operators form a canon- ical conjugate pair, obeying the commutation relation 1 2 1 2 L = TC − UL = CΦ˙ − Φ . (7) [φ, n] = i. We note that the factor 4 in front of the 2 2L charging energy EC is solely a historical artifact, namely, From the Lagrangian in Eq. (7), we can further derive that this energy scale was first defined for single-electron the Hamiltonian using the Legendre transformation, for systems and then adopted to two-electron Cooper-pair which we need to calculate the momentum conjugate to systems. the flux, which in this case is the charge on the capacitor The Hamiltonian in Eq. (13) is identical to the one de- scribing a particle in a one-dimensional quadratic poten- tial, a quantum harmonic oscillator (QHO). We can treat ∂L Q = = CΦ˙ . (8) φ as the generalized position coordinate, so that the first ∂Φ˙ term is the kinetic energy and the second term is the po- tential energy. We emphasize that the functional form The Hamiltonian of the system is now defined as of the potential energy influences the eigensolutions. For 2 example, the fact that this term is quadratic (UL ∝ φ ) Q2 Φ2 1 1 in Eq. (13) gives rise to the shape of the potential in Fig. H = QΦ˙ − L = + ≡ CV 2 + LI2, (9) 2C 2L 2 2 1(b). The solution to this eigenvalue problem gives an in- finite series of eigenstates |ki,(k = 0, 1, 2,...), whose cor- as one would expect for an electrical LC circuit. Note responding eigenenergies Ek are all√ equidistantly spaced,√ that this Hamiltonian is analogous to that of a mechan- i.e. Ek+1 −Ek = ~ωr, where ωr = 8ELEC /~ = 1/ LC ical harmonic oscillator,√ with mass m = C, and reso- denotes the resonant frequency of the system, see Fig. nant frequency ω = 1/ LC, which expressed in posi- 1(b). We may represent these results in a more compact tion, x, and momentum, p, coordinates takes the form form (second quantization) for the quantum harmonic 2 2 2 H = p /2m + mω x /2. oscillator (QHO) Hamiltonian The Hamiltonian described above is classical. In order to proceed to a quantum-mechanical description of the system, we need to promote the charge and flux coordi- † 1 H = ~ωr a a + , (14) nates to quantum operators. And, whereas the classical 2 coordinates satisfy the Poisson bracket: where a†(a) is the creation (annihilation) operator of a δf δg δg δf single excitation of the resonator. The Hamiltonian in {f, g} = − (10) δΦ δQ δΦ δQ Eq. (14) is written as an energy. It is, however, often δΦ δQ δQ δΦ preferred to divide by ~ so that the expression has units → {Φ,Q} = − = 1 − 0 = 1, (11) of radian frequency, since we will later resonantly drive δΦ δQ δΦ δQ transitions at a particular frequency or reference the rate the quantum operators similarly satisfy a commutation at which two systems interact with one another. There- relation: fore, from here on, ~ will be omitted. The original charge number and phase operators can † † [Φˆ, Qˆ] = ΦˆQˆ − QˆΦˆ = i~, (12) be expressed as n = nzpf×i(a−a ) and φ = φzpf×(a+a ), 1/4 1/4 where nzpf = [EL/(32EC )] and φzpf = (2EC /EL) where the operators are indicated by hats. From this are the zero-point fluctuations of the charge and phase point forward, however, the hats on operators will be variables, respectively. Quantum mechanically, the quan- omitted for simplicity. tum states are represented as wavefunctions that are gen- In a simple LC resonant circuit [Fig. 1(a)], both the erally distributed over a range of values of n and φ and, inductor L and the capacitor C are linear circuit ele- consequently, the wavefunctions have non-zero standard ments. Defining the reduced flux φ ≡ 2πΦ/Φ0 and the deviations. Such wavefunction distributions are referred reduced charge n = Q/2e, we can write down the follow- to as “quantum fluctuations,” and they exist, even in the ing quantum-mechanical Hamiltonian for the circuit, ground state, where they are called “zero-point fluctua- tions”.
2 1 2 The linear characteristics of the QHO has a natural H = 4EC n + ELφ , (13) 2 limitation in its applications for processing quantum in- formation. Before the system can be used as a qubit, we 2 where EC = e /(2C) is the charging energy required need to be able to define a computational subspace con- to add each electron of the Cooper-pair to the island sisting of only two energy states (usually the two-lowest 2 and EL = (Φ0/2π) /L is the inductive energy, where energy eigenstates) in between which transitions can be Φ0 = h/(2e) is the superconducting magnetic flux quan- driven without also exciting other levels in the system. tum. Moreover, the quantum operator n is the ex- Since many gate operations, such as single-qubit gates cess number of Cooper-pairs on the island, and φ – the (Sec. IV), depend on frequency selectivity, the equidis- reduced flux – is denoted the “gauge-invariant phase” tant level-spacing of the QHO, illustrated in Fig. 1(b), 5
(a) i (c)
+ L ~ dφ L C v J C I = Ic sin(φ),V = , (15) r r C s e dt - J 2 resulting in a modified Hamiltonian (b) (d) 5 5 Transmon 4 4 2 H = 4EC n − EJ cos(φ), (16) 3 3 2 where EC = e /(2CΣ), CΣ = Cs + CJ is the total ca- 2 2 pacitance, including both shunt capacitance Cs and the Energy [ ] [ Energy 1 ] [ Energy 1 self-capacitance of the junction CJ , and EJ = IcΦ0/2π is
QHO Comp. the Josephson energy, with Ic being the critical current 0 0 subspace ‡ - - /2 0 /2 - - /2 0 /2 of the junction . After introducing the Josephson junc- Superconducting phase, Superconducting phase, tion in the circuit, the potential energy no longer takes a manifestly parabolic form (from which the harmonic FIG. 1. (a) Circuit for a parallel LC-oscillator (quantum har- spectrum originates), but rather features a cosinusoidal monic oscillator, QHO), with inductance L in parallel with form, see the second term in Eq. (16), which makes the capacitance, C. The superconducting phase on the island is energy spectrum non-degenerate. Therefore, the Joseph- denoted φ, referencing ground as zero. (b) Energy potential son junction is the key ingredient that makes the oscilla- for the QHO, where energy levels are equidistantly spaced tor anharmonic and thus allows us to identify a uniquely ~ωr apart. (c) Josephson qubit circuit, where the nonlinear addressable quantum two-level system, see Fig. 1(d). inductance LJ (represented with the Josephson-subcircuit in Once the nonlinearity has been added, the system dy- the dashed orange box) is shunted by a capacitance, Cs. (d) namics is governed by the dominant energy in Eq. (16), The Josephson inductance reshapes the quadratic energy po- reflected in the EJ /EC ratio. Over time, the super- tential (dashed red) into sinusoidal (solid blue), which yields non-equidistant energy levels. This allows us to isolate the conducting qubit community has converged towards cir- two lowest energy levels |0i and |1i, forming a computational cuit designs with EJ EC . In the opposite case when EJ ≤ EC , the qubit becomes highly sensitive to charge subspace with an energy separation ~ω01, which is different than ~ω12. noise, which has proven more challenging to mitigate than flux noise, making it very hard to achieve high co- herence. Another motivation is that current technologies allow for more flexibility in engineering the inductive (or † poses a practical limitation . potential) part of the Hamiltonian. Therefore, working To mitigate the problem of unwanted dynamics in- in the EJ ≤ EC limit, makes the system more sensitive volving non-computational states, we need to add anhar- to the change in the potential Hamiltonian. Therefore, monicity (or nonlinearity) into our system. In short, we we will focus here on the state-of-the-art qubit modalities 0→1 1→2 require the transition frequencies ωq and ωq be suffi- that fall in the regime EJ EC . For readers who are ciently different to be individually adressable. In general, interested in the physics in the EJ ≤ EC regime, such the larger the anharmonicity the better. In practise, the as the earlier Cooper-pair box charge qubit, we refer to amount of anharmonicity sets a limit on how short the Refs. 48–51. pulses used to drive the qubit can be. This is discussed To access the EJ EC regime, one preferred approach in detail in Sec. IV D 3. is to make the charging EC small by shunting the junction To introduce the nonlinearity required to modify the with a large capacitor, Cs CJ , effectively making the harmonic potential, we use the Josephson junction – a qubit less sensitive to charge noise – a circuit commonly nonlinear, dissipationless circuit element that forms the known as the transmon qubit52. In this limit, the super- 46,47 backbone in superconducting circuits . By replacing conducting phase φ is a good quantum number, i.e. the the linear inductor of the QHO with a Josephson junc- spread (or quantum fluctuation) of φ values represented tion, playing the role of a nonlinear inductor, we can by the quantum wavefunction is small. The low-energy modify the functional form of the potential energy. The eigenstates are therefore, to a good approximation, local- potential energy of the Josephson junction can be derived ized states in the potential well, see Fig. 1(d). We may from Eq. (3) and the two Josephson relations gain more insight by expanding the potential term of Eq. (16) into a power series (since φ is small), that is
†Even though linear resonant systems cannot be addressed properly, their long coherence times have proven them useful as quantum ac- cess memories for storing quantum information, where a nonlinear ‡The critical current is the maximum supercurrent that the junction ancilla system is used as a quantum controller for feeding and ex- can support before it switches to the resistive state with non-zero tracting excitations to/from the resonant cavity modes45. voltage. 6
B. Qubit Hamiltonian engineering
1 2 1 4 6 EJ cos(φ) = EJ φ − EJ φ + O(φ ). (17) 1. Tunable qubit: split transmon 2 24
The leading quadratic term in Eq. (17) alone will To implement fast gate operations with high-fidelity, result in a QHO, recall Eq. (13). The second term, as needed to implement quantum logic, many (though however, is quartic which modifies the eigensolution and not all63) of the quantum processor architectures imple- disrupts the otherwise harmonic energy structure. Note mented today feature tunable qubit frequencies64–67. For that, the negative coefficient of the quartic term indicates instance, in some cases, we need to bring two qubits into 1→2 0→1 that the anharmonicity α = ωq −ωq is negative and resonance to exchange (swap) energy, while we also need its limit in magnitude thus cannot be made arbitrarily the capability of separating them during idling periods large. For the case of the transmon, α = −EC is usually to minimize their interactions. To do this, we need an designed to be 100 − 300 MHz, as√ required to maintain external parameter which allows us to access one of the a desirable qubit frequency ωq = ( 8EJ EC − EC )/~ = degrees of freedom of the system in a controllable fashion. 3−6 GHz, while keeping an energy ratio sufficiently large 52 One widely-used technique is to replace the single (EJ /EC ≥ 50) to suppress charge sensitivity . Fortu- Josephson junction with a loop interupted by two iden- nately, the charge sensitivity is exponentially suppressed tical junctions – forming a dc superconducting quantum for increased EJ /EC , while the reduction in anharmonic- interference device (dc-SQUID)68. Due to the interfer- ity only scales as a weak power law, leading to a workable ence between the two arms of the SQUID, the effective device. critical current of the two parallel junctions can be de- Including terms up to fourth order and using the QHO creased by applying a magnetic flux threading the loop, eigenbases, the system Hamiltonian resembles that of a see Fig. 2(a). Due to the fluxoid quantization condition, Duffing oscillator the algebraic sum of branch flux of all of the inductive elements along the loop plus the externally applied flux † α † † equal an integer number of superconducting flux quanta, H = ωqa a + a a aa. (18) 2 that is
Since |α| ωq, we can see that the transmon qubit is basically a weakly anharmonic oscillator (AHO). If exci- ϕ1 − ϕ2 + 2ϕe = 2πk, (20) tation to higher non-computational states is suppressed over any gate operations, either due to a large enough |α| where ϕe = πΦext/Φ0. Using this condition, we can elim- or due to robust control techniques such as the DRAG inate one degree of freedom and treat the SQUID-loop pulse, see Sec. IV D 3, we may effectively treat the AHO as a single junction, but with the important modification as a quantum two-level system, simplifying the Hamilto- that EJ is tunable (via the SQUID critical current) by nian to means of the external flux Φext. The effective Hamilto- nian of the so-called split transmon (ignoring the con- σz stant) is H = ωq , (19) 2 where σz is the Pauli-z operator. However, one should 2 H = 4EC n − 2EJ |cos (ϕe)| cos(φ). (21) always keep in mind that the higher levels physically | {z } 53 E0 (ϕ ) exist . Their influence on system dynamics should be J e taken into account when designing the system and its control processes. In fact, there are many cases where We can see that Eq. (21) is analogous to Eq. (16), with 0 the higher levels have proven useful to implement more EJ replaced by EJ (ϕe) = 2EJ |cos (ϕe)|. The magnitude 54 0 efficient gate operations . of the net, effective Josephson energy EJ has a period In addition to reducing the charge dispersion, the use of Φ0 in applied flux and spans from 0 to its maximum of a large shunt capacitor also enables us to engineer value 2EJ . Therefore, the qubit frequency can be tuned the electric field distribution of the quantum system, and periodically with Φext, see Fig. 2(b). thus the participation of surface loss mechanisms. In the While the split transmon enables frequency tunabil- development of the 3D transmon55, e.g. a 2D transmon ity by the externally applied magnetic field, it also in- coupled to a 3D cavity, it was demonstrated that by mak- troduces sensitivity to random flux fluctuations, known ing the gap between the two lateral capacitor plates large as flux noise. At any working point, the slope of the (compared to the film thickness) the coherence time in- qubit spectrum, ∂ωq/∂Φext, indicates to first order how creases since a smaller portion of the electric field in- strongly this flux noise affects the qubit frequency. The teracts with the lossy interfaces, e.g. metal-substrate sensitivity is generally non-zero, except at multiples of and substrate-vacuum interfaces, which has been stud- the flux quantum, Φext = kΦ0, where k is an integer, 56–61 ied extensively . where ∂ωq/∂Φext = 0. 7
(a) Symmetric transmon (c) Asymmetric transmon(e) C-shunted Flux qubit(g) C-shunted Fluxonium
φ1 φ1 φ2 φ1 φ2 φ1 φ2 CJ Ic Ic Ic φ2 φN
(b) (d) (f) (h) 8 8 8 8 Ec /h = 0.3 GHz EJ /h = 15 GHz 6 6 6 6
4 4 4 4
2 2 2 2
0 0 0 0 - - /2 0 /2 - - /2 0 /2 - - /2 0 /2 - - /2 0 /2
FIG. 2. Modular qubit circuit representations for capacitively shunted qubit modalities (orange box Fig. 1c) and corresponding qubit transition frequencies for the two lowest energy states as a function of applied magnetic flux in units of Φ0. (a-b) Symmetric transmon qubit, with Josephson energy EJ are shunted with a capacitor yielding a charging energy EC . (c-d) Asymmetric transmon qubit, with junction asymmetry γ = EJ2/EJ1 = 2.5. (e-f) Capacitively shunted flux qubit, where a small principle junction (red) is shunted with two larger junctions (orange). Parameters are the same as Yan et al.62. (g-h) C-shunted fluxonium qubit, where the small junction is inductively shunted with a large array of N junctions.
One recent development has focused on reducing the resonance gate is optimized with certain frequency de- qubit sensitivity to flux noise, while maintaining suf- tuning between two qubits70. Therefore, by using an ficient tunability to operate our quantum gates. The asymmetric transmon, a small frequency-tuning range idea is to make the two junctions in the split transmon is introduced that is sufficient to compensate for fabri- asymmetric69, see Fig. 2(c). This yields the following cation variations, without introducing unnecessary large Hamiltonian susceptibility to flux noise and thus maintaining high co- herence. For another example, a surface code scheme q based on the adiabatic CPHASE-gate requires specific fre- 2 2 2 2 H = 4EC n − EJΣ cos (ϕe) + d sin (ϕe) cos(φ), (22) quency configuration among qubits in order to avoid fre- | {z } quency crowding issues, and asymmetric transmons fit E0 (ϕ ) J e well with its well-defined frequency range71. In general, where EJΣ = EJ1 + EJ2 and d = (γ − 1)/(γ + 1) is as the quantum processors scale up and fabrication im- the junction asymmetry parameter, with γ = EJ2/EJ1. proves, asymmetric transmons are likely to be found in Again, we can treat the two junctions as a single-junction wider applications in the future. 0 transmon, with an effective Josephson energy EJ (ϕe). In particular, we can recognize the two special cases; for d = 0, the Hamiltonian in Eq. (22) reduces to the sym- 2. Towards larger anharmonicity: flux qubit and fluxonium 0 metric case with EJ (ϕe) = EJΣ |cos(ϕe)|, as in Eq. (21) with EJΣ = 2EJ . In the other limit, when |d| → 1, We see that split transmon qubits, be it symmetric or 0 EJ (ϕe) → EJΣ and the flux-tunability of the Josephson not, still share the same topology as the single junction energy vanishes, which is equivalent to the single junction version, yielding a sinusoidal potential. Therefore, the case, recall Eq. (16). degree to which the properties of these qubits can be en- From the discussion above we see that going from sym- gineered has not fundamentally changed. In particular, metric to asymmetric transmons does not change the cir- the limited anharmonicity in transmon-type qubits in- cuit topology. This seemingly trivial modification, how- trinsically causes significant residual excitation to higher- ever, has profound impact for practical applications. As energy states, undermining performance of gate opera- we can see from the qubit spectra, Fig. 2(d), the flux tions. To go beyond this, it is necessary to introduce sensitivity is suppressed across the entire tunable fre- additional complexity into the circuit. quency range. For example, the performance of the cross- One outstanding development in this regard is the in- 8 vention of the flux qubit72,73, where the qubit loop is in- the introduction of the capacitive shunt, similar to the terrupted by three (or four) junctions, see Fig. 2(e). On modified Cooper-pair box leading to the transmon qubit. one branch is one smaller junction; on the other branch The double-well case obtained for γ < 2 was demon- are two identical junctions, both a factor γ larger in size strated and investigated much earlier72,73. The intuitive compared to the small junction. The addition of one picture based on circulating current states – so it gets more junction as compared to the split transmon is non- the name persisting-current flux qubit (PCFQ) – gives trivial, as it changes the circuit topology and reshapes a satisfying physical description of the qubit degrees of the potential energy profile. freedom. However, from the perspective of a quantum Each junction is associated with a phase variable, and engineer, the qubit properties are of more interest, even the fluxoid quantization condition again allows us to if sometimes we may lose physical intuition about the eliminate one degree of freedom. Consequently, we have system in certain regimes; such as when γ ≈ 2 and there a two-dimensional potential landscape, which in compar- are no clear circulating current states. The most impor- ison to the simpler topology of the transmon, compli- tant feature of the PCFQ is that its anharmonicity can cates the problem both conceptually and computation- be much greater than the transmon and CSFQ and the ally. Fortunately, under the assumed setting that the ar- transition matrix elements |h1|nˆ|0i|, |h1|φˆ|0i| become con- ray junctions are larger in size (γ > 1), it is usually a good siderably smaller given equivalent EJ /EC . Therefore, a approximation to treat the problem as a particle moving longer relaxation time can be expected. These features in a quasi-1D potential, which also helps us gain more have been demonstrated even more prominently in its insight and intuition about the system and draw qualita- close relative, the fluxonium qubit75. tive conclusions. The Hamiltonian under this quasi-1D The flux qubit is a striking example that illustrates approximation reads, how one dramatically can engineer the qubit properties through the choice of various circuit parameters. The in- troduction of array junctions and consequent biharmonic 2 H ≈ 4EC n − EJ cos(2φ + ϕe) − 2γEJ cos(φ). (23) profile generates rich dynamics as well as broad applica- tions. An extention of this idea is the fluxonium qubit, Note that the phase variable in Eq. (23) is the sum which generated substantial interest recently, due partly of the branch phases across the two array junctions, φ = to its capability of engineering the transition matrix el- ements to achieve millisecond T1 time, and due partly (ϕ1 + ϕ2)/2, assuming the same current direction across to the invention of novel gate schemes applicable to such ϕ1 and ϕ2. The external magnetic flux is denoted ϕe = well-protected qubits76,77. 2πΦext/Φ0. The second term in Eq. (23) is contributed Compared to flux qubits, which usually contain two by the small junction with Josephson energy EJ , whereas 78 the third term takes into account the two array junctions, or three array junctions , the number of array junctions in the fluxonium qubit is dramatically increased75,79, in together with Josephson energy 2γEJ . Clearly, the sum of these two terms no longer has the characteristics of a some cases, to the order of 100, see Fig. 2(g). Following simple cosinusoid, and the final potential profile as well the same quasi-1D approximation as for the flux qubit, −NγE φ/N as the corresponding eigenstates depends on both the the last term in Eq. (23) becomes J cos( ), where N denotes the number of array junctions. For external flux ϕe and the junction area ratio γ. large N, the argument in the cosine term φ/N becomes The most common working point for this system is sufficiently small that a second order expansion is a good when ϕe = π + 2πk, where k is an integer – that is approximation. This results in the fluxonium Hamilto- when half a superconducting flux quantum threads the nian, qubit loop. At this flux bias point, the qubit spectrum reaches its minimum, and the qubit frequency is first- 2 1 2 order insensitive to flux noise, see Fig. 2(f). This point H ≈ 4EC n − EJ cos(φ + ϕe) + ELφ , (24) is often referred to as the flux degeneracy point, where 2 flux qubits tend to have the optimal coherence time. where EL = (γ/N)EJ is the inductive energy of the At this operation point, the potential energy may as- effective inductance contributed by the junction ar- sume a single-well (γ ≥ 2) or a double-well (γ < 2) ray – often known as superinductance due to its large profile. The single-well case shares some simularities value79–81. Therefore, we can treat the potential energy with the transmon qubit, where the quadratic and quar- as a quadratic term modulated by a sinusoidal term, sim- tic terms of the Hamiltonian determines the harmonic- ilar to that of an rf-SQUID type flux qubit82. However, ity and anharmonicity, respectively. The capacitively- the kinetic inductance of the Josephson junction array is shunted flux qubit (CSFQ)62,74 was explored in this in general much larger than the geometric inductance of regime, demonstrating long coherence and decently high the wire in an rf-SQUID. anharmonicity. Note that as opposed to the transmon Depending on the relative magnitude of EJ and EL, qubit, the anharmonicity of the CSFQ is positive (α > 0). the fluxonium system could involve plasmon states (in While the improvement in anharmonicity can be associ- the same well) and fluxon states (in different wells). ated with reshaping the energy potential, the improved There are a variety of schemes to utilize them for quan- coherence over the first flux qubits can be attributed to tum information processing. Generally, the spectrum of 9 the transition between the lowest energy states is similar (a) Direct capacitive coupling to that of the flux qubit, see Fig. 2(h). Both long coher- g12 ence and high anharmonicity can be expected at the flux V1 V2 sweet spot. Cg Lastly, we want to point out a further extension – the 0 − π qubit – which has even stronger topological protec- IC1 C1 IC2 C2 tion from noise83,84. However, the strongly suppressed sensitivity to external fluctuations also makes it hard to manipulate. (b) Capacitive coupling via coupler gr1 gr2 C. Interaction Hamiltonian engineering Cg1 Cg2 To generate entanglement between individual quantum IC1 C1 Lr Cr IC2 C2 systems – it is necessary to engineer an interaction Hamil- tonian that connects degrees of freedom in those indi- vidual systems. In this section, we discuss the physical coupling mechanism and its representation in the qubit (c) Direct inductive coupling eigenbasis. The use of coupling to form 2-qubit gates is I1 M12 I2 discussed in Sec. IV.
IC1 L1 L2 IC2 1. Physical coupling: capacitive and inductive Φe1 Φe2
The Hamiltonian of two coupled systems takes a generic form (d) Inductive coupling via coupler ICC M1C M2C
H = H1 + H2 + Hint, (25) IC1 L1 L2 IC2 Φe1 ΦeC Φe2 where H1 and H2 denote the Hamiltonians of the individ- ual quantum systems, which could be any combination of the qubit circuits mentioned in Sec. II A and II B. The H FIG. 3. Schematic of capacitive and inductive coupling last term, int, is the interaction Hamiltonian, which schemes between two superconducting qubits, labeled 1 and couples variables of both systems. In superconducting 2. (a) Direct capacitive coupling, where the voltage nodes circuits, the physical form of the coupling energy is either of two qubits V1 and V2 are connected by a capacitance Cg. an electric or magnetic field (or a combination thereof). (b) Capacitive coupling via a coupler in form of a linear res- To achieve capacitive coupling, a capacitor is placed onator. (c) Direct inductive coupling, where the two qubits between the voltage nodes of the two participating cir- are coupled via mutual inductance, M12. (d) Inductive cou- cuits, yielding an interaction Hamiltonian Hint of the pling via mutual inductances M1C and M2C to a frequency- form tunable coupler.
Hint = CgV1V2, (26) where the expressions in brackets are the two Hamiltoni- C V V where g is the coupling capacitance and 1( 2) is the ans of the individual qubits, [see Eq. (16)], and we take voltage operator of the corresponding voltage node being Vi = (2e/Ci)ni in Eq. (26). From Eq. (27), we see that connected. Fig. 3(a) illustrates a realistic example of a the coupling energy depends on the coupling capacitance direct capacitive coupling between the top nodes of two as well as the matrix elements of the voltage operators. transmon qubits. Circuit quantization in the limit of The dependencies are bilinear in the perturbative limit Cg C1,C2 yields (Cg C1,C2). To implement the coupling capacitance, one only need X 2 bring the edges of the capacitor pads into close proxim- H = 4EC,ini − EJ,i cos(φi) i=1,2 ity, as has been demonstrated in state-of-the-art planar designs85. The coupling capacitance is determined by the planar capacitor geometry as well as the surround- 2 Cg ing environment, such as the dielectric constant of the +4e n1n2, (27) C1C2 substrate and the ground plane proximity. 10
In the case of inductive coupling, a mutual inductance this simplifying approximation is only exact in the patho- shared by two loops is the coupling mechanism, yielding logical limit of no coupling. an interaction Hamiltonian that is of the intuitive form To realize a mutual inductance, two looped circuits are brought into close proximity to one another, or, to Hint = M12I1I2, (28) make it stronger, overlap with each other88, and even may share the same wire or Josephson junction induc- M 89–92 where 12 denotes the mutual inductance between the tor . In the case of a Josephson junction, and for I I loops and 1 and 2 are the current operators for the certain metals, the inductance is dominated by kinetic currents through the inductors. A typical example com- inductance contributions, rather than solely geometric prises two closely positioned (rf-SQUID type) flux qubits, inductance93,94. Kinetic inductance arises from the me- as illustrated in Fig. 3(c). The system Hamiltonian can chanical, inertial mass of the charge carriers, but is only be expressed as (see Refs. 86 and 87): practically witnessed in very high-conductance materi- als like superconductors. A primary feature of kinetic 2 inductance is that its values can vastly exceed those of X 2 1 ΦLi H = 4EC,ini − EJ,i cos(φi) + conventional geometric inductances, which are generally L − K2 79 i=1,2 2 i(1 ) limited by electromagnetic considerations .
2 ΦL1 ΦL2 − M12(1 − K ) , (29) L − K2 L − K2 1(1 ) 2(1 ) 2. Coupling axis: transverse and longitudinal where the first two terms are the energies associated with the Josephson junctions, the third term captures the in- Regardless of its physical realization, the effect of a ductor energies, and the fourth term is the mutual cou- coupling on system dynamics is determined by its form pling energy of the form MI1I2 (note that in general as represented in the eigenbasis of the individual systems. Φ = LI). ΦL1 and ΦL2 are the magnetic fluxes asso- That is, how Hint appears in the representation spanned ciated with the currents flowing through the respective by the eigenbasis of H1 ⊗ H2. inductors, and K is the√ unitless mutual coupling factor Let us start with the previous example of two capaci- defined by M12 = K L1L2. Importantly, note that in tively coupled transmon qubits [Fig. 3(a)]. Using second 2 Eq. (29), K renormalizes L1, L2 and M12, essentially quantization, the system Hamiltonian in Eq. (27) can be capturing the loading effect on the circuit due to their expressed as presence, and is found by inverting an inductance matrix (see Refs. 86 and 87). X h † αi † † i H = ωia ai + a a aiai As we did with the charge degree of freedom, we will i 2 i i normalize the inductor and externally applied magnetic i∈1,2 fluxes – in this case, by the reduced quantum unit of flux Φ0/2π – to define phases φ ≡ Φ/(Φ0/2π) that are † † −g a1 − a a2 − a , (31) on the same footing as the junction phases φ1,2. Due 1 2 to fluxoid quantization around the closed loop, these phases must sum to zero or an integer multiple of 2π. where the expression within brackets represent the Duff- g For current directions entering the top of the junctions ing oscillator Hamiltonian for the qubits and is the coupling energy. Since we define V ∝ n ∝ i(a − a†), and [with I1 counterclockwise and I2 as shown in Fig. 3(c)], † consequently I ∝ φ ∝ (a + a ), the original n1n2-term φi −φLi = φei +2πn with n = 0, ±1,... (see, for example, becomes what is shown in Eq. (31). Such a coupling is Ref. 72). Replacing φL1(L2) with φ1(2) − φ1(e2) yields: called transverse, because the coupling Hamiltonian has 2 non-zero matrix elements only at off-diagonal positions " Φ0 2 # 1 (φi − φei) † X 2 2π with respect to both oscillators, i.e. ihk|ai − ai |kii = 0 = 4EC,ini − EJ,i cos(φi) + 2 2 Li(1 − K ) for any integer k and for i ∈ 1, 2 and in this case i=1,2 † ihk ± 1|ai − a |kii =6 0. Φ0 φ − φ Φ0 φ − φ i 2 2π ( 1 e1) 2π ( 2 e2) If we can ignore higher energy levels (k ≥ 2) either − M12(1 − K ) 2 2 . (30) L1(1 − K ) L2(1 − K ) because of sufficient anharmonicity or through careful control protocols that ensure these levels never have in- M I I The coupling is of the form 1,2 1 2, with both the mu- fluence, we may truncate the Hamiltonian in Eq. (31) tual coupling and the circulating currents (via the L1 2 to and L2) “renormalized” by the factor (1 − K ). To capture the renormalization explicitly, the Hamiltonian is generally simulated using phase operators (Eq. 30), X 1 H = ωiσz,i + gσy,1σy,2. (32) rather than current operators. In the weak-coupling limit i∈1,2 2 2 2 K 1 (equivalently, M L1L2), the coupling term may be approximated by the Josephson currents with This is a Hamiltonian of two spins, coupled by an ex- M12I1I2 ≈ M12Ic1 sin φ1Ic2 sin φ2. Note, however, that change interaction. As we will see in Sec. IV D 1, such 11 a Hamiltonian is most commonly used in contemporary and approximation – be treated as an isolated system, implementations and can generate various types of two- and the composite system simplified to two transversely qubit entangling gates. Note that, more often, we see coupled qubits, see Eq. (32). that the interaction term is expressed in σxσx instead of We now turn to the previous example of two induc- σyσy. The choice in the context here is arbitrary and does tively coupled flux qubits, see Fig. 3(c). Assume that not change the dynamics. However, when both capaci- the double-well potential [Fig. 2(g)] has a relatively tive and inductive couplings are present in the system, high inter-well barrier, which leads to an exponentially both σxσx and σyσy may be needed. In this case, the small qubit transition frequency at the energy degener- † voltage operator V ∝ i(a − a ) (reduced to σy after two- acy point, (Φe = π). Around this degeneracy point, the level approximation in the lab frame) is transversal to the off-diagonal matrix element of sin(φ) is zero, i.e. the † current operator I ∝ (a+a ) (reduced to σx) and both of ground and excited states are localized in different wells them may be transverse to the qubit. A similar example and h1| sin(φ)|1i − h0| sin(φ)|0i= 6 0. We can then rewrite is demonstrated between a qubit and a resonator by Lu the Hamiltonian in Eq. (??) as et al.95 Transverse coupling can be engineered between a qubit X 1 and a harmonic oscillator, see Fig. 3(b). In this case, the H = ωiσzi + gσz1σz2. (35) Hamiltonian becomes i=1,2 2
Now, the coupling axis is the same as the qubit quan- 1 † † H = ωqσz + ωra a + g(σ+a + σ−a ), (33) tization axes and therefore termed longitudinal coupling. 2 Note, however, that the physical σxσx and σzσz couplings where ωq and ωr denote the qubit and resonator fre- can change in the qubit frame. quencies, and σ+ = |0ih1| and σ− = |1ih0| describes Longitudinal coupling is an important type of inter- the processes of exciting and de-exciting the qubit, re- action, because it can generate entanglement without spectively. Here, we have assumed that the coupling is energy exchange. Moreover, it is found a necessary in- in the dispersive limit, i.e. g ωq, ωr, hence ignoring gredient in the application of quantum annealing, where † the double (de)excitation terms proportional to σ+a and certain hard combinatorial optimization problems can be σ−a, which under typical operation regimes oscillate suf- modeled by the Ising Hamiltonian in Eq. (35) and finding ficiently fast to average to zero. The Hamiltonian in Eq. its ground state would solve this problem. (33), is the standard model used for describing how a two- An intermediate qubit mode may also be used as a cou- level atom interacts with a resonant cavity that houses pler in the longitudinal case. In Fig. 3(d), an additional it. Such a structure is also known as cavity quantum rf-SQUID is used to mediate the coupling. The coupling electrodynamics (cQED), and it is extended to the cir- strength can be tuned by the flux bias of the coupler cuit version here. It has many useful applications in su- SQUID103. Note that a tunable coupler may also be re- perconducting quantum information architectures, such alized in a structure with capacitive couplings63. A tun- 96 97 as high-fidelity readout , see Sec. V, cavity buses , able coupler is useful because it provides a wide range of 98,99 quantum memory , quantum computation with cat coupling strengths87, a high on-off ratio104 for reducing 100–102 states , etc. gate error-rates, and more ways of achieving high-fidelity Here, we briefly mention the use of a cavity or res- entangling gates67,105–107. The trade-off is an additional onator to mediate coupling between two qubits, which control line. may be physically well-separated (≈ 1 cm). Since most In addition to the pure transversal and longitudinal superconducting resonators are in the GHz frequency qubit-qubit interactions, there are also examples of mixed range, they can be made much longer than any dimen- types of interaction Hamiltonians108 sion of a qubit circuit (≈ 1 mm). One can use such a resonator to mediate coupling between two or more oth- erwise non-interacting qubits. An example is shown in 1 † † H = ωqσz + ωra a + gσz(a + a ), (36) Fig. 3(b), where two transmon qubits are both capac- 2 itively coupled to the center resonator. The two-level system Hamiltonian is: which are longitudinal with respect to a qubit, but trans- verse with respect to a harmonic oscillator in a qubit- resonator system. Such a model is called longitudinal X † αi † † † H = ωiai ai + ai ai aiai + ωrarar but one should note that it is only longitudinal to one i=1,2 2 participating system. It is hard to engineer physically longitudinal coupling with respect to a harmonic oscil- lator, since either the E-field (V ) or the B-field (I) is † † † † +g1r a1ar + a1ar + g2r a2ar + a2ar . (34) transverse with respect to the eigen field of the harmonic oscillator. Note, however, that a transversal model such It can be shown that in the dispersive limit, i.e. gir as in Eq. (33) may be transformed into a longitudinal |ωi −ωr|, the resonator can – after proper transformation one in certain operating regimes, see Sec. V. 12
In some applications, such as for quantum annealing, 2. Stochastic noise both longitudinal and transverse couplings are desired (σzσz coupling for mapping the problem and σxσx cou- The second type of noise is stochastic noise, arising pling for enhancing the annealing performance) and re- from random fluctuations of parameters that are coupled quire independent control. to our qubit109. For example, thermal noise of a 50Ω re- sistor in the control lines leading to the qubit will have voltage and current fluctuations – Johnson noise – with a III. NOISE, DECOHERENCE, AND ERROR noise power that is proportional to both temperature and MITIGATION bandwidth. Or, the oscillator that provides the carrier for a qubit control pulse may have amplitude or phase Random, uncontrollable physical processes in the qubit fluctuations. Additionally, randomly fluctuating electric control and measurement equipment or in the local envi- and magnetic fields in the local qubit environment – e.g., ronment surrounding the quantum processor are sources on the metal surface, on the substrate surface, at the of noise that lead to decoherence and reduce the opera- metal-substrate interface, or inside the substrate – can tional fidelity of the qubits. In this section, we introduce couple to the qubit. This creates unknown and uncon- the basics of noise leading to decoherence in supercon- trolled fluctuations of one or more qubit parameters, and ducting circuits, and we discuss coherent control methods this leads to qubit decoherence. to mitigate certain types of noise.
3. Noise strength and qubit susceptibility A. Types of noise
The degree to which a qubit is affected by noise is In a closed system, the dynamical evolution of a qubit related to the amount of noise impinging on the qubit, state is deterministic. That is, if we know the starting and the qubit’s susceptibility to that noise. The former state of the qubit and its Hamiltonian, then we can pre- is often a question of materials science and fabrication; dict the state of the qubit at any time in the future. that is, can we make devices with lower levels of noise. However, in open systems, the situation changes. The Or, it may be related to the quality of the control elec- qubit now interacts with uncontrolled degrees of freedom tronics and cryogenic engineering to limit the levels of in its environment, which we refer to as fluctuations or noise on the control lines that necessarily connect to the noise. In the presence of noise, as time progresses, the qubits to control them. The latter – qubit susceptibility qubit state looks less and less like the state we would – is a question of qubit design. Qubits can be designed to have predicted and, eventually, the state is lost. There trade off sensitivity to one type of noise at the expense of are many different sources of noise that affect quantum increased sensitivity to other types of noise. Thus, mate- systems, and they can be categorized into two primary rials science, fabrication engineering, electronics design, types: systematic noise and stochastic noise. cryogenic engineering, and qubit design all play a role in creating devices with high coherence. In general, one should strive to eliminate the sources of noise, and then 1. Systematic noise design qubits that are insensitive to the residual noise. The qubit response to noise depends on how the noise Systematic noise arises from a process that is trace- couples to it – either through a longitudinal or a trans- able to a fixed control or readout error. For example, we verse coupling as referenced to the qubit quantization apply a microwave pulse to the qubit that we believe axis. This can be visualized using a Bloch Sphere picture will impart a 180-degree rotation. However, the con- of the qubit state, as illustrated in Fig. 4 and discussed trol field is not tuned properly and, rather than rotating in detail in Section III B. the qubit 180 degrees, the pulse slightly over-rotates or under-rotates the qubit by a fixed amount. The underly- ing error is systematic, and it therefore leads to the same B. Modeling noise and decoherence rotation error each time it is applied. However, when such erroneous pulses are used in practice in a variety of control sequences, the observed results may appear 1. Bloch sphere representation to be influenced by random noise. This is because the pulse is generally not applied in the same way for each The Bloch sphere is a unit sphere used to represent the experiment: it could be applied a different number of quantum state of a two-level system (qubit). Fig. 4(a) times, interspersed with different pulses in different or- shows a Bloch sphere with a Bloch vector representing ders, and therefore generally differs from experiment to the state |ψi = α|0i + β|1i. If we visualize the Bloch experiment. However, once systematic errors are iden- sphere as the planet Earth, then by convention, the north tified, they can generally be corrected through proper pole represents state |0i and the south pole state |1i. For calibration or the use of improved hardware. pure quantum states such as |ψi, the Bloch vector is of 13
(a)Bloch sphere (b)Longitudinal relaxation (c) Pure dephasing (d) Transverse relaxation z z z Longitudinal z (Longitudinal) noise |0〉 |0〉 |0〉 |0〉 |Ψ〉 = α|0〉 + β|1〉 Pure Excitation Relaxation Dephasing Decoherence θ Γ1 Γ1 Γφ Γ2 Γ1 y y y y φ (Transverse) Transverse x x x Γφ x Γφ (Transverse) noise Transverse |1〉 noise |1〉 |1〉 |1〉
FIG. 4. Transverse and longitudinal noise represented on the Bloch sphere. (a) Bloch sphere representation of the quantum state |ψi = α|0i + β|1i. The qubit quantization axis – the z axis – is longitudinal in the qubit frame, corresponding to σz terms in the qubit Hamiltonian. The x-y plane is transverse in the qubit frame, corresponding to σx and σy terms in the qubit Hamiltonian. (b) Longitudinal relaxation results from energy exchange between the qubit and its environment, due to transverse noise that couples to the qubit in the x − y plane and drives transitions |0i ↔ |1i. A qubit in state |1i emits energy to the environment and relaxes to |0i with a rate Γ1↓ (blue arched arrow). Similarly, a qubit in state |0i absorbs energy from the environment, exciting it to |1i with a rate Γ1↑ (orange arched arrow). In the typical operating regime kBT ~ωq, the up-rate is suppressed, leading to the overall decay rate Γ1 ≈ Γ1↓. (c) Pure dephasing in the transverse plane arises from longitudinal noise along the z axis that fluctuates the qubit frequency. A Bloch vector along the x-axis will diffuse clockwise or counterclockwise around the equator due to the stochastic frequency fluctuations, depolarizing the azimuthal phase with a rate Γφ. (d) Transverse relaxation results in a loss of coherence at a rate Γ2 = Γ1/2 + Γφ,√ due to a combination of energy relaxation and pure dephasing. Pure dephasing leads to decoherence of the quantum state (1/ 2)(|0i + |1i), initially pointed along the x-axis. Additionally, the excited state component of the superposition state may relax to the ground state, a phase-breaking process that loses the orientation of the vector in the x-y plane. unit length, |α|2 + |β|2 = 1, connecting the center of the |ψihψ| for a pure state |ψi is equivalently sphere to any point on its surface.
The z-axis connects the north and south poles. It −iφ 1 1 1 + cos θ e sin θ is called the longitudinal axis, since it represents the ρ = (I + ~a · ~σ) = (38) eiφ sin θ 1 + sin θ qubit quantization axis for the states |0i and |1i in the 2 2 2 θ −iφ θ θ qubit eigenbasis. In turn, the x-y plane is the trans- cos 2 e cos 2 sin 2 = iφ θ θ 2 θ verse plane with transverse axes x and y. In this e cos 2 sin 2 sin 2 Cartesian coordinate system, the unit Bloch vector ~a = (39) θ φ, θ φ, θ (sin cos sin sin cos ) is represented using the polar |α|2 αβ∗ ≤ θ ≤ π ≤ φ < π angle 0 and the azimuthal angle 0 2 , as = α∗β |β|2 (40) illustrated in Fig. 4 (a). Following our convention, state |0i at the north pole is associated with +1, and state |1i (the south pole) with −1. We can similarly represent the where I is the identity matrix, and ~σ = [σx, σy, σz] is a quantum state using the angles θ and φ, vector of Pauli matrices. If the Bloch vector ~a is a unit vector, then ρ represents a pure state ψ and Tr(ρ2) = 1. More generally, the Bloch sphere can be used to represent θ iφ θ |ψi = α|0i + β|1i = cos |0i + e sin |1i. (37) mixed states, for which |~a| < 1; in this case, the Bloch 2 2 vector terminates at points inside the unit sphere, and 0 ≤ Tr(ρ2) < 1. To summarize, the surface of the unit The Bloch vector is stationary on the Bloch sphere in sphere represents pure states, and its interior represents the rotating frame picture. If state |1i has a higher en- mixed states6. ergy than state |0i (as it generally does in superconduct- ing qubits), then in a stationary frame, the Bloch vector would precess around the z-axis at the qubit frequency (E1 −E0)/~. Without loss of generality (and much easier to visualize), we instead choose to view the Bloch sphere in a reference frame where the x and y-axes also rotate 2. Bloch-Redfield model of decoherence around the z-axis at the qubit frequency. In this rotating frame, the Bloch vector appears stationary as written in Within the standard Bloch-Redfield110–112 picture of Eq. (37). The rotating frame will be described in detail two-level system dynamics, noise sources weakly coupled in Section IV D 1 in the context of single-qubit gates. to the qubits have short correlation times with respect For completeness, we note that the density matrix ρ = to the system dynamics. In this case, the relaxation pro- 14 cesses are characterized by two rates (see Fig. 4): ground state (|0i) at the north pole, p = −1 is entirely in the excited state (|1i) at the south pole, and p = 0 is 1 longitudinal relaxation rate: Γ1 ≡ (41) a completely depolarized mixed state at the center of the T1 Bloch sphere. 1 Γ1 As illustrated in Fig. 4(b), longitudinal relaxation is transverse relaxation rate: Γ2 ≡ = + Γϕ T2 2 caused by transverse noise, via the x- or y-axis, with (42) the intuition that off-diagonal elements of an interaction Hamiltonian are needed to connect and drive transitions which contains the pure dephasing rate Γϕ. We note between states |0i and |1i. that the definition of Γ2 as a sum of rates presumes that Depolarization occurs due to energy exchange with an the individual decay functions are exponential, which oc- environment, generally leading to both an “up transition curs for Lorentzian noise spectra (centered at ω = 0) rate” Γ1↑ (excitation from |0i to |1i), and a “down tran- such as white noise (short correlation times) with a high- sition rate” Γ1↓ (relaxation from |1i to |0i). Together, frequency cutoff. these form the longitudinal relaxation rate Γ1: The impact of noise on the qubit can be visualized on the Bloch sphere in Fig. 4(a). For an initial state (t = 0) 1 Γ1 ≡ = Γ1↓ + Γ1↑. (45) T1 |ψi = α|0i + β|1i, (43) T1 is the 1/e decay time in the exponential decay func- the Bloch-Redfield density matrix ρBR for the qubit is tion in Eq. (44), and it is the characteristic time scale written113,114, over which qubit population will relax to its steady-
2 −Γ t ∗ iδωt −Γ t state value. For superconducting qubits, this steady- 1 + (|α| − 1)e 1 αβ e e 2 state value is generally the ground state, due to Boltz- ρBR = ∗ −iδωt −Γ t 2 −Γ t . (44) α βe e 2 |β| e 1 mann statistics and typical operating conditions. Boltz- mann equilibrium statistics lead to the “detailed bal- There are a few important distinctions between Eq. (44) ance” relationship Γ1↑ = exp(−~ωq/kBT )Γ1↓, where T and Eq. (40), which we list here and then describe in is the temperature and kB is the Boltzmann constant, more detail in subsequent sections. with an equilibrium qubit polarization approaching p = • First, we have introduced the longitudinal decay tanh(~ωq/2kBT ). Typical qubits are designed at fre- quency ωq/2π ≈ 5 GHz and are operated at dilution function exp(−Γ1t), which accounts for longitudi- nal relaxation of the qubit. refrigerator temperatures T ≈ 20 mK. In this limit, the up-rate Γ1↑ is exponentially suppressed by the Boltz- • Second, we introduced the transverse decay func- mann factor exp(−~ωq/kBT ), and so only the down- tion exp(−Γ2t), which accounts for transverse de- rate Γ1↓ contributes significantly, relaxing the popula- cay of the qubit. tion to the ground state. Thus, qubits generally spon- taneously lose energy to their cold environment, but the • Third, we have introduced an explicit phase ac- environment rarely introduces a qubit excitation. As a crual exp(iδωt), where δω = ωq − ωd, which gen- result, the equilibrium polarization approaches unity [see eralizes the Bloch sphere picture to account for Eq. (44)]118,119. cases where the qubit frequency ωq differs from the Only noise at the qubit frequency mediates qubit tran- rotating-frame frequency ωd, as we will see later sitions, whether absorption or emission, and this noise is when discussing measurements of T2 using Ramsey 115,116 generally “well behaved” (short correlation time, many interferometry , and in Section IV D 1 in the modes weakly coupled to qubit, no divergences) around context of single-qubit gates. the qubit frequency for superconducting qubits. The in- • And, fourth, we have constructed the matrix such tuition is that qubit-transition linewidths are relatively narrow in frequency, and so the noise generally does not that for t (T1,T2), the upper-left matrix ele- ment will approach unit value, indicating that all vary much over this narrow frequency range. Although there are a few notable exceptions, for example, qubit population relaxes to the ground state, while the 120–122 other three matrix elements decay to zero. This is decay in the presence of hot quasiparticles , which related to the assumption that the environmental can lead to non-exponential decay functions, longitudinal temperature is low enough that thermal excitations depolarization measurements generally present exponen- of the qubit from the ground to excited state rarely tial decay functions consistent with the Bloch-Redfield occur. picture. An example of a T1 measurement is shown in Fig. 5(a). Longitudinal relaxation The qubit is prepared in its excited state using an Xπ- The longitudinal relaxation rate Γ1 describes depolar- pulse, and then left to spontaneously decay to the ground ization along the qubit quantization axis, often referred state for a time τ, after which the qubit is measured. A to as “energy decay” or “energy relaxation.” In this lan- single measurement will project the quantum state into guage, a qubit with polarization p = 1 is entirely in the either state |0i or state |1i, with probabilities that cor- 15
(a) (c)
Xπ Readout Xπ Xπ/2 τ/2 τ/2 Xπ/2 Relaxation: τ Echo: t t
(b) (d)
Xπ/2 τ Xπ/2 Ramsey: t
117 FIG. 5. Characterizing longitudinal (T1) and transverse (T2) relaxation times of a transmon qubit . (a) Longitudinal relaxation (energy relaxation) measurement. The qubit is prepared in the excited state using an Xπ-pulse and measured after a waiting time τ. For each value τ, this procedure is repeated to obtain an ensemble average of the qubit polarization: +1 corresponding to |0i, and −1 corresponding to |1i. The resulting exponential decay function has a characteristic time T1 = 85 µs. (b) Transverse relaxation (decoherence) measurement via Ramsey interferometry. The qubit is prepared on the equator using an Xπ/2-pulse, intentionally detuned from the qubit frequency by δω, causing the Bloch vector to precess in the rotating frame at a rate δω around the z-axis. After a time τ, a second Xπ/2 pulse then projects the Bloch vector back on to the z axis, effectively mapping its former position on the equator to a position on the z axis. The oscillations decay with an ∗ approximately (but not exactly) exponential decay function, with a characteristic time T2 = 95 µs. (c) Transverse relaxation (decoherence) measurement via a Hahn echo experiment116. The qubit is prepared and measured in the same manner as the Ramsey interfometry experiment, except that a single Xπ pulse is applied midway through the free-evolution time τ. The decay function is approximately exponential, with a characteristic time T2E = 120 µs. The coherence improvement using the Hahn echo over panel (b) indicates that some low-frequency dephasing noise has been mitigated; however, a small amount remains since T2E has not yet reached the 2T1 limit. (d) Coherence function incorporating T1 loss and Gaussian dephasing components of the Ramsey interferometry data in panel (b). The Gaussian-distributed 1/f noise spectrum of magnetic flux 2 2 noise leads to a decay function exp(−t/2T1) exp(−χN ) = exp(−t/2T1) exp(−t /Tϕ,G) in Eq. (46). These two decay functions together match well the Ramsey data in panel (b). respond to the qubit polarization. To make an estimate in the x − y plane of the Bloch sphere. It is referred to of this polarization, one needs to identically prepare the as “pure dephasing,” to distinguish it from other phase- qubit and repeat the experiment many times. This is breaking processes such as energy excitation or relax- analogous to flipping a coin: any single flip will yield ation. heads or tails, but the probability of obtaining a heads or tails can be estimated by flipping the coin many times As illustrated in Fig. 4(c), pure dephasing is caused and taking the ensemble average. The resulting expo- by longitudinal noise that couples to the qubit via the z- nential decay has a characteristic time T1 = 85 µs. axis. Such longitudinal noise causes the qubit frequency ωq to fluctuate, such that it is no longer equal to the Pure dephasing rotating frame frequency ωd, and causes the Bloch vec- The pure dephasing rate Γφ describes depolarization tor to precess forward or backward in the rotating frame. 16
Intuitively, we can imagine identically preparing several Bloch vector back on to the z-axis. Repeated measure- instances of the Bloch vector along the x-axis. For each ments are made to take an ensemble averaged estimate of instance, the stochastic fluctuations of qubit frequency the qubit polarization, as a function of τ. The resulting will result in a different precession frequency, resulting oscillations in Fig. 5(b) feature an approximately expo- ∗ in a net fanout of the Bloch vector in the x − y plane. nential decay function with time T2 = 98 µs. The “*” This eventually leads to a complete depolarization of the indicates that the Ramsey experiment is sensitive to in- azimuthal angle φ. Note that this stochastic effect will be homogeneous broadening. That is, it is highly sensitive captured in the transverse relaxation rate Γ2 (next sec- to quasi-static, low-frequency fluctuations that are con- tion); it is not the deterministic term exp(±iδωt) that stant within one experimental trial, but vary from trial appears in Eq. (44), which represents intentional detun- to trial, e.g., due to 1/f-type noise. This sensitivity to ing of the qubit reference frame. quasi-static noise is related to the corresponding N = 0 There are a few important distinctions between pure noise filter function shown in Fig. 5(d) being centered dephasing and energy relaxation. First, in contrast to at at zero-frequency, as described in more detail in Sec- energy relaxation, pure dephasing is not a resonant phe- tion III D 2. nomenon; noise at any frequency can modify the qubit The Hahn echo shown in Fig. 5(c) is an experiment frequency and cause dephasing. Thus, qubit dephasing that is less sensitive to quasi-static noise. By placing is subject to broadband noise. Second, since pure de- a Yπ pulse at the center of a Ramsey interferometry ex- phasing is elastic (there is no energy exchange with the periment, the quasi-static contributions to dephasing can environment), it is in principle reversible. That is, the be “refocused,” leaving an estimate T2E that is less sen- dephasing can be “undone” – with quantum informa- sitive to inhomogeneous broadening mechanisms. The tion being preserved – through the application of unitary pulses are generally chosen to be resonant with the qubit operations, e.g., dynamical decoupling pulses78, see Sec. transition for a Hahn echo, since any frequency detuning III D 2. would be nominally refocused anyway. The resulting de- The degree to which the quantum information can be cay function in Fig. 5(c) is essentially exponential with retained depends on many factors, including the band- time T2E = 120 µs. width of the noise, the rate of dephasing, the rate at With the known T1 and T2 times, one can infer the pure which unitary operations can be performed, etc. This dephasing time Tϕ from Eq. (42), provided the decay should be contrasted with spontaneous energy relaxation, functions are exponential. In superconducting qubits, which is an irreversible process. Intuitively, once the however, the broadband dephasing noise (e.g., flux noise, qubit emits energy to the environment and its myriad charge noise, critical-current noise, ...) tends to exhibit uncontrollable modes, the quantum information is essen- a 1/f-like power spectrum. Such noise is singular near tially lost with no hope for its recovery and reconstitution ω = 0, has long correlation times, and generally does back into the qubit. not fall within the Bloch-Redfield description. The de- Transverse relaxation cay function of the off-diagonal terms in Eq. (44) are The transverse relaxation rate Γ2 = Γ1/2 + Γϕ describes generally non-exponential, and for such cases, the simple expression in Eq. (42) is not applicable. the√ loss of coherence of a superposition state, for example (1/ 2)(|0i+|1i), pointed along the x-axis on the equator of the Bloch sphere as illustrated in Fig. 4(d). Decoher- ence is caused in part by longitudinal noise, which fluctu- 3. Modification due to 1/f-type noise ates the qubit frequency and leads to pure dephasing Γϕ (red). It is also caused by transverse noise, which leads If we assume that the qubit is coupled to many in- to energy relaxation of the excited-state component of dependent fluctuators, then, regardless of their individ- the superposition state at a rate Γ1 (blue). Such a relax- ual statistics, they will in concert generate noise with ation event is also a phase-breaking process, because once a Gaussian distribution due to the central limit theo- it occurs, the Bloch vector points to the north pole, |0i, rem. We therefore say that the longitudinal fluctuations and there is no longer any knowledge of which direction exhibit Gaussian-distributed 1/f noise123,124. For 1/f the Bloch vector had been pointing along the equator; noise spectra, the phase decay function is itself a Gaus- the relative phase of the superposition state is lost. 2 sian exp −(t/Tϕ,G) , where we write Tϕ,G to distinguish Transverse relaxation T2 can be measured using Ram- it from Tϕ used in Eq. (42). Furthermore, this function sey interferometry, as shown and described in Fig. 5(b). is separable from the T1-type exponential decay, because The protocol positions the Bloch vector on the equator the T1-noise remains regular at the qubit frequency. The using a Xπ/2-pulse. Typically, the carrier frequency of density matrix in Eq. (44) becomes, following Refs. 78 this pulse is slightly detuned from the qubit frequency and 113, by an amount δω. As a result, the Bloch vector will pre- cess around the z-axis at a rate δω. This is done for 2 −Γ t ∗ iδωt − Γ1 t −χ (t)! 1 + (|α| − 1)e 1 αβ e e 2 e N convenience sake, so that the resulting Ramsey measure- ρ Γ , = ∗ −iδωt − 1 t −χ (t) 2 −Γ t ment will oscillate, making it easier to analyze. After α βe e 2 e N |β| e 1 precessing for a time τ, a second Xπ/2-pulse projects the (46) 17
For completeness, in addition to 1/f dephasing mech- anisms, we note that there are also “white” pure dephas- Qubit absorption Qubit emission ing mechanisms, which give rise to an exponential de- from environment to environment cay function for the dephasing component of T2. One common example is dephasing due to the shot noise of Thermal 1/f Johnson- residual photons in the readout resonator coupled to su- (Johnson) noise Nyquist perconducting qubits, as we discuss in Section III C 3. noise Quantum noise (Nyquist) noise 4. Noise power spectral density (PSD)
The frequency distribution of the noise power for a sta- tionary noise source λ is characterized by its PSD Sλ(ω) Z ∞ −iωτ “Low-frequency” Sλ(ω) = dτ hλ(τ)λ(0)ie . (47) noise −∞
FIG. 6. Examples of symmetric and asymmetric noise spec- The Wiener-Khintchine theorem states that the PSD tral densities. Noise at positive (negative) frequencies cor- is the Fourier transform of the autocorrelation function responds to the qubit emitting (absorbing) energy to (from) cλ(τ) = hλ(τ)λ(0)i of the noise source λ. Since the in- its environment. Thermal noise is proportional to temper- tegration limits are (−∞, ∞), this is the bilateral PSD. ature T and carries essentially a white noise spectrum. As Symmetrizing the PSD allows one to consider only posi- it represents a classical fluctuating parameter, such as elec- tive frequencies, which is termed a unilateral PSD. Both tric current, the noise power spectral density is symmetric unilateral and bilateral PSDs are used, often with the in frequency. When resonant with the qubit, it will drive same notation, and so one needs to know how the PSD both stimulated emission and absorption processes. The qubit is defined, keep track of the factors of 2 and π, and also may also spontaneously emit energy to its environment, rep- be aware of the implications for quantum systems. resented as Nyquist noise125, a quantum mechanical effect that is not symmetric in frequency. At sufficiently low tem- For classical systems, the noise power spectral density is symmetric. This is because the autocorrelation func- peratures or high frequencies, ~ω > 2kBT , the Nyquist noise dominates thermal noise. Another common example is 1/f tion of real signals is itself a real function, and the Fourier noise, which is also a classical noise fluctuation and symmet- transform of a real temporal function is symmetric in the ric in frequency. frequency domain. Dephasing noise is caused by real, fluctuating fields, and so its PSD is generally symmetric. Examples of such classical noise include thermal (John- son) noise and 1/f noise126 (see Fig. 6). where the decay function hexp(−χN (t))i contains the co- In turn, the inverse Fourier transform of the PSD will herence function χN (t), which generalizes pure dephas- yield the autocorrelation function: ing to include non-exponential decay functions. As we shall see later, the subscript N labeling the decay func- Z ∞ 1 iωτ π cλ(τ) = dω Sλ(ω)e . (48) tion refers to the number of -pulses used to refocus the 2π low-frequency noise, which impacts the form of the de- −∞ cay function. Because the function is no longer purely This implies that integrating the noise power spectral exponential, we cannot formally write the transverse re- density with τ = 0 yields the second moment of the noise, laxation decay function as exp(−t/T2). However, an ex- or, for zero-mean fluctuations, the variance. ponential decay remains a practically reasonable approx- However, the autocorrelation function for a quantum imation for Tϕ & T1. We also note that the energy decay system may be complex-valued due to the fact that quan- component of the transverse relaxation is exp(−t/2T1), tum operators generally do not commute at different and so T2 can never be larger than 2T1. In the absence times. This means that time-ordering of the operators of pure dephasing, the maximum T2 = 2T1 is reached. matters, and the PSD need not be symmetric in fre- As an example, consider the Ramsey interferometry quency. This is generally the case for transverse noise data in Fig. 5(b). Since the dephasing is relatively weak, causing longitudinal energy relaxation. Noise at a posi- the transverse relaxation function as exp(−t/T2) is a tive frequency S(ωq) corresponds to energy transfer from reasonable fit and yields T2 = 95 µs. However, us- the qubit to the environment, including both stimulated ing the value T1 = 85 µs from Fig. 5(a) and dividing and spontaneous emission, associated with the down-rate out exp(−t/2T1) from the data in Fig. 5(b), the remain- Γ1↓. Noise at a negative frequency S(−ωq) corresponds ing pure dephasing decay function is shown in Fig. 5(d) to energy transfer to the qubit from the environment, as- and assumes a Gaussian envelope hexp(−χN (t))i = sociated with the up-rate Γ1↑. For a detailed discussion, 2 exp −(t/Tϕ,Gt) , with Tϕ,G = 98 µs. The Hahn echo see Refs. 127 and 128. Spontaneous emission to a cold data in Fig. 5(c) may be treated similarly. environment or electromagnetic vacuum, represented by 18
2 2 Nyquist noise in Fig. 6, is an example of an asymmetric SQ(ω) = BQ[ω/(2π ×1Hz)], where the noise strength BQ noise PSD125. at 1 Hz can assume a range of values depending on the In general, making a connection between Sλ(ω) and level of dissipation in the system. Likewise, the cross-over the measured qubit decay functions is the basis for noise from 1/f-like behavior to f-like behavior generally occurs spectroscopy up to second-order statistics78,129–132. The at around 1 GHz, but will vary higher or lower between search for higher-order spectra related to non-Gaussian samples depending on the degree of dissipation62,136. noise is a current topic of active research133.
C. Common examples of noise 2. Magnetic flux noise There are many sources of stochastic noise in super- conducting qubits, and we refer the reader to Ref. 40 for Another commonly observed noise in solid-state de- a review. Here, we briefly present several of the most vices is magnetic flux noise. The origin of this noise is common types of noise, their affect on coherence, and understood to arise from the stochastic flipping of spins refer the reader to the references for a more detailed dis- (magnetic dipoles) that reside on the surfaces of the su- cussion. perconducting metals comprising the qubit137, resulting in random fluctuations of the effective magnetic field that biases flux-tunable qubits. 1. Charge noise For example, in the case of the split transmon, the external magnetic field threading the loop couples longi- Charge noise is ubiquitous in solid-state devices. It tudinally to the qubit and modulates the transition fre- arises from charged fluctuators present in the defects or quency via the Josephson energy EJ (except at ϕe = 0, charge traps that reside in interfacial dielectrics, the junc- where the qubit is first-order insensitive to magnetic- tion tunnel barrier, and in the substrate itself. These field fluctuations). Because the flux noise is longitudinal are often modeled as an ensemble of fluctuating two-level to the transmon, it contributes to pure dephasing (Tϕ). 134,135 systems or as bulk dielectric loss . For example, in However, in the case of the flux qubit, and depending on the case of a transmon qubit, the electric field between the flux-bias point, the flux noise may be either longitu- the capacitor plates traverses and couples to dielectric dinal – causing dephasing Tϕ – or it may couple trans- 62,78 defects residing on the metal surfaces of the plates (for versely and thus contribute to T1 relaxation . The lateral-plate-type capacitors) or the capacitor dielectric noise power spectrum of these fluctuations generally ex- between the plates (for parallel-plate-type capacitors). hibits a “quasi-universal” dependence, The electric field variable is transverse with respect to the quantization axis of the transmon qubit, which means γΦ that this noise is mainly responsible for energy relaxation 2 2π × 1Hz SΦ(ω) = AΦ , (50) (T1). Additionally, if the EJ /EC ratio of the transmon is ω not made sufficiently large (smaller than around 60), the qubit frequency itself will also be sensitive to broadband 2 2 with γΦ ≈ 0.8 − 1.0 and AΦ ≈ (1 µΦ0) /Hz, and has charge fluctuations. In this case, low-frequency charge been shown to extend from less than millihertz to beyond noise couples longitudinally to the transmon and causes gigahertz frequencies78,131,132,138,139. pure dephasing (Tϕ). Charge noise is modeled primarily as a combination of The large, low-frequency weighting of the 1/f power inverse-frequency noise and Nyquist noise, also referred distribution enables the use of engineered error miti- ohmic gation techniques – such as dynamical decoupling – to to as noise. At lower frequencies, the spectral 78,140–142 density takes the form achieve better coherence and for improving sin- gle and two-qubit gate fidelity143. It was recently demon- γQ /f T 2 2π × 1Hz strated that 1 flux noise is also a 1-mechanism when SQ(ω) = A , (49) 62 Q ω extended out to the qubit frequency , and one similarly expects a crossover to ohmic flux noise at high enough 2 −3 2 frequencies144. with quasi-universal values AQ = (10 e) /Hz at 1 Hz, and γQ ≈ 1. In addition to large 1/f fluctuations, early Although much is known about the statistics and charge qubits often witnessed discrete, charge offsets rem- number of the defects presumed responsible for flux iniscent of random telegraph noise. Together, these two noise, their precise physical manifestation remains uncer- mechanisms severely limited the utility of charge qubits, tain137,145. The fact that the 1/f noise is quasi-universal and served as a strong motivation to move to capaci- and largely independent of device, strongly suggests a tively shunted charge qubits (transmons), which greatly common origin for the noise. Recent studies suggest that reduced the qubit longitudinal sensitivity to charge noise. adsorbed molecular oxygen may be responsible for flux- At higher frequencies, the power spectrum takes the form noise145,146. 19
3. Photon number fluctuations pair149,150. Although this quasiparticle generation mech- anism is not yet well understood, it has been shown that In the circuit QED architecture, resonator pho- quasiparticles can be transiently pumped away, improv- 122 ton number fluctuation is another major decoherence ing T1 times and reducing T1 temporal variation . source147. Residual microwave fields in the cavity have photon-number fluctuations that in the dispersive regime D. Operator form of qubit-environment interaction impact the qubit through an interaction term χσzn, see Sec. II C 2, leading to a frequency shift ∆Stark = 2ηχn¯, wheren ¯ is the average photon number, and η = κ2/(κ2 + Similar to the way that two qubits are coupled, a qubit 4χ2) effectively scales the photon population seen by the may couple and interact with uncontrolled degrees of qubit due to the interplay between the qubit-induced dis- freedom (DOF) in its environment (the noise sources). persive shift of the resonator frequency (χ) and the res- The interaction Hamiltonian between the qubit DOF onator decay rate (κ). (Oˆq) and those of the noise source (λˆ) may be expressed In the dispersive limit, the noise is longitudinally cou- in a general form pled to the qubit and leads to pure dephasing at a rate, 2 4χ Hˆint = νOˆqλˆ (53) Γφ = η n.¯ (51) κ where ν denotes the coupling strength – which is related The fluctuations originate from residual photons in the to the sensitivity of the qubit to environmental fluctua- resonator, typically due to radiation from higher tem- tions ∂Hˆq/∂λ – and we assume that Oˆq is a qubit opera- perature stages in the dilution refrigerator107,148. The tor within the qubit Hamiltonian Hˆq. The noisy environ- corresponding noise spectral density is of a Lorentzian λˆ type, ment represented by the operator produces fluctuations δλ. Note that we retained the hats in this section to re- 2 2ηnκ¯ mind us that these are quantum operators. S(ω) = 4χ , (52) ω2 + κ2 which exhibits an essentially white noise spectrum up to 1. Connecting T to S(ω) a 3dB cutoff frequency ω = κ set by the resonator decay 1 rate κ, see Ref. 62. If the coupling is transverse to the qubit, e.g. Oˆq † is of the type σx or (a + a ) – see the related case of 4. Quasiparticles qubit-qubit coupling treated in Sec. II C – then noise at the qubit frequency can cause transitions between Quasiparticles, i.e. unpaired electrons, are another im- the qubit eigenstates. Since this is a stochastic process, portant noise source for superconducting devices121. The the ensemble-average manifests itself as a decay (usually tunneling of quasiparticles through a qubit junction may exponential) of the qubit population towards a certain equilibrium value (usually the qubit ground state |0i for lead to both T1 relaxation and pure dephasing Tϕ, de- pending on the type of qubit, the bias point, and the kBT ~ωq). Again, this process is equivalently referred junction through which the tunneling event occurs120,122. to as “T1 relaxation”, “energy relaxation”, or “longitudi- Quasiparticles are naturally excited due to thermo- nal relaxation”. As stated above, T1 is the characteristic dynamics, and the quasiparticle density in equilibrium time scale of the decay. Its inverse, Γ1 = 1/T1 is called superconductors should be exponentially suppressed as the relaxation rate and depends on the power spectral temperature decreases. However, below about 150 mK, density of the noise S(ω) at the transition frequency of the quasiparticle density observed in superconducting de- the qubit ω = ωq: vices – generally in the range 10−8−10−6 per Cooper pair – is much higher than BCS theory would predict for a su- 2 1 ∂Hˆq perconductor in equilibrium with its cryogenic environ- h | | i S ω , Γ1 = 2 0 1 λ( q) (54) ment at 10 mK. The reason for this excess quasiparticle ~ ∂λ population is unclear, but it is very likely related to the presence of additional, non-thermal mechanisms that in- where ∂Hˆq/∂λ is the qubit transverse susceptibility to crease the generation rates, “bottleneck effects” that oc- fluctuations δλ, such that |δλ|2 is the ensemble average cur at millikelvin temperatures to reduce recombination value of the environmental noise sources as seen by the rates, or a combination of both. qubit. Eq. (54) is equivalent to Fermi’s Golden Rule, It has been shown that the observed T1 and excess in which the qubit’s transverse susceptibility to noise is excited-state population measured in today’s state-of- driven by the noise power spectral density. The qubit the-art high-coherence transmon are self-consistent with transverse susceptibility can be used to calculate the pref- excess “hot” nonequilibrium quasiparticles at the quasi- actors; for example, for fluctuations δλ = δn, the rele- universal density of around 10−7 − 10−6 per Cooper vant term in the transmon Hamiltonian in Eq. (16) is 20
(a) CP / CPMG cally modulate the transition frequency of the qubit and 2 2 thereby introduce a stochastic phase evolution of a qubit 2N N N N N N 2N superposition state. This gradually leads to a loss of phase information, and it is therefore called pure dephas- N -pulses ( N > 1) t ing (time constant Tϕ). Unlike T1 relaxation, which is 1 generally an irreversible (incoherent) error, pure dephas- s ing Tϕ is in principle reversible (a coherent error). The N = 0 CPMG degree of pure dephasing depends on the control pulse sequence applied while the qubit is subject to the noise 1 process. 2 ( ) 6 10 Consider the relative phase ϕ of a superposition state N
g undergoing free evolution in the presence of noise. The S~1/f superposition state’s accumulated phase, Z t 0 ϕ(t) = ωqdt = hωqit + δϕ(t) (55) 0 0 2 4 6 8 0 Frequency, f (MHz) (b) diffuses due to adiabatic fluctuations of the transition frequency, 2T1 Z t ∂ωq 10 δϕ(t) = δλ(t0)dt0, (56) CPMG simulation ∂λ 0
s) where ∂ωq/∂λ = (1/~)|h∂Hˆq/∂λi| is the qubit’s longi- CPMG tudinal sensitivity to λ-noise. For noise generated by a ( 2 T 1 large number of fluctuators that are weakly coupled to the qubit, its statistics are Gaussian. Ensemble averaging over all realizations of the Gaussian-distributed stochas- Ramsey ( N=0) tic process δλ(t), the dephasing is
1 2 i δϕ(t) − hδϕ (t)i −χN (t) 0.1 he i = e 2 ≡ e , (57) 1 10 100 1000 Number of pulses, N leading to a coherence decay function,
FIG. 7. Dynamical error suppression. (a) Carr-Purcell- τ 2 ∂ω Z ∞ Meiboom-Gill (CPMG) pulse sequence applies N equally −χN (τ) q he i = exp − gN (ω, τ)S(ω)dω , spaced π pulses within an otherwise free-evolution time τ. 2 ∂λ −∞ Pulses in the time domain correspond to bandpass filters in (58) the frequency domain (lower panel) which serve to shape the where g(ω, τ) is a dimensionless weighting function. noise power spectrum seen by the qubit. The centroid of the The function gN (ω, τ) can be viewed as a frequency- bandpass filter shifts to higher frequencies as N is increased. domain filter of the noise Sλ(ω) [see Fig. 7(a)]. In gen- For noise that decreases with frequency, such as 1/f noise, eral, its filter properties depend on the number N and larger N corresponds to less integrated noise impinging on distribution of applied pulses. For example, considering the qubit. (b) CPMG pulse sequence applied to a flux qubit 78,154–158 biased at a point that is highly sensitive to 1/f flux noise. The sequences of π-pulses , Ramsey (N = 0) time is approximately 300 ns, and the Hahn echo (N = 1) time is approximately 1.5 µs. Increasing the 1 1+N gN (ω, τ) = 1 + (−1) exp(iωτ)+ number of CPMG pulses continues to increase the effective (ωτ)2 T time towards the 2T limit. Adapted from Ref. 78. 2 1 N 2 X j 2 (−1) exp(iωδjτ) cos(ωτπ/2) , (59) j=1 2 4EC (ˆn − ng) , where we allow for an offset charge ng, and the susceptibility is given by 8EC nˆ. We refer the where δj ∈ [0, 1] is the normalized position of the centre reader to Refs. 151–153 for more details. of the jth π-pulse between the two π/2-pulses, τ is the total free-induction time, and τπ is the length of each 157,158 π-pulse , yielding a total sequence length τ + Nτπ. 2. Connecting Tϕ to S(ω) As the number of pulses increases for fixed τ, the filter function’s peak shifts to higher frequencies, leading to a If the coupling to the qubit is instead longitudinal, reduction in the net integrated noise for 1/f α-type noise † e.g. Hˆq is of the type σz or a a, the noise will stochasti- spectra with α > 0. Similarly, for a fixed N, the filter 21 function will shift in frequency with τ. Additionally, for a during free evolution – periods of time for which no con- fixed time separation τ 0 = τ/N (valid for N ≥ 1), the fil- trol is applied to the qubit, except for very short dynami- ter sharpens and asymptotically peaks at ω0/2π = 1/2τ 0 cal decoupling pulses – and during driven evolution – pe- as more pulses are added. gN (ω, τ) is thus called the “fil- riods of time during which the control fields are applied ter function”78,156, and it depends on the pulse sequences to the qubit. Both free-evolution and driven-evolution being applied. From Eq. (58), the pure dephasing de- noise is important to characterize, as the noise PSD may cay arises from a noise spectral density that is “shaped” differ for these two types of evolution, and both are uti- or “filtered” by the sequence-specific filter function. By lized in the context of universal quantum computation. choosing the number of pulses, their rotation axes, and We refer the reader to Ref. 132 for a summary of noise their arrangement in time, we can design filter functions spectroscopy during both types of evolution. that minimize the net noise power for a given noise spec- The Ramsey frequency itself is sensitive to longitudinal tral density within the experimental constraints of the ex- noise, and monitoring its fluctuations is one means to periment (e.g., pulse-modulation bandwidth of the elec- map out the noise spectral density over the sub-millihertz tronics used to control the qubits). to ∼ 100 Hz range131,163. To give a standard example, we compare the coherence integral for two cases: a Ramsey pulse sequence and a At higher frequencies, the CPMG dynamical decou- Hahn echo pulse sequence. Both sequences involve two pling sequence can be used to create narrow-band fil- π/2 pulses separated by a time τ, during which free evo- ters that “sample” the noise at different frequencies as a τ lution of the qubit occurs in the presence of low-frequency function of the free-evolution time and the number of pulses N. This has been used to map out the noise PSD dephasing noise. The distinction is that the Hahn echo 78 will place a single π pulse (N = 1) in the middle of the in the range 0.1 - 300 MHz . One must be careful of free-evolution period, whereas the Ramsey does not use the additional small peaks at higher-frequencies, which all contribute to the dephasing used to perform the noise any additional pulses (N = 0). The resulting filter func- 164 tions are: spectroscopy . In fact, using pulse envelopes such as Slepians165 – which are designed to have concentrated frequency re- 2 ωτ g0(ω, τ) = sinc (60) sponse – to perform noise spectroscopy is one means to 2 reduce such errors157. 2 ωτ 2 ωτ g1(ω, τ) = sin sinc (61) 4 4 At even higher frequencies, measurements of T1 can be used in conjunction with Fermi’s golden rule to map out where the subscript N = 0 and N = 1 indicate the num- the transverse noise spectrum above 1 GHz62,78,166. ber of π-pulses applied for the Ramsey and Hahn echo The aforementioned are all examples of noise spec- experiments, respectively. The filter function g0(ω, τ) for troscopy during free evolution. Noise spectroscopy dur- the Ramsey case is a sinc-function centered at ω = 0. For ing driven evolution was also demonstrated using a “spin- noise that decreases with frequency, e.g., 1/f flux noise locking” technique, where a strong drive along x or y in superconducting qubits, the Ramsey experiment win- axes defines a new qubit quantization axis, whose Rabi dows through the noise in S(ω) where it has its highest frequency is the new qubit frequency in the spin-locking value. This is the worst choice of filter function for 1/f frame. The spin-locking frame is then used to infer the noise. In contrast, the Hahn echo filter function has a noise spectrum while the qubit is continually subject to centroid that is peaked at a higher frequency, away from a driving field. For more information, we refer the reader ω = 0. In fact, it has zero value at ω = 0. For noise to Ref. 132. that decreases with frequency, such as 1/f noise, this is advantageous. This concept extends to larger numbers N of π pulses, and is called a Carr-Purcell-Meiboom-Gill 159,160 (CPMG) sequence . In Fig. 7(b), the T2 time of a qubit under the influence of strong dephasing noise is in- E. Engineering noise mitigation creased toward the 2T1 limit using a CPMG dynamical error-suppression pulse sequence with an increasing num- Here, we briefly review a few examples of techniques N ber of pulses, . We refer the reader to Refs. 78, 161, and that have been developed to reduce noise or reduce its im- 162, where these experiments were performed with super- pact on decoherence (sensitivity). We stress that improv- conducting qubits. ing gate fidelity is a comprehensive optimization task, one that is full of trade-offs. It is thus important to identify what the limiting factors are, what price we have to pay 3. Noise spectroscopy to diminish these limiting factors, and what advantage we can achieve until reaching a better trade-off. These all The qubit is highly sensitive to its noisy environment, require an accurate understanding the limitations on the and this feature can be used to map out the noise power gate fidelity, the sources of decoherence, the properties spectral density. In general, one can map the noise PSD of the noise, and how it affects the system performance. 22
1. Materials and fabrication improvements Returning to excess quasiparticles, it has been shown that quasiparticles can be stochastically pumped away from the qubit region, resulting in longer, and more sta- Numerous efforts have been undertaken to reduce 122 noise-induced defects due to materials and fabrica- ble T1 times . Although the pumping technique uses tion40,167. In the case of charge noise, significant efforts a series of π-pulses, this technique differs from dynami- have been made to reduce the number of defects, such cal error suppression of coherent errors in that pulses are as substrate cleaning59,168, substrate annealing169, and stochastically applied, and that it addresses incoherent trenching41,61. In the case of flux noise, several groups errors (T1). have performed experiments to characterize the behavior and properties of magnetic-flux defects137,170,171. More recently, a number of groups have tried optical surface 4. Cryogenic engineering treatments to remove these defects145. In the context of residual quasiparticles, it has been In the case of photon shot-noise, in addition to ap- shown that adding quasiparticle traps to the circuit de- plying dynamical decoupling techniques, there have been sign can reduce the quasiparticle number, particularly in several recent works aimed at reducing the thermal pho- devices that create excess quasiparticles, such as classi- ton flux that reaches the device. This include optimizing cal digital logic or operation in the presence of thermal 107,149,173 172 the attenuation of the cryogenic setup , remak- radiation ing the cryogenic attenuators with more efficient heat sinking148, adding absorptive “black” material to absorb stray thermal photons174,175, and adding additional cav- 2. Design improvements ity filters for thermalization176.
Another strategy is to reduce qubit sensitivity to the noise by design. A qubit can only lose energy to defects IV. QUBIT CONTROL if it couples to them. It has been demonstrated that altering the capacitor geometry to increase the electric- In this section, we will introduce how superconduct- field mode volume reduces the electric field density in the ing qubits are manipulated to implement quantum al- thin dielectric regions that cause loss. This effectively gorithms. Since the transmon-like variety of supercon- reduces the “participation” of the defects and makes the 55,62,134 ducting qubits has so far been the most widely deployed qubits less senstivie to these noise sources. . modality for implementing quantum programs, the dis- In another example, the split transmons built using cussion throughout this section will be focused on modern asymmetric junctions have lower sensitivity to flux noise techniques for transmons. Nonetheless, the techniques than their symmetric counterparts at the expense of de- 69 introduced here are applicable to all types of supercon- creased frequency tunability . This is a good trade-off ducting qubits. to make, because generally one is interested in tuning the We start with a brief review of the gates used in clas- qubit frequency over a somewhat restricted range (typ- sical computing as well as quantum computing, and the ically around 1 GHz) about the qubit frequency. When concept of universality. Subsequently we discuss the most such asymmetric transmons are used in a gate scheme common technique of driving single qubit gates via a CPHASE 65 such as the adiabatic -gate , (see Sec.IV F) the capacitive coupling of a microwave line, coupled to the qubit is less sensitive to flux noise, has a lower dephasing qubit. We introduce the notion of “virtual Z gates” and rate, and this should improve the gate fidelity in general. “DRAG” pulsing. In the latter part of this section, we re- view some of the most common implementations of two- qubit gates in both tunable and fixed-frequeny transmon 3. Dynamical error suppression qubits. The single-qubit and two-qubit operations to- gether form the basis of many of the medium-scale su- As introduced in the previous section, it is advan- perconducting quantum processors that exist today. tageous to leverage the 1/ω distribution of flux noise, Throughout this section, we write everything in the wherein a considerable amount of the noise power re- computational basis {|0i, |1i} where |0i is the +1 eigen- sides at low frequencies, and so the noise is “quasi-static”. state of σz and |1i is the −1 eigenstate. We use capi- The spin-echo technique116, which disrupts the free evo- talized serif-fonts to indicate the rotation operator of a lution by a π-pulse, is extremely effective in mitigating qubit state, e.g. rotations around the x-axis by an angle the pure dephasing by refocusing the coherent phase dis- θ is written as persion due to low-frequency noise. The more advanced θ −i σx versions, such as the CPMG-sequence, use multiple π- Xθ = RX (θ) = e 2 = cos(θ/2)1 − i sin(θ/2)σx (62) pulses to interrupt the system more frequently, pushing the filter band to even higher frequencies – a technique and we use the shorthand notation ‘X’ for a full π rotation 78 known as dynamical decoupling . about the x axis (and similarly for Y := Yπ and Z := Zπ). 23
CIRCUIT A. Boolean logic gates used in classical computers GATE SYMBOL TRUTH TABLE
Universal boolean logic can be implemented on classi- NOT The output is 1 Input Output cal computers using a small set of single-bit and two-bit when the input is 0 and 0 when 0 1 gates. Several common classical logic gates are shown in the input is 1. 1 0 Fig. 8 along with their truth tables. In classical boolean logic, bits can take on one of two values: state 0 or state 1. The state 0 represents logical FALSE, and state 1 rep- resents logical TRUE. Input Output Beyond the trivial “identity operation,” which simply AND The output is 1 0 0 0 only when both passes a boolean bit unchanged, the only other possible inputs are 1, 0 1 0 single-bit boolean logic gate is the NOT gate. As shown otherwise the 1 0 0 in Fig. 8, the NOT gate flips the bit: 0 → 1 and 1 → 0. output is 0. 1 1 1 This gate is reversible, because it is trivial to determine the input bit value given the output bit values. As we Input Output will see, for two-bit gates, this is not the case. OR The output is 0 0 0 0 only when both There are several two-bit gates shown in Fig. 8. A inputs are 0, 0 1 1 two-bit gate takes two bits as inputs, and it passes as otherwise the 1 0 1 an output the result of a boolean operation. One com- output is 1. 1 1 1 mon example is the AND gate, for which the output is 1 if and only if both inputs are 1; otherwise, the output Input Output is 0. The AND gate, and the other two-bit gates shown NAND The output is 0 0 0 1 in Fig. 8, are all examples of irreversible gates; that is, only when both 0 1 1 the input bit values cannot be inferred from the output inputs are 1, values. For example, for the AND gate, an output of log- otherwise the 1 0 1 output is 1. 1 1 0 ical 1 uniquely identifies the input 11, but an output of 0 could be associated with 00, 01, or 10. Once the opera- tion is performed, in general, it cannot be “undone” and Input Output the input information is lost. There are several variants NOR The output is 1 0 0 1 only when both of two-bit gates, including, 0 1 0 inputs are 0, otherwise the 1 0 0 • AND and OR; output is 0. 1 1 0 • NAND (a combination of NOT and AND) and NOR Input Output XOR The output is 1 (a combination of NOT and OR); only when the two 0 0 0 inputs have 0 1 1 • XOR (exclusive OR) and NXOR (NOT XOR). di!erent value, 1 0 1 otherwise the 1 1 0 output is 0. The XOR gate is interesting, because it is a parity gate. That is, it returns a logical 0 if the two inputs are the Input Output XNOR The output is 1 same values (i.e., they have the same parity), and it re- only when the two 0 0 1 turns a logical 1 if the two inputs have different values inputs have the 0 1 0 (i.e., different parity). Still, the XOR and NXOR gates same value, 1 0 0 otherwise the are not reversible, because knowledge of the output does output is 0. 1 1 1 not allow one to uniquely identify the input bit values. The concept of universality refers to the ability to per- FIG. 8. Classical single-bit and two-bit boolean logic gates. form any boolean logic algorithm using a small set of For each gate, the name, a short description, circuit repre- single-bit and two-bit gates. A universal gate set can in sentation, and input/output truth tables are presented. The principle transform any state to any other state in the numerical values in the truth table correspond to the classical state space represented by the classical bits. The set of bit values 0 and 1. Adapted from Ref. 177. gates which enable universal computation is not unique, and may be represented by a small set of gates. For ex- ample, the NOT gate and the AND gate together form a universal gate set. Similarly, the NAND gate itself is universal, as is the NOR gate. The efficiency with which one can implement arbitrary boolean logic, of course, de- pends on the choice of the gate set. 24
CIRCUIT MATRIX TRUTH BLOCH GATE REPRESENTATION REPRESENTATION TABLE SPHERE
z I Identity-gate: 1 0 Input Output I = no rotation is I 0 1 performed. ( ) |0〉 |0〉 |1〉 |1〉 y x
z X gate: 0 1 Input Output rotates the X = qubit state by X ( 1 0 ) |0〉 |1〉 π radians |1〉 |0〉 (180º) about the x-axis. 180º y x
z Y gate: 0 i Input Output rotates the Y = qubit state by Y ( i 0 ) |0〉 i |1〉 π radians 180º (180º) about |1〉 !i |0〉 the y-axis. y x
z Z gate: 180º 1 0 Input Output rotates the Z = qubit state by Z ( 0 1) |0〉 |0〉 π radians |1〉 !|1〉 (180º) about the z-axis. y x
90º z S gate: Input Output rotates the 1 0 π S = π S (0 e i 2) qubit state by 2 |0〉 |0〉 radians (90º) π |1〉 i 2 |1〉 about the e z-axis. y x
45º z T gate: Input Output rotates the 1 0 π T = π T (0 e i 4) qubit state by 4 |0〉 |0〉 radians (45º) π |1〉 e i 4 |1〉 about the z-axis. y x
H gate: Input Output rotates the 1 1 1 z H = |0〉 |0〉 + |1〉 qubit state by H 2 ( 1 1 ) π radians (180º) 2 180º about an axis diagonal in the |1〉 |0 〉 |1 〉 x-z plane. This is 2 y equivalent to an X-gate followed x π by a 2 rotation about the y-axis.
FIG. 9. Quantum single-qubit gates. For each gate, the name, a short description, circuit representation, matrix representation, input/output truth tables, and Bloch sphere represenation are presented. Matrices are defined in the basis spanned by the state vectors |0i ≡ [1 0]T and |1i ≡ [0 1]T . The numerical values in the truth table correspond to the quantum states |0i and |1i. Adapted from Ref. 177.
B. Quantum logic gates used in quantum computers course assume the classical states |0i and |1i, at the north
Quantum logic can similarly be performed by a small set of single-qubit and two-qubit gates. Qubits can of 25
CIRCUIT MATRIX TRUTH GATE REPRESENTATION REPRESENTATION TABLE
Input Output Controlled-NOT gate: 1 0 0 0 |00 〉 |00 〉 apply an X-gate to the 0 1 0 0 CNOT = |01 |01 target qubit if the 0 0 0 1 〉 〉 control qubit is in state 0 0 1 0 |10 〉 |11 〉 |1 〉 |11 〉 |10 〉
Input Output Controlled-phase gate: 1 0 0 0 |00 〉 |00 〉 apply a Z-gate to the 0 1 0 0 target qubit if the CPHASE = |01 〉 |01 〉 0 0 1 0 control qubit is in |10 〉 |10 〉 state |1 Z 0 0 0 1 〉 |11 〉 |11 〉
FIG. 10. Quantum two-qubit gates: the controlled NOT (CNOT) gate and the controlled phase (CPHASE or CZ). For each gate, the name, a short description, circuit representation, matrix representation, and input/output truth tables are presented. Matrices are defined in the basis spanned by the two-qubit state vectors |00i ≡ [1 0 0 0]T , |01i ≡ [0 1 0 0]T , |10i ≡ [0 0 1 0]T , and |11i ≡ [0 0 0 1]T , where the first qubit is the control qubit, and the second qubit is the target qubit. The CNOT gate flips the state of the target qubit conditioned on the control qubit being in state |1i. The CPHASE gate applies a Z gate to the target qubit conditioned on the control qubit being in state |1i. Adapted from Ref. 177. pole and south pole of the Bloch sphere, but they can also can be built from the CNOT-gate and single-qubit gates. assume arbitrary superpositions α|0i + β|1i, correspond- ing to any other position on the sphere. Single-qubit operations translate an arbitrary quan- tum state from one point on the Bloch sphere to another point by rotating the Bloch vector (spin) a certain an- gle about a particular axis. As shown in Fig. 9, there are several single-qubit operations, each represented by a matrix that describes the quantum operation in the computational basis represented by the eigenvectors of T T the σz operator, i.e. |0i ≡ [1 0] and |1i ≡ [0 1] . For example, the identity gate performs no rotation on the state of the qubit. This is represented by a two-by- two identity matrix. The X-gate performs a π rotation about the x axis. Similarly, the Y-gate and Z-gate per- form a π rotation about the y axis and z axis, respec- tively. The S-gate performs a π/2 rotation about the z axis, and the T-gate performs a rotation of π/4 about the z axis. The Hadamard gate H is also a common single-qubit gate the performs a π rotation about an axis diagonal in the x-z plane, see Fig. 9. Two-qubit quantum-logic gates are generally condi- tional gates and take two qubits as inputs. Typically, the first qubit is the control qubit, and the second is the target qubit. A unitary operator is applied to the target qubit, dependent on the state of the control qubit. The two common examples shown in Fig. 10 are the controlled NOT (CNOT-gate) and controlled phase (CZ or CPHASE gate). The CNOT-gate flips the state of the target qubit conditioned on the control qubit being in state |1i. The CPHASE-gate applies a Z gate to the target qubit, condi- tioned on the control qubit being in state |1i. As we will shown later, the iSWAP gate – another two-qubit gate – 26
CNOT (a) The unitary operator of the gate can be written Classical NOT gate (b) Quantum X gate in a useful way, highlighting that it applies an X depend- ing on the state of the control qubit. b b
1 0 0 0 in out 0 1 0 0 UCNOT = = |0ih0| ⊗ 1 + |1ih1| ⊗ X (63) 0 0 0 1 b b 0 0 1 0 and similarly for the CPHASE gate, FIG. 11. Comparison of the classical inverter (NOT) gate and quantum bit flip (X) gate. (a) The classical NOT gate 1 0 0 0 that inverts the state of a classical bit. (b) The quantum X 0 1 0 0 gate, which flips the amplitudes of the two components of a UCPHASE = = |0ih0| ⊗ 1 + |1ih1| ⊗ Z (64) 0 0 1 0 quantum bit. 0 0 0 −1
Comparing the last equality above with the unitary for A universal set of single-qubit and two-qubit gates is the CNOT [Eq. (63)], it is clear that the two gates are sufficient to implement arbitrary quantum logic. This closely related. Indeed, a CNOT can be generated from means that this gate set can in principle reach any state a CPHASE by applying two Hadamard gates, in the multi-qubit state-space. How efficiently this is done depends on the choice of quantum gates that com- 1 1 UCNOT = ( ⊗ H)UCPHASE( ⊗ H), (65) prise the gate set. We also note that each of the single- qubit and two-qubit gates is reversible, that is, given the HZH X CPHASE since = . Due to the form of Eq. (64), the output state, one can uniquely determine the input state. CZ gate is also denoted the gate, since it applies a con- As we discuss further, this distinction between classical Z CNOT trolled operator, by analogy with (a controlled and quantum gates arises, because quantum gates are X application of operator). Inspection of the definition of based on unitary operations U. If a unitary operation U CPHASE in Fig. 10 makes no distinction between which is a particular gate applied to a qubit, then its hermitian qubit acts as the target and which as the control and, conjugate U † can be applied to recover the original state, consequently, the circuit-diagram is sometimes drawn in since U †U = I resolves an identity operation. a symmetric fashion • CPHASE = • (66) C. Comparing classical and quantum gates The CNOT in terms of CPHASE can then be realized as The gate-sequences used to represent quantum algo- • rithms have certain similarities to those used in classi- CNOT = (67) cal computing, with a few striking differences. As an H • H example, we consider first the classical NOT gate (dis- cussed previously), and the related quantum circuit ver- Some two-qubit gates such as CNOT and CPHASE are sion, shown in Fig. 11. also called entangling gates, because they can take prod- While the classic bit-flip gate inverts any input state, uct states as inputs and output entangled states. They the quantum bit-flip does not in general produce the an- are thus an indispensable component of a universal gate tipodal state (when viewed on the Bloch sphere), but set for quantum logic. For example, consider two qubits rather exchange the prefactors of the wavefunction writ- A and B in the following state: ten in the computational basis. The X operator is some- 1 times referred to as ‘the quantum NOT’ (or ‘quantum |ψi = √ (|0i + |1i)A |0iB. (68) bit-flip‘), but we note that X only acts similar to the 2 classical NOT gate in the case of classical data stored in If we perform a CNOT gate, UCNOT, on this state, with the quantum bit, i.e. X|gi = |g¯i for g ∈ {0, 1}. qubit A the control qubit, and qubit B the target qubit, As briefly mentioned in Sec. IV B, all quantum gates the resulting state is (see the truth table in Fig. 10): are reversible, due to the underlying unitary nature of the operators implementing the logical operations. Cer- 1 tain other processes used in quantum information pro- UCNOT|ψi = √ (|0iA|0iB + |1iA|1iB) =6 (...)A(...)B, 2 cessing, however, are irreversible. Namely, measurements (69) (see Sec. V for detailed discussion) and energy loss to the which is a state that cannot be factored into an isolated the environment (if the resulting state of the environment qubit-A component and a qubit-B component. This is is not known). Here, we will not consider how these pro- one of the two-qubit entangled Bell states, a manifestly cesses are modeled, but refer the interested reader to e.g. quantum mechanical state. Ref. 178, and will only consider unitary control oper- 27 ations throughout the rest of this section. Finally, we as possible, and one wants to use as many of the native note that quantum circuits are written left-to-right (in gates as possible, to reduce the amount of time spent syn- order of application), while the calculation of the result thesising. Moreover, running a quantum algorithm also of a gate-sequences, e.g the circuit depends on the qubit connectivity of the device. The process of designing a quantum gate sequence that effi- |ψini U0 U1 ··· Un |ψouti (70) ciently implements a specific algorithm, while taking into account the considerations outlined above is known as is performed right-to-left, i.e. gate synthesis and gate compilation, respectively. A full discussion of this large research effort is outside the scope |ψouti = Un ··· U1U0|ψini. (71) of this review, but the interested reader may consult e.g. Refs. 183–185 and references therein as a starting point. As discussed in Sec. IV A, the NOR and NAND gates are As a concrete (and trivial) example of how gate identities each individually universal gates for classical computing. can be used, in Eq. (74) we illustrate how the Hadamard Since both of these gates have no direct quantum ana- gate from G1 can be generated by two single-qubit gates logue (because they are not reversible), it is natural to (from G0) and an overall phase gate, ask which gates are needed to build a universal quan- tum computer. It turns out that the ability to rotate about arbitrary axes on the Bloch-sphere (i.e. a complete 1 1 −1 −i 0 1 1 1 H = Ph π Y π Zπ = i√ = √ single-qubit gate set), supplemented with any entangling 2 2 2 1 1 0 i 2 1 −1 178,179 2-qubit operation will suffice for universality . By (74) using what is known as the ‘Krauss-Cirac decomposition’, any two-qubit gate can be decomposed into a series of As we show in Sec. IV D 1, the gates Xθ, Yθ and Zθ CNOT operations178,180. are all natively available in a superconducting quantum processor. We now address the question of how single qubit ro- 1. Gate sets and gate synthesis tations and two-qubit operations are implemented in transmon-based superconducting quantum processors. A common universal quantum gate set is
G0 = {Xθ, Yθ, Zθ, Phθ, CNOT} (72) 2. Addressing superconducting qubits
iθ where Phθ = e 1 applies an overall phase θ to a single The modes of addressing transmon-like superconduct- qubit. For completeness we mention another universal ing qubits can roughly be split into two main categories: gate set which is of particular interest from a theoretical i) Capacitive coupling between a resonator (or a feedline) perspective, namely and the superconducting qubit dipole-field allows for mi- crowave control to implement single-qubit rotations (see G1 = {H, S, T, CNOT}, (73) Sec. IV D) as well as certain two-qubit gates (see Sec- As a technical aside, we mention that the restriction to tions IV G and IV G 4). ii) For flux-tunable qubits, local a discrete gate set still gives rise to universality. This magnetic fields can be used to tune the frequency of indi- fact relies on using the so-called Solovay-Kitaev181,182 vidual qubits. This allows the implementation of z-axis theorem, which (roughly) states that any other single- single-qubit rotation as well as multiple two-qubit gates qubit gate can be approximated to an error using only (see Sections IV E, IV F and IV H). c O(log (1/)) (where c > 0) single-qubit gates from G1. The gate-set G1 is typically referred to as the ‘Clifford + T ’ set, where H, S and CNOT are all Clifford gates. D. Single-qubit gates Each quantum computing architecture will have cer- tain gates that are simpler to implement at the hard- In this section we will review the steps necessary to ware level than others (sometimes referred to as ’na- demonstrate that capacitive coupling of microwaves to a tive’ gates of the architecture). These are typically the superconducting circuit can be used to drive single-qubit gates for which the Hamiltonian governing the gate- gates. To this end we consider coupling a superconduct- implementation gives rise to a unitary propagator that ing qubit to a microwave source (sometimes referred to corresponds to the gate itself. We will show several exam- as a ‘qubit drive’) as shown in Fig. 12(a). A full circuit ples of this in Sections IV E, IV F, and IV G. Regardless analysis of the circuit in Fig. 12(a) is beyond the scope of which gates are natively available, as long as one has a of this review, so here we settle for highlighting the steps complete gate set, one can use the Solovay-Kitaev theo- that elucidate the physics of the qubit/drive coupling. rem to synthesize any other set efficiently. In general one The interested reader may consult a number of lectures wants to keep the overall number of time steps in which notes and pertinent theses (e.g. Refs. 44, 163, 186–188). gates are applied (denoted the depth of a circuit) as low Here we follow Ref. 163. 28
V (t) To elucidate the role of the drive, we move into a frame d Rw Cd rotating with the qubit at frequency ωq (also denoted ‘the rotating frame’ or the ‘the interaction frame’). To see the usefulness√ of this rotating frame, consider a state T |ψ0i = (1 1) / 2. By the time-dependent Schr¨odinger C equation this state evolves according to
iω t/2 room 1 e q |ψ0(t)i = UH0 |ψ0i = √ −iω t/2 , (80) temperature wiring on-chip 2 e q
FIG. 12. Circuit diagram of capacitive coupling of a mi- where UH0 is the propagator corresponding to H0. By crowave drive line (characterized by a time-dependent voltage calculating e.g. hψ0|σx|ψ0i = cos(ωqt) it is evident Vd(t)) to a generic transmon-like superconducting qubit. that the phase is winding with a frequency of ωq due to the σz term. By going into a frame rotating with the qubit at frequency ωq, the action of the drive can 1. Capacitive coupling for X,Y control be more clearly appreciated. To this end we define † U eiH0t U rf = = H0 and the new state in the rotating We start by modeling the qubit as an harmonic os- frame is |ψrf(t)i = Urf|ψ0i. The time-evolution in this cillator, for which the (classical) circuit Hamiltonian can new frame is again found from the Schr¨odingerequation be calculated by circuit quantization techniques, starting (using the shorthand ∂t = ∂/∂t), from Kirchoffs laws, and is given by163 i∂ |ψ t i i ∂ U |ψ i iU ∂ |ψ i Q˜ t 2 2 C t rf( ) = ( t rf) 0 + rf ( t 0 ) (81) ( ) Φ d ˜ H = + + Vd(t)Q, (75) ˙ † 2CΣ 2L CΣ = iUrfUrf|ψrfi + UrfH0|ψ0i (82) † † where CΣ = C + Cd is the total capacitance to ground = iU˙ rfU + UrfH0U |ψrfi (83) ˜ ˙ rf rf and Q = CΣΦ − CdVd(t) is a renormalized charge vari- | {z } able for the circuit. We can now promote the flux and He0 charge variables to quantum operators and assume weak coupling to the drive-line, so that Q˜ ≈ Qˆ, and arrive at We can think of the term He0 in the parentheses in C Eq. (83) as the form of H0 in the rotating frame. Sim- d ˆ H = HLC + Vd(t)Q, (76) ple insertion shows that He0 = 0 as expected (the ro- CΣ tating frame should take care of the time-dependence). 2 2 where HLC = Qˆ /(2C) + Φˆ /(2L) and we have kept However, one could also think of the term in brackets only terms that couple to the dynamic variables. Similar in Eq. (83) as a prescription for calculating the form of to the momentum operator for a harmonic oscillator in any Hamiltonian in the rotating frame given by Urf, by (x, p)–space, we can express the charge variable in terms replacing H0 with some other H. In general, we will not of raising and lowering operators, as done in Sec. II find He = 0. † Returning to Eq. (79), the form of Hd in the rotating Qˆ = −iQzpf a − a (77) frame is found to be p where Qzpf = ~/2Z is the zero-point charge fluctations p Hed = ΩVd(t) (cos(ωqt)σy − sin(ωqt)σx) . (84) and Z = L/C is the impedance of the circuit to ground. Thus, the LC oscillator capacitively coupled to a drive We can in general assume that the time-dependent part line can be written as, of the voltage (Vd(t) = V0v(t)) has the generic form † 1 Cd † H = ω a a + − Vd(t)iQzpf a − a . (78) v(t) = s(t) sin(ωdt + φ) (85) 2 CΣ = s(t) (cos(φ) sin(ωdt) + sin(φ) cos(ωdt)) , (86) Finally, by truncating to the lowest transition of the oscil- lator we can make the replacement a → σ− and a† → σ+ where s(t) is a dimensionless envelope function, so that throughout and arrive at the amplitude of the drive is set by V0s(t). Adopting the ωq definitions H = − σz + ΩVd(t)σy (79) 2 | {z } | {z } I = cos(φ) (the ‘in-phase’ component) (87) H H 0 d Q = sin(φ) (the ‘out-of-phase’ component) (88) § where Ω = (Cd/CΣ)Qzpf and ωq = (E1 − E0)/~.
§ 1 Starting from a generic qubit Hamiltonian, H0 = E0|0ih0| + H0 = ((E0 + E1)/2) − ((E1 − E0)/2)σz. In the main text we E1|1ih1|, we can rewrite in terms of Pauli matrices, and get have ignored the constant offset term. 29 the driving Hamiltonian in the rotating frame takes the (a) (b) form
LO AWG I I Q baseband Hed = ΩV0s(t)(I sin(ωdt) − Q cos(ωdt)) carrier Q pulses × (cos(ωqt)σy − sin(ωqt)σx) (89) Amplitude Performing the multiplication and dropping fast rotating ω s(t) LO LO terms that will average to zero (i.e. terms with ωq + ωd), ωAWG I Q (c) time known as the rotating wave approximation (RWA), we RF V t are left with d( ) I s(t) -x(- ) 1 -y(-Q) y(Q) Hed = ΩV0s(t) (−I cos(δωt) + Q sin(δωt)) σx x(I) 2 ω ω +ω d= LO AWG + (I sin(δωt) − Q cos(δωt)) σy (90) to qubits
FIG. 13. (a) Schematic of a typical qubit drive setup. A mi- where δω = ωq − ωd. Finally, by re-using the definitions crowave source supplies a high-frequency signal (ωLO), while from Eq. (86), the driving Hamiltonian in the rotating an arbitrary waveform generator (AWG) supplies a pulse- frame using the RWA can be written as envelope (s(t)), sometimes with a low frequency component, ωAWG, generated by the AWG. The IQ-mixer combines the two i(δωt+φ) Ω 0 e signals to generate a shaped waveform Vd(t) with a frequency Hed = − V0s(t) −i(δωt+φ) . (91) 2 e 0 ωd = ωLO ± ωAWG, typically resonant with the qubit. (b) Ex- ample of how a gate sequence is translated into a waveform Equation (91) is a powerful tool for understanding single- generated by the AWG. Colors indicate I and Q components. qubit gates in superconducting qubits. As a concrete (c) The action of a Xπ/2 pulse on a |0i state to produce the |−ii = √1 (|0i − i|1i) state. example, assume that we apply a pulse at the qubit fre- 2 quency, so that δω = 0, then
Ω Hed = − V0s(t)(Iσx + Qσy) , (92) 2 sequence of pulses (see Fig. 13(a)) Θk, Θk−1, ...Θ0 is con- verted to a sequence of gates operating on a qubit as showing that an in-phase pulse (φ = 0, i.e. the I- component) corresponds to rotations around the x-axis, k Y i [− 2 Θn(t)(Inσx+Qnσy )] while an out-of-phase pulse (φ = π/2, i.e. the Q- Uk ··· U1U0 = T e , (95) component), corresponds to rotations about the y-axis. n=0 As a concrete example of an in-phase pulse, writing out the unitary operator yields where T is an operator that ensures the pulses are gen- erated in the time-ordered sequence corresponding to Z t U ··· U U φ=0 i 0 0 k 1 0. Urf,d (t) = exp ΩV0 s(t )dt σx , (93) In Fig. 13 we outline the typical IQ modulation setup 2 0 used to generate the pulses used in Eq. (95). Fig. 13(a) which depends only on the macroscopic design parame- shows how a pulse at frequency ωd is generated using a ters of the circuit as well as the envelope of the baseband low phase-noise microwave generator (typically denoted pulse s(t) and amplitude V0, which can both be controlled ‘the local oscillator (LO)’), while the pulse is shaped by using arbitrary waveform generators (AWGs). Equation combining the LO signal in an IQ mixer with pulses gen- (93) is known as Rabi driving and can serve as a use- erated in an AWG. To allow for frequency multiplexing, ful tool for engineering the circuit parameters needed for the AWG signal will typically be generated with a low- efficient gate operation (subject to the available output frequency component, ωAWG, and the LO signal will be voltage V0). To see this we define the shorthand offset, so that ωLO+ωAWG = ωd. By mixing in more than one frequency ωAWG1, ωAWG2, ... it is possible to address Z t 0 0 multiple qubits (or readout resonators) simultaneously, Θ(t) = −ΩV0 s(t )dt (94) via the superposition of individual drives. 0 The I (Q) input of the IQ mixer will multiply the base- which is the angle by which a state is rotated given the band signal to the in-phase (out-of-phase) component of capacitive couplings, the impedance of the circuit, the the LO. In Fig. 13(b) we schematically show the com- magnitude V0, and the waveform envelope, s(t). This parison between XY gates in a quantum circuit and the means that to implement a π-pulse on the x-axis one corresponding waveforms generated in the AWG (omit- would solve the equation Θ(t) = π and output the sig- ting for clarity the frequency ωAWG component). The nal in-phase with the qubit drive. In this language, a inset in Fig. 13(b) shows an example of a gate on the 30
Bloch sphere, with indication of (I,Q) axes. More so- Finally we mention one more salient feature of the phisticated and compact approaches exist to reduce the virtual-Z gates. As shown in Ref.63, any single-qubit hardware needed for XY qubit control, relative to the operation (up to a global phase) can be written as setup shown in Fig. 13, see e.g.189–191. U θ, φ, λ Z π X π Z X π Z π , ( ) = φ− 2 2 π−θ 2 λ− 2 (100)
2. Virtual Z gate for appropriate choice of angles θ, φ, λ. This means X π that access to a single physical 2 combined with the virtual-Z gives access to a complete single qubit gate As we saw in Sec. IV D, the distinction between x– set! An explicit example of Eq. (100) in action is the and y–rotations was merely a choice of phase on the mi- Hadamard gate, which can be written as H = Z π X π Z π , crowave signals, and the angle to be rotated is given by 2 2 2 but since the Z’s can be virtual, it is possible to imple- Θ(t), both of which are generated using an AWG. Since ment Hadamards as an effective single pulse operation in the choice of phase φ has an arbitrary starting point, we superconducting qubits. could consider φ → φ + π/2. This would lead to I → Q and Q → −I. Therefore, changing the phase effectively changes rotations around x to rotations around y (and 3. The DRAG scheme vice-versa, with a change of sign). This is reminiscent of the result of applying a Zπ rotation to x– and y–rotations, In going from Eq. (78) to Eq. (79) we assumed we where ZπXπ = iYπ and ZπYπ = −iXπ. This analogy be- tween shifting a phase of an AWG-generated signal and could ignore the higher levels of the qubit. However, applying Z rotations can be utilized to implement virtual for weakly anharmonic qubits, such as the transmon (see 192 Sec. II), this may not be a justified assumption, since Z gates . As shown by McKay et al., this intuition 1→2 0→1 ωq only differs from ωq(≡ ωq ) by the anharmonicity, can be formalized via the following example: consider 1→2 the case of applying a pulse with an angle θ on the I α = ωq − ωq, which is negative and typically around channel (i.e. a Xθ) followed by another θ pulse on the I 200 to 300 MHz. This situation is sketched in Fig. 14(a- channel, but with a phase φ0 relative to the first pulse c), where we illustrate how Gaussian pulses with stan- (φ0) dard deviations σ = {1, 2, 5} ns have spectral content (denoted Xθ , where X indicates we still use the I chan- 1→2 that leads to non-zero overlaps with the ωq = ωq − |α| nel, but the rotation axis is now an angle φ0 away from 1 the x-axis). Using Eq. (95) this corresponds to a pulse frequency. This leads to two deleterious effects: ( ) leak- sequence age errors which take the qubit out of the computational subspace, and (2 ) phase errors. Effect 1 can occur be- θ (φ0) −i (cos(φ )σ +sin(φ )σ ) 2 0 x 0 y cause a qubit in the state |1i may be excited to |2i as a π Xθ Xθ = e Xθ (96) pulse is applied, or be excited directly from the |0i, since = Z−φ XθZφ Xθ (97) 0 0 the qubit spends some amount of time in the |1i state from which we see that the effect of the offset phase φ0 during the π pulse. Effect 2 occurs because the presence | i | i is to apply Zφ0 . The equality above can be verified with of the drive results in a repulsion between the 1 and 2 ω0→1 a little trigonometric footwork. The final Z−φ0 is due to levels, in turn changing q as the pulse is applied. This the rotation being in the frame of reference of the qubit. leads to the accumulation of a relative phase between 194 195–197 However, since readout is along z-axis (see Sec. V), a |0i and |1i . The so-called DRAG procedure final phase rotation about z will not change the mea- (Derivative Reduction by Adiabatic Gate) seeks to com- surement outcome. Thus, if one wants to to implement bat these two effects by applying an extra signal in the the gate sequence out-of-phase component. The trick is to modify the wave- form envelope s(t) according to ··· U Zθ Ui+1 Zθ Ui+2 ··· (98) i 0 1 s(t) on I 0 where Ui’s are arbitrary gates, this can be done by re- s(t) → s (t) = s˙(t) , (101) vising the gate sequence (in the control software for the λ on Q α AWG) and changing the phase of subsequent pulses where λ is a dimensionless scaling parameter, and λ = (θ0) (θ0+θ1) ··· Ui Ui+1 Ui+2 ··· (99) 0 correspond to no DRAG pulse ands ˙(t) is the time derivative of s(t). The theoretically optimal choice for which reduces the number of overall gates. Moreover, reducing dephasing error is λ = 0.5 and an optimal choice the virtual-Z gates are “perfect”, in the sense that no for reducing leakage error is λ = 1196,198. Interchanging additional pulses are required, and the gate takes “zero I and Q in Eq. (101) corresponds to DRAG pulsing for time”, and thus the gate fidelity is nominally unity. As the Q component. we show in Sections IV E and IV F, operation of two- In practice there can be a deviation from these two qubit gates can incur additional single-qubit phases. Us- optimal values, often due to pulse distortions in the lines ing the virtual-Z strategy, these phases can be cancelled leading to the qubits. Typically, randomized benchmark- out, leaving a pure two-qubit interaction. ing experiments combined with single-shot measurements 31
(a) (b) σ (ns) quency), i.e. 2E1 α 1 s0 (t) = s0(t)ei2πδft, (102) E 2 δf 2 5 and choosing λ to minimize leakage errors, then phase 199
Ampl. (a.u.) errors can be reduced simultaneously . Similarly, by a judicious use of the virtual-Z gate, it is also -20 -10 0 10 20 possible to reduce phase errors in combination with E1 192 (c) Time (ns) DRAG pulsing to reduce leakage . Modern single-qubit FFT gates using DRAG pulsing now routinely reach fidelities 65,67,199,202–205 F1qb & 0.99 . Other techniques also exist for operating single-qubit gates in a spectrally crowded device206,207. 01 α Ampl. (a.u.) ωq E0 E. The iSWAP two-qubit gate in tunable qubits -0.5 -0.25 0 0.25 0.5 Frequency (GHz) (d) (e) time As briefly mentioned in Sec. IV C, single-qubit gates λ = 0 I supplemented with an entangling two-qubit gate can Q form the gate set required for universal quantum compu- tation. The two-qubit gates available in the transmon- like superconducting qubit architecture can roughly be
Ampl. (a.u.) split into two broad families as outlined previously: one group requiring local magnetic fields to tune the transition frequency of qubits and one group consist- (f) time (ns) (g) ing of all-microwave control. There exist several hy- λ = 0.5 I time brid schemes that combine various aspects of these two Q categories and, in particular, the notions of tunable coupling and parametric driving are proving to be im- portant ingredients in modern superconducting qubit processors63,67,91,104,107,208–214. In this section, however, Ampl. (a.u.) we start by introducing the iSWAP gate, and then review the CPHASE (controlled-phase) in Section IV F and the time (ns) CR (cross-resonance) in Section IV G. We briefly review a few other two-qubit gates and discuss their merits in FIG. 14. (a) Schematic level diagram of a weakly anharmonic Sections IV G 4 and IV H. transmon qubit subjected to a drive at transition frequency ωd = ωq. (b) Gaussian waveform with standard deviation σ. (c) Fourier transform of (b) showing how the short pulse 1. Deriving the iSWAP unitary 1→2 lengths lead to significant overlap with the ωq transition, separated from ωq by the anharmonicity α. (d) Waveform of a Xπ pulse without DRAG modulation. (e) Effect of the As we saw in Sec. II, Eq.32 the interaction term be- waveform from (d) on a qubit initialized in the |0i state with tween two capacitively coupled qubits (in the two-level α = −200 MHz and ωq = 4 GHz. The dephasing error is approximation) is given by visible as a deviation from the |1i after the pulse. (f) Wave- Hqq = gσy1 ⊗ σy2, (103) form of a Xπ pulse with DRAG modulation for a qubit with anharmonicity α = −200 MHz and DRAG parameter λ = 0.5 where g is the coupling strength and ⊗ is used to em- to cancel dephasing errors (see text for details). (g) Effect of phasize the tensor product. If the capacitive coupling is the waveform from (f) on the same qubit as (e). Calculated 215,216 using mesolve in the software package QuTiP193. mediated through a bus resonator, then
g1g2(∆1 + ∆2) g → gq-r-q = , (104) 2∆1∆2
where gi is the resonator coupling to qubit i (depen- | i (see Sec. V) of the 2 state is used to determine the op- dent on the qubit-resonator coupling capacitance Cqir) λ λ { . , } timal value of . The = 0 5 1 tradeoff was demon- and ∆i = ωqi − ωr is the detuning of qubit i to the res- strated explicitly in192,199. However, by extending the 200,201 onator. In the simpler case where the qubits are directly original DRAG pulse implementation , it is is possi- coupled217, ble to reduce both errors simultaneously. By introducing a frequency detuning parameter δf to the waveform196 1√ Cq-q g → gq-q = ωq1ωq2 p p , (105) (defined such that δf = 0 corresponds to qubit fre- 2 Cq-q + C1 Cq-q + C2 32 where Cq-q is the qubit-qubit coupling capacitance and (a) initial point Ci is the capacitance of qubit i. Throughout this sec- final point tion, we will assume a direct capacitive coupling between qubits of the flux-tunable transmon type, so that g = gq-q 5.0 and ωqi → ωqi(Φi). For simplicity, we suppress the ex- (GHz) plicit flux dependence of the ωqi’s and simply refer to the coupling as g. Equation (103) can be rewritten as