Normal-Metal Quasiparticle Traps For Superconducting Qubits: Modeling, Optimization, and Proximity Effect
Von der Fakultät für Mathematik, Informatik und Naturwissenschaften der RWTH Aachen University zur Erlangung des akademischen Grades eines Doktors der Naturwissenschaften genehmigte Dissertation
vorgelegt von Amin Hosseinkhani, M.Sc.
Berichter: Universitätsprofessor Dr. David DiVincenzo, Universitätsprofessorin Dr. Kristel Michielsen
Tag der mündlichen Prüfung: March 01, 2018
Diese Dissertation ist auf den Internetseiten der Universitätsbibliothek online verfügbar.
Metallische Quasiteilchenfallen für supraleitende Qubits: Modellierung, Optimisierung, und Proximity-Effekt
Kurzfassung: Bogoliubov Quasiteilchen stören viele Abläufe in supraleitenden Elementen. In supraleitenden Qubits wechselwirken diese Quasiteilchen beim Tunneln durch den Josephson- Kontakt mit dem Phasenfreiheitsgrad, was zu einer Relaxation des Qubits führt. Für Tempera- turen im Millikelvinbereich gibt es substantielle Hinweise für die Präsenz von Nichtgleichgewicht- squasiteilchen. Während deren Entstehung noch nicht einstimmig geklärt ist, besteht dennoch die Möglichkeit die von Quasiteilchen induzierte Relaxation einzudämmen indem man die Qu- asiteilchen von den aktiven Bereichen des Qubits fernhält. In dieser Doktorarbeit studieren wir Quasiteilchenfallen, welche durch einen Kontakt eines normalen Metalls (N) mit der supraleit- enden Elektrode (S) eines Qubits definiert sind. Wir entwickeln ein Modell, das den Einfluss der Falle auf die Quasiteilchendynamik beschreibt, wenn überschüssige Quasiteilchen in ein Trans- monqubit injiziert werden. Dieses Modell ermöglicht es, unter Berücksichtigung der Fallenpa- rameter die Zeitskala zu bestimmen, in der die überschüssigen Quasiteilchen aus dem Transmon evakuiert werden. Wir zeigen, dass die Evakuierungsdauer monoton mit der Fallengröße ansteigt und letztlich auf einen Grenzwert zuläuft, der von der Quasiteilchen-Diffusionskonstante und von der Qubitgeometrie abhängt. Wir errechnen die charakteristische Fallengröße, bei welcher dieser Grenzwert erreicht wird. Wie es sich herausstellt, ist der limitierende Faktor für die Einfangrate der Falle durch die langsame Quasiteilchenrelaxation im normalen Metall gegeben; diese Relaxation ist jedoch nur schwer kontrollierbar. Um das Einfangen von Quasiteilchen zu optimieren, studieren wir den Einfluss von Größe, Anzahl und räumlicher Anordnung der Fallen. Diese Faktoren sind insbesondere wichtig, wenn die Falle die charakteristische Größe überschreitet. Wir diskutieren für einige experimentell rel- evante Beispiele wie die Evakuierungsdauer der überschüssigen Quasiteilchen optimiert werden kann. Darüberhinaus zeigen wir, dass eine Falle nahe des Josephson-Kontaktes die stationäre Quasiteilchendichte an demselben Kontakt unterdrückt und den Einfluss von Fluktuationen der Quasiteilchenerzeugung reduziert. Wenn metallische Elemente an ein supraleitendes Material gekoppelt sind, können Cooper- paare ins Metall entweichen. Mit dem Usadelformalismus greifen wir zunächst den Proximity- Effekt von gleichförmigen NS-Doppelschichten wieder auf; trotz der bereits langjährigen Er- forschung dieses Problems erlangen wir zu neuen Erkenntnissen über die Zustandsdichte. Wir verallgemeinern unsere Resultate danach für das ungleichförmige Problem in der Nähe der Fal- lenkante. Durch die Kombination dieser Resultate mit dem davor entwickelten Modell zur Unterdrückung der Quasiteilchendichte finden wir einen optimalen Abstand zwischen Falle und Josephson-Kontakt in einem Transmonqubit, welcher zu einer Minimierung der Qubitrelaxation führt. Dieser optimale Abstand, der die 4- bis 20-fache Kohärenzlänge beträgt, resultiert aus dem Wechselspiel zwischen Proximity-Effekt und Unterdrückung der Quasiteilchendichte. Wir schließen daraus, dass der schädliche Einfluss des Proximity-Effekts umgangen werden kann solange die Entfernung zwischen Falle und Kontakt größer als der optimale Abstand ist.
Normal-Metal Quasiparticle Traps for Superconducting Qubits: Modeling, Optimization, and Proximity Effect
Abstract: Bogoliubov quasiparticle excitations are detrimental for the operation of many su- perconducting devices. In superconducting qubits, quasiparticles interact with the qubit degree of freedom when tunneling through a Josephson junction, and this interaction can lead to qubit relaxation. At millikelvin temperatures, there is substantial evidence of nonequilibrium quasi- particles. While there is no agreed upon explanation for the origin of these excess quasiparticles, it is nevertheless possible to limit the quasiparticle-induced relaxation by steering quasiparticles away from qubit active elements. In this thesis, we study quasiparticle traps that are formed by a normal-metal in tunnel contact with the superconducting electrode of a qubit. We develop a model to explain how a trap can influence the dynamics of the excess quasiparticles injected in a transmon-type qubit. This model makes it possible to find the time it takes to evacuate the injected quasiparticles from the transmon as a function of trap parameters. We show when the trap size is increased, the evacuation time decreases monotonically and saturates at a level that depends on the quasiparticles diffusion constant and the qubit geometry. We find the charac- teristic trap size needed for the evacuation time to approach the saturation value. It turns out that the bottleneck limiting the trapping rate is the slow quasiparticle energy relaxation inside the normal-metal trap, a quantity that is very hard to control. In order to optimize normal-metal quasiparticle trapping, we study the effects of trap size, number, and placement. These factors become important when the trap size increases beyond the characteristic length. We discuss for some experimentally relevant examples how to shorten the evacuation time of the excess quasiparticle density. Moreover, we show that a trap in the vicinity of a Josephson junction can significantly suppress the steady-state quasiparticle density near that junction and reduce the impact of fluctuations in the generation rate of quasiparticles. When such normal-metal elements are connected to a superconducting material, Cooper- pairs can leak into the normal-metal trap. This modifies the superconductor properties and, in turn, affects the qubit coherence. Using the Usadel formalism, we first revisit the proximity effect in uniform NS bilayers; despite the long history of this problem, we present novel findings for the density of states. We then extend our results to describe a non-uniform system in the vicinity of a trap edge. Using these results together with the previously developed model for the suppression of the quasiparticle density due to the trap, we find in a transmon qubit an optimum trap-junction distance at which the qubit relaxation rate is minimized. This optimum distance, of the order of 4 to 20 coherence lengths, originates from the competition between proximity effect and quasiparticle density suppression. We conclude that the harmful influence of the proximity effect can be avoided so long as the trap is farther away from the junction than this optimum.
Acknowledgements
Before starting the main part of the thesis, I would like to take the opportunity to express my deepest gratitude to several people who helped and supported me during my doctoral studies. First of all, I am greatly thankful to my PhD supervisor Dr. Gianluigi Catelani for his kind and constant support. I am indebted for countless enjoyable discussions with him during which I learned to focus on Physics behind the mathematics. I also like to thank him for sending me to a lot of conferences and schools where I got a chance to expand my knowledge and to meet and discuss with a lot of physicists. I also enjoyed a lot by collaborating with experimental scientists at Yale University, which was made possible by Gianluigi. I am grateful to my friend Dr. Roman-Pascal Riwar for a lot of scientific as well as every- day-life discussions that we had. He also kindly helped me in translating the abstract of this thesis into (Swiss) German. I like to express my deep gratitude to Prof. David DiVincenzo for reviewing my thesis, his support for my postdoc applications, and the friendly and relaxed atmosphere at the JARA- Institute of Quantum Information. I would like to thank the second reviewer of my thesis Prof. Kristel Michielsen and also other members of my PhD committee, Prof. Thomas Schäpers and Prof. Christoph Stampfer, for the time they spent on reading my thesis and their fruitful comments. I would also like to thank all of my colleagues in JARA-Institute of Quantum Information at Forschungszentrum Jülich and RWTH Aachen University; few of which include Dr. Daniel Zeuch for a lot of discussions we had almost every day and also helping me in the German abstract of the thesis, Dr. Sbastian Mehl for helping me with a lot of paper works at the time I just had started my doctoral study in Germany, Alessandro Ciani for kindly preparing my PhD graduation hat, and Alwin van Steensel for a lot of discussions and his kind PhD gift. I like to particularly thank Dr. Mohammad H. Ansari for lots of interesting discussions that we had during my doctoral studies and also his kind support for my postdoc application to his research group. I am very thankful to Ms. Luise Snyders who helped me very much for handling official paper works at Forschungszentrum Jülich. I also like to thanks Ms. Helene Barton for helping me in doing paper works at RWTH Aachen University. I like to thank Dr. Hamed Saberi who supervised me during my Master’s program at Shahid Beheshti University, Tehran, Iran and supporting my PhD applications. I would like to thank Prof. Ali Rezakhani for co-advising my master’s thesis and also for inviting me to visit Institute for Research in Fundamental Sciences (IPM), Tehran, Iran in September 2016. I also thank Dr. Sahar Alipour for scheduling that visit. I would like to thank Prof. Maksym Serbyn for inviting me to visit his group at Institute of Science and Technology Austria in February 2018. I also like to thank Prof. Johannes Fink and Dr. Shabir Barzanjeh for a lot of interesting discussions that I had with them and their hospitality during my visit to IST Austria. I am thankful to Dr. Frank Deppe for inviting me to visit Walther-Meißner-Institute for Low-Temperature Research in March 2018. I greatly appreciate his hospitality and a lot of interesting discussions that we had. I am grateful to Prof. Ahmad Ghodsi Mahmoudzadeh, my advisor during my bachelor studies at Ferdowsi University of Mashhad, Iran, for his constant support and hospitality. vi
I also like to thank all of my friends at Forschungszentrum Jülich whom I share a lot of sweet memories: Esmaeel, Amin, Keyvan, Ali, Vahid, Masood and Davood. I like to especially thank my friends at Aachen: Mojataba and his wife Shima, Alireza and his wife Zeinab, Sina and Mehrdad. Our very nice memories made my life so colourful while I was living far from my family. Finally and most importantly, I would like to express my deepest gratitude to my beloved parents who helped me at each second of my life. I am greatly thankful to my brother Yasin who has always supported me in my life and my education. At the very last days of finalizing this thesis, I was utterly joyed by born of my beloved nephew, Omid. I wish him all the best in his life, and a happy family forever. Contents
1 Introduction 1 1.1 Overview...... 1 1.2 Outline...... 2
2 Quantum Coherent Superconducting Devices5 2.1 Superconductivity...... 5 2.1.1 Bogoliubov approach to BCS superconductivity...... 6 2.1.2 Quasiparticle density in thermal equilibrium...... 7 2.1.3 Superconducting gap...... 8 2.2 Josephson Effect...... 8 2.3 Superconducting Qubits...... 10 2.3.1 Cooper-pair box...... 10 2.3.2 Transmon qubit...... 13 2.4 Qubit-Quasiparticle Interaction...... 13 2.4.1 Energy relaxation induced by quasiparticle tunneling...... 14
3 Normal-Metal Quasiparticle Traps 19 3.1 Modeling...... 19 3.1.1 Introduction...... 19 3.1.2 The diffusion and trapping model...... 20 3.1.3 Quasiparticle dynamics during injection and trapping...... 24 3.1.4 Experimental data...... 29 3.2 Optimization...... 32 3.2.1 Introduction...... 32 3.2.2 Enhancing the decay rate of the density...... 33 3.2.3 Suppression of steady-state density and its fluctuations...... 40 3.3 Summary and Conclusions...... 45
4 Quasiclassical Theory of Superconductivity 47 4.1 Gor’kov Equations...... 47 4.2 Dyson Equation in Keldysh-Nambu Space...... 50 4.3 Eilenberger Equations...... 52 4.4 The Dirty Limit...... 56 4.4.1 Boundary conditions for proximity systems...... 58 4.4.2 Usadel equations for normal-superconducting hybrids...... 59
5 Proximity Effect in Normal-Metal Quasiparticle Traps 61 5.1 Introduction...... 61 5.2 Qubit relaxation due to quasiparticles...... 62 5.3 Proximity effect in thin films...... 63 5.3.1 Uniform NS bilayers...... 63 viii Contents
5.3.2 Proximity effect near a trap edge...... 66 5.4 Qubit relaxation with a trap near the junction...... 70 5.4.1 Thermal equilibrium...... 70 5.4.2 Suppressed quasiparticle density...... 72 5.5 Summary...... 77
6 Summary and Conclusions 79
A Tunneling rate equations 81
B Derivation of effective trapping rate 83 B.1 Thin normal metal...... 83 B.2 Effective trapping rate...... 85 B.3 Effective trapping rate integrated over energy...... 87
C Comparison with vortex trapping 89
D Finite-size trap 91
E Quasi-degenerate modes and their observability 95
F Effective length 99 F.1 Effective length due to the pad...... 99 F.2 Effective length due to the gap capacitor...... 100
G Quasiparticle Decay Rate and Steady-State Density 101 G.1 Slowest Quasiparticle Decay Rate Due To Trap...... 101 G.1.1 Single trap...... 101 G.1.2 Multiple side traps...... 103 G.2 Suppression of Quasiparticle Steady-State Density...... 105
H Traps in the Xmon geometry 107
I Proximity effect in uniform NS bilayers 109 I.1 Weak-coupling limit...... 109 I.2 Strong-coupling limit...... 111
J Numerical solution of the self-consistent equation for the order parameter 115
K Spatial evolution of single-particle density of states and pair amplitude 119
L Spectral function in the presence of a trap 123 L.0.1 Thermal equilibrium...... 123 L.0.2 Suppressed quasiparticle density...... 125
Bibliography 127 Chapter 1 Introduction
1.1 Overview
In quantum world, a physical system can be in different eigenstates simultaneously. In partic- ular, one can imagine a quantum system in a superposition of two states,
θ θ |Ψi = cos |1i + sin eiφ|0i, (1.1) 2 2 where we refer such quantum system as qubit; the qubit state can be conveniently illustrated by the Bloch sphere drawn in figure (1.1). In 1994, Peter Shor devised a quantum algorithm running on a set of qubits that enables to factorize an integer N to prime factors in an exponentially faster time compared with classical algorithms using binary logic [1]. Such an algorithm can, for example, be used to break public key cryptography protocols. Since this discovery, there has been extraordinary efforts to exploit the potentials of quantum systems for developing future technologies and also in finding physical systems promising to be controlled as a qubit and feasible to scale up. As any system interacts with its surrounding environment, a qubit that is initially prepared in the superposition of two states will eventually decohere into a mixed state that manifests itself in vanishing of the non-diagonal elements of qubit density matrix. In addition to qubit decoherence, faulty quantum preparation, faulty measurements and faulty quantum gates impose an obstacle for realizing quantum computation. For these reasons, performing quantum error correction is at the heart of any scheme for realizing a universal quantum computer, that in turn requires the coherence time of physical qubits compared with gate time to exceed a threshold [2]. Further improvement of the qubit coherence above the threshold is always favorable as it reduces the computational overhead due to quantum error correction. Indeed, maintaining the qubit coherence, denoted by T2, for long enough times
Figure 1.1: Qubit state on a Bloch sphere. All points on the surface correspond to a superpo- sition between states |0i and |1i. 2 Chapter 1. Introduction is a key issue in all aspects of quantum technologies. There are generally two processes that contribute to the qubit decoherence. First, energy relaxation denoted by T1 that is an irreversible process giving energy from or to the qubit which results in qubit state transition to the ground or excited state. Second, pure dephasing denoted by Tφ, is due to perturbations that do not induce qubit state change, but randomly modulate the qubit phase. These two process combine to give, 1 1 1 = + . (1.2) T2 2T1 Tφ A large number of different systems have been proposed to realize physical qubits, each of which utilize a specific physical property of the system to encode quantum information; for example, intrinsic spin degree of freedom, photon polarization or phase difference across a superconducting tunnel junction. Depending on what physical property is used to build a qubit, there are different decoherence mechanisms that are relevant for the system. In this thesis, we particularly focus on superconducting qubits that are among the most promising candidates for realization of quantum computation. There has been an extensive research in the community to control and suppress various decoherence mechanisms leading to nearly six orders of magnitude improvement of coherence time in the past 20 years [3]. In particular, designing qubits that operate in a regime robust against charge noise [4,5] together with improved understanding and control of dielectric losses [6] and the Purcell effect [7] have made possible to reaching coherence time more than 100µs. As the gate operation time for superconducting qubits is of order tens of nanoseconds [8], such long coherence time has enabled to perform about 104 operations per error [3]. While the mentioned decoherence factors are imposed on the qubit from the surrounding environment, there is also an intrinsic decoherence channel that originates from the supercon- ductor itself: tunneling of unpaired electrons or quasiparticles across the Josephson junction. The theory of quasiparticle-induced decoherence is already developed and quasiparticle effects on relaxation, dephasing and parameter renormalization has been studied both theoretically [30, 31, 32, 33, 34, 35, 36, 37] and experimentally [18, 38, 39, 40, 41, 42, 43, 44, 45, 46]. While the generation mechanism of non-equilibrium quasiparticles has yet remained a mystery in the field, in this thesis we study ways to suppress quasiparticle density in order to minimize their detrimental effects. We illustrate how planting normal metals over some parts of the qubit can trap quasiparticles and show how such trapping process can be optimized. Moreover, we study inverse superconducting proximity effect and discuss how its detrimental consequences to the qubit operation can be avoided. While we mostly consider a 3D transmon qubit and give some consideration of an Xmon qubit to explain and discuss our proposal, the idea of controlling quasiparticle population is also important in other superconducting devices such a Cooper-pair pumps, turnstiles and possible topological qubits based on Majorana zero modes where quasi- particle poisoning is a major obstacle. We therefore hope that our work will be of importance and finds applications in those communities as well.
1.2 Outline
The thesis is structured as follows: In chapter2 we review the necessary backgrounds; section 2.1 is devoted to remind the reader about formation of superconductivity and showing how Chapter 1. Introduction 3 quasiparticles appear in the formalism. Section 2.2 uses perturbation theory to discuss about the Josephson effect. In section 2.3 we explain about two types of superconducting qubits, Cooper- pair box and transmon. We finally discuss in section 2.4 how quasiparticle tunneling results in qubit energy relaxation. Chapter3 contains two sections that covers our published papers cited in references [47] and [68]. In section 3.1 we develop a phenomenological model governing the dynamics and steady state density of quasiparticles. Considering a 3D-transmon qubit, this model allows one to evaluate the time it takes to evacuate the injected quasiparticles from the transmon as a function of trap parameters. With the increase of the trap size, this time decreases monotonically, saturating at the level determined by the quasiparticles diffusion constant and the qubit geometry. We determine the characteristic trap size needed for the evacuation time to approach that saturation value. We also present experimental data (obtained by our colleagues) that support our theoretical findings. In section 3.2 we discuss how normal-metal quasiparticle traps can be optimized. We consider some experimentally relevant examples and find optimum trap configurations that maximize the decay rate of excess quasiparticle density. Moreover, we show that a trap in the vicinity of a Josephson junction can significantly reduce the steady-state quasiparticle density near that junction, thus suppressing the quasiparticle-induced relaxation rate of the qubit. Such a trap also reduces the impact of fluctuations in the generation rate of quasiparticles, rendering the qubit more stable in time. We then turn our attention to study how a normal metal in contact with a superconductor can modify superconducting properties and what are the related consequences for the qubit coherence. In chapter4, we review the quasicalssical theory of superconductivity in terms of Green’s functions in Keldysh-Nambu space and derive the Eilenberger equation. We then consider the dirty limit and derive the Usadel equation for a normal-superconducting hybrid that is our starting point to find the influence of a normal-metal trap on the quasiparticle density of states. Chapter5 contains our published paper cited in reference [78]. Here, we first apply Usadel theory to a uniform bilayer and find new analytical corrections to the density of states and minigap energy. We then consider a non- uniform hybrid relevant for physical realization of normal-metal trap for superconducting qubits. We find both theoretically and numerically the density of states and pair amplitude as a function of distance from trap edge. Moreover, building on the phenomenological diffusion equation that we develop in chapter3, we model the effect of the trap on the quasiparticle distribution function. This enables us to calculate different contributions to the qubit relaxation induced by quasiparticle tunneling and pair processes. We find an optimum trap-junction distance that minimizes the qubit relaxation rate. Placing the trap further from the junction than this optimum distance ensures that inverse proximity effect does not harm qubit coherence. We summarize and conclude our work in the last chapter. We have included a number of Appendices as well that complement the main text and present some details of calculations.
Chapter 2 Quantum Coherent Superconducting Devices
In this chapter we review the Bogoliubov approach to conventional superconductivity followed by discussing the Josephson effect and then two specific types of superconducting qubits, Cooper-pair box and transmon. We then consider qubit-quasiparticle interaction and present how quasiparticle tunneling results in qubit energy relaxation.
2.1 Superconductivity
Superconductivity was discovered in 1911 by H.K. Onnes when he witnessed electric resistance of solid mercury suddenly dropping below any measurable value by cooling down to 4.2 K. Later, it was discovered by Meissner and Ochsenfeld that superconductors are diamagnets as well [9] so that the electromagnetic field is expelled from a bulk superconductor beyond the material- and temperature-dependent penetration length. The microscopic explanation about superconductivity remained a challenge until 1957 when Bardeen, Cooper and Schrieffer pro- posed a model (BCS) that successfully explains this phenomenon [10]. The key idea of this model is that electrons condensate into a coherent state of pairs. The challenging question here is how two electrons in a lattice can overcome the repulsive Coulomb interaction between them. This can be explained by taking into account the motion of ion cores or phonons; the first elec- tron polarizes its surrounding medium by attracting ion cores; the resulting positive ions can then attract the second electron. While the importance of this electron-lattice interaction was first pointed out in 1950 by Fröhlich [11], in 1956 Cooper showed however weak the (phonon- mediated) attractive interaction between two electrons is, the Fermi surface is unstable against formation of a Cooper pair [12]. There are different approaches to BSC superconductivity, these include variational method, the Bogoliubov approach [13] and the Gor’kov approach [14]. The original BSC paper uses the variational approach for which a trial wave function for the ground state of a superconductor up to the global phase is taken as:
Y iφ † † |ΨBCSi = (|uk| + e |vk|ck↑c−k↓)|0i (2.1) k
† where ckσ is the electron creation operator with momentum k and spin σ, |0i is the vacuum, the product is taken over all one-electron states and the coherence factors, uk and vk, are complex numbers that satisfy normalization condition:
2 2 |uk| + |vk| = 1. (2.2) 6 Chapter 2. Quantum Coherent Superconducting Devices
Indeed, the BSC ground state is clearly a coherent superposition of all states with even number of electrons from zero to infinity with an arbitrary phase factor eiφ. The model Hamiltonian in BCS theory is, X † X † † HBCS = ξkckσc−kσ + λkk0 ck↑c−k↓c−k0↓ck0↑, (2.3) kσ kk0 in which, the energy ξk is measured from the Fermi energy. Since the attractive interaction is mediated by phonons, the coupling strength can be taken constant, λkk0 = −λ, when both the 0 scattering-in and out electrons possess energies less than Debye frequency, |ξk|, |ξk| < ωD, and is zero otherwise. This model Hamiltonian describes the interaction between Cooper pairs and leads to formation of Cooper pair condensate as the ground state for a superconductor. The idea of variational approach to BSC superconductivity is to minimize the expectation value of the model Hamiltonian using the trial ground state given by Eq. (2.1). This enables us to find the coherence factors in a straightforward calculation that we do not present here. Rather, in this section we review the Bogoliubov approach that also gives the coherence factors and provides a handy framework to deal with superconducting excitations as well. In chapter4 we start from Gor’kov approach and present quasiclassical theory of superconductivity to study modifications in superconducting properties when the superconductor is in proximity to a normal metal.
2.1.1 Bogoliubov approach to BCS superconductivity To begin with, we note that since the ground state of a superconductor is formed from Cooper pair condensate including a macroscopic number of Cooper pairs, adding or removing an extra Cooper pair does not really matter. In other words, the value for anomalous average defined as hc−k0↓ck0↑i is finite (averaging is taken with respect to the ground state of a superconductor) and fluctuations around this finite value are small . This enables us to arrive to the mean-field BCS Hamiltonian: |∆2| HMF = X ξ c† c − X(∆c† c† + ∆∗c c ) − . (2.4) BCS k kσ −kσ k↑ −k↓ −k↓ k↑ λ kσ k Here we only consider s-wave superconductivity where the order parameter is everywhere sym- metric in k space and is given by, X ∆ = λhc−k0↓ck0↑i. (2.5) k0 We now introduce the Bogoliubov-quasiparticle operator defined as ∗ † γk↑ = ukck↑ + vkc−k↓, (2.6) ∗ † γ−k↓ = −vkck↑ + ukc−k↓, (2.7) where the normalization condition for coherence factors, Eq. (2.2), is derived by enforcing † fermionic anticommutation relations for Bogoliubov operators, {γkσ, γk0σ0 } = δσ,σ0 δk,k0 . The mean-field BCS Hamiltonian is diagonal in basis of Bogoliubov operators provided the following relations for coherence factors hold, 2 1 ξk |uk| = 1 + , (2.8a) 2 εk 2 1 ξk |vk| = 1 − , (2.8b) 2 εk Chapter 2. Quantum Coherent Superconducting Devices 7 where q 2 2 εk = ξk + |∆ |. (2.9) In finding these relations, it also turns out that the phase of superconducting order parameter ∆ is equal to phase of vk relative to uk so that the order parameter has the same phase as the BSC ground state. The Hamiltonian then becomes,
MF HBCS = HG + Hqp, (2.10) for which we defined,
|∆2| H = X(ξ − ε ) − , (2.11) G k k λ k and,
X † † Hqp = εk(γk↑γk↑ + γ−k↓γ−k↓). (2.12) k
Assuming a normal state at T = 0, we have ∆ = 0 and εk = |ξk|. The first term in Eq. (2.10) then differs from the corresponding one in the normal phase by,
|∆2| H − HN (T = 0) = 2 X (ξ − ε ) − . (2.13) G G k k λ k>kf By changing the summation to an integration, it is easy to simplify this energy difference to 1 2 − 2 N0|∆ | that is the condensation energy and expresses the energy gain by forming supercon- ductivity. What is important for us is the second term in Eq. (2.10), Hqp, that illustrates the energy increase corresponding to quasiparticle excitation above the Cooper pair condensate. Indeed, one can directly check from Eq. (2.1) that the superconducting ground state is vacuum state for quasiparticle excitations, γk|ΨBCSi = 0, that are gapped from the ground state condensate by the value determined by the order parameter. Later in this chapter we will show that the density of quasiparticle excitations has an important role in energy relaxation of superconducting qubits. In the following we find this density in thermal equilibrium.
2.1.2 Quasiparticle density in thermal equilibrium The Bogoliubov transformation makes it clear that there is a one-to-one correspondence between electronic and quasiparticle excitations. We can therefore write,
Nqp(ε)dε = Ne(ξ)dξ, (2.14) where the quasiparticle density of states is denoted by Nqp(ε) and electronic density of states by Ne(ξ). As we are interested in energies close to Fermi level, we take Ne(ξ) ' Ne(ξF ) ≡ N0 and find the normalized quasiparticle density of states for a bulk superconductor, N (ε) dξ ε n(ε) = qp = = Re √ sgn(ε). (2.15) N0 dε ε2 − ∆2 8 Chapter 2. Quantum Coherent Superconducting Devices
In thermal equilibrium, we use Fermi-Dirac distribution function to find the density of quasi- particle excitations relative to Cooper pair density, s Z ∞ eq 2 eq 2πT −∆/T xqp = Nqp(ε)f (ε)dε = e . (2.16) N0∆ ∆ ∆ This predicts that the quasiparticle density can be arbitrarily suppressed by reducing the tem- perature and is essentially negligible if T ∆.
2.1.3 Superconducting gap In order to find the superconducting energy gap, we use the Bogoliubov transformation and find from Eq. (2.5),
X ∗ h † † i ∆ = λ uk0 vk0 hγ−k0↓γ−k0↓i − hγ−k0↓γ−k0↓i . (2.17) k0 As Bogoliubov quasiparticles are fermionic excitations, their occupation probability is deter- mined by the usual Fermi-Dirac distribution so that we can write,
† † hγ−k0↓γ−k0↓i − hγ−k0↓γ−k0↓i = 1 − 2n(εk0 ) = tanh(εk0 /2T ). (2.18)
We now change the summation in Eq. (2.17) to integration and given the coherence factors, Eqs. (2.8), we find
Z ωD q 1 1 2 2 = p tanh( ξ + ∆ /2T )dξ. (2.19) λN0 0 ξ2 + ∆2
In the limit where temperature approaches zero, we find 1 = sinh−1 ωD that in the weak λN0 ∆ coupling limit, λN0 1, results in,
−1/λN0 ∆ ' 2ωDe ωD. (2.20)
This indicates that the superconducting order parameter cannot be derived by treating the coupling strength in a perturbative way. We can alternatively express Eq. (2.19) in terms of quasiparticle energy, √ ω2 −∆2 1 Z D 1 = √ tanh(ε/2T )dε. (2.21) λN0 ∆ ε2 − ∆2 We use this relation in chapter5 to calculate the order parameter for a proximitized supercon- ductor relative to a bulk superconductor.
2.2 Josephson Effect
In a Josephson junction formed by two superconductors that are interrupted by an insulating tunnel barrier, a supercurrent flows through the device even in absence of an external bias. This phenomenon is due to the phase difference between the two superconducting electrodes forming the junction. In this subsection, we use perturbation theory to study the Josephson effect. Chapter 2. Quantum Coherent Superconducting Devices 9
One can alternatively use quasiparticle bound states to find the same Josephson equations [15]. Let us consider the following Hamiltonian that expresses single electron tunneling across the junction, ˜X X †R L HT = t (cmσcnσ + H.c.), (2.22) m,n σ where we assumed a constant tunneling matrix element t˜, and R and L are labeling right and left sides of the junction, respectively. The electron tunneling operator in terms of Bogoliubov operators reads,
†R L R† L R† L i(φR−φL) cmσcnσ =umunγnσ γmσ + vmvnγnσ γmσe R† L† R† L† iφR R L R L −iφL + unvm γn↓ γm↑ − γn↑ γm↓ e + umvn γn↑γm↓ − γn↓γm↑ e , (2.23) where the coherence factors, u and v, are taken real as Eq. (2.23) explicitly accounts for the phase difference across the junction. As the ground state of a superconductor is vacuum state for quasiparticles, the tunneling Hamiltonian results in zero expectation value. However, the tunneling Hamiltonian taken to the second-order perturbation theory given by, 1 H(2) = X H H , (2.24) T T ε T i i has a finite value in the ground states. Here εi is the energy of intermediate states and the Hamiltonian has terms that transfer two electrons to the right, two to the left, and with no net electron transfer. The latter leads to a constant value in the expectation value that has no physical effect. The terms with a net transfer to the right gives,
R R† L L† R R† L L† (2) γn↑γn↑ γm↓γm↓ + γn↓γn↓ γm↑γm↑ ˜2 X L R −iφL iφR L R hHT i = − t hΨBCS, ΨBCS|umvne L R unvme |ΨBCS, ΨBCSi n,m εn + εm 1 ˜2 i(φR−φL) X = − 2t e umvnunvm L R n,m εn + εm Z ∞ Z ∞ ∆ ∆ 1 ˜2 i(φR−φL) L R L R = − 2t e N0 N0 dξ dξ L R L R −∞ −∞ ε ε ε + ε Z ∞ Z ∞ 1 ˜2 i(φR−φL) L R = − 2t ∆e N0 N0 dθL dθR −∞ −∞ cosh θL + cosh θR 1 gT = − ∆ei(φR−φL), (2.25) 16 gK where the final equality is obtained by change of variables u = (θL + θR)/2 and v = (θL − θR)/2 in the last integral that, up to prefactors, results in complete elliptic integral of the first kind L R at zero, K(0) = π/2. Here, N0 and N0 are the density of states per spin at the Fermi energy 2 L R˜2 2 of the left and right electrode, gT = 4πe N0 N0 t is the junction conductance, gK = e /2π is the conductance quantum and we assumed equal gap for both sides of the junction. A similar calculation for the net transfer to the left gives the complex conjugate of Eq. (2.25). The sum of these two terms give the energy gain by electron pair tunneling,
U = −EJ cos φ. (2.26) 10 Chapter 2. Quantum Coherent Superconducting Devices
1 gT for which EJ = ∆ is Josephson energy and φ = φR − φL is the phase difference across the 8 gK junction. This energy is associated with a supercurrent that is driven by the phase difference across the junction which reads,
2π ∂U π ∆ IJ = = gT sin φ, (2.27) Φ0 ∂φ 2 e where Φ0 = h/2e denotes the superconducting flux quantum. As we have just shown, this supercurrent solely originates from the phase difference between the ground states of the super- conducting leads; therefore, it is a dissipationless current. If an external voltage V is imposed to the junction, the phase difference evolves in time according to the AC-Josephson effect,
Φ dφ V = 0 (2.28) 2π dt
Hence, the inductance of the Josephson junction LJ has a non-linear relation with the phase difference,
dIJ 1 1 LJ = V/ = (2.29) dt π∆gT cos φ
This non-linearity together with the ultra-low dissipation provided by superconductivity makes Josephson junctions promising candidates to build qubits.
2.3 Superconducting Qubits
In the previous section we ignored the fact that as supercurrent flows through the junction, charges build up on the islands of the junction and consequently Coulomb interactions become important as well. These repulsive interactions give another energy scale for the system that 2 is the charging energy. For a single electron transfer to the island, it becomes Ec = e /2C where C is the total capacitance that the island makes with the environment. Once we take into account the charging energy, it becomes clear the Josephson-junction-based devices can act like an artificial atom, as we explain in the following. In this section, we consider two types of superconducting qubits: Cooper-pair box and transmon qubit. The former is the earliest type of the superconducting qubits while the latter was realized some years later and is one the most promising ones in terms of the coherence time and scalability. In writing the Hamiltonian of the qubit, for the moment we neglect the presence of quasiparticles while just trying to give a brief qualitative explanation of their effect. In the next section we explicitly consider quasiparticles and study how their tunneling across the junction result in qubit energy relaxation.
2.3.1 Cooper-pair box
In a pioneering experiment by Nakamura and co-authors [16], the Cooper-pair box was the first superconducting device used to demonstrate quantum Rabi oscillations. As schematically depicted in the left panel of figure (2.1), this qubit simply consists of a Josephson junction where one of its islands is used to store charges, and the other island is to provide these charges. There is also a gate electrode enabling to shift the electrostatic potential of the island with respect to Chapter 2. Quantum Coherent Superconducting Devices 11
Figure 2.1: Left panel: Circuit diagram of a Cooper-pair box. The island is isolated by the Josephson junction and a capacitor. Tuning the gate voltage Vg enables us to control the number of extra Cooper pairs on the island. This voltage is sensitive to fluctuations in the charges that are surrounding the island. Right panel: Circuit diagram of a single-junction transmon qubit. The Junction is shunted by a large capacitance to increase the ration of EJ /EC that makes the qubit robust against the charge noise. the bulk electrode in order to tune the number of charges on the island. The Hamiltonian of this qubit reads, 2 ˆ H = Ec(Nˆ − Ng) − EJ cos φ, (2.30) where the operator Nˆ counts the number of single electrons that are tunneling-in or out of the island,
Nˆ|Ni = N|Ni, (2.31) and is conjugate to the Josephson phase operator, [φ,ˆ N/ˆ 2] = i. The offset charge Ng = CgVg/e is a continuous variable expressing the polarization charge on the island induced by the gate voltage Vg. Important feature of this qubit is that the island is made small enough such that the acces- sible thermal energy at millikelvin temperatures (where the qubit is operating) is much smaller than the charging energy, Ec kBT . Moreover, the charging energy also dominates the Joseph- son energy, Ec EJ ; in this condition, the number of extra charges on the island becomes a well defined variable. The qubit Hamiltonian in charge basis reads, E H = X E (Nˆ − N )2|NihN| − J (|NihN + 2| + |N + 2ihN|) . (2.32) c g 2 N
The charging energy as a function of the offset charge, Ng, gives a set of parabolas associated with single-electron charges, N, present at the island. On the other hand, the Josephson energy connects the nearby charge states with the same parity. By tuning the gate voltage such that the offset charge is close to the values where these parabolas cross each other, only the two crossing states remain important and the effective Hamiltonian become a 2 × 2 matrix. In particular, assuming initially there is no single electron present at the island, the qubit working point is at Ng = 1 that results in qubit states to be in superposition of |Ni and |N + 2i charge states, |Ni + |N + 2i |Ni − |N + 2i |0i = , and |1i = , (2.33) 2 2 12 Chapter 2. Quantum Coherent Superconducting Devices
Figure 2.2: Energy diagram for the first two low-laying states with even and odd parity. The zero point energy in each panel is chosen at the bottom of ground state and energies are ¯ 1 even odd normalized to average energy E1 = 2 (E0 + E1 ). Panel (a): for Cooper-pair box, where EJ /EC 1, energy levels have high charge dispersion. The dashed line in the figure points the offset charge at the qubit working point; any deviation from this point changes the qubit frequency. In addition, a transition from even to odd charge states induced by a single electron tunneling destroys the qubit state that is a superposition of charge states with same parity. Panels (b) and (c): as the ratio of EJ /EC is increased, the energy levels become less sensitive to the offset charge. Panel (d): in the transmon regime, where EJ /EC 1, the energy levels become insensitive to the offset charge. Moreover, the even and odd charge sectors contribute equally to the qubit logical state.
while the qubit frequency becomes ω01 = EJ . Panel (a) of figure (2.2) illustrates the eigenenergies of the two low-lying states for the even and odd sectors of the qubit Hamiltonian, Eq. (2.32). The figure makes it clear that the Cooper- pair box sufferers from two major drawbacks that limit its coherence time: First, the high charge dispersion of the energy levels makes qubit vulnerable to the charge noise. Indeed, fluctuations of the charges in the surrounding environment causes the offset charge deviate from the working point; this in turn modulates the qubit frequency and leads to qubit dephasing. Second, as it is shown in Eq. (2.33), the qubit states consist of symmetric superposition of charge states with Chapter 2. Quantum Coherent Superconducting Devices 13 equal parity; if an unpaired electron tunnels to the island, it poisons the device by changing the charge parity that brings the qubit out of its computational subspace. These issues limited coherence time of the Cooper-pair box to about 10−9 s that is more than 5 orders of magnitude less than the nowadays state-of-the-art qubits, [6, 17].
2.3.2 Transmon qubit Reducing the charge dispersion makes the qubit frequency less sensitive to the charge noise. This can be achieved, for example, by going into the transmon regime where the ratio of EJ /EC is large. In this case, the quantum fluctuations of the phase is relatively small, while the uncertainty of charge in the qubit state is significant. Going to this regime, however, would also reduce the anharmonicity of the energy levels, but remarkably, while the charge dispersion decreases exponentially in EJ /EC , the anharmonicity is suppressed algebraically with a slow power law in EJ /EC [4]. Indeed, there is a range for the ratio between energy scales EJ and EC for which the charge dispersion is flattened, rendering the qubit robust against charge noise, while enough anharmonicity is kept, thus avoiding the excitation of higher-level states. As schematically illustrated in right panel of figure 2.1, in order to realize the high ratio of EJ /EC , the junction in transmon qubit is shunted by a large capacitor; this lowers the total charging energy, EC , and makes EJ /EC large. Panel(d) of figure (2.2) shows the two energy levels of transmon qubit with even and odd charge parity, while panels (b) and (c) illustrate the crossover from charging regime to transmon regime. The picture clearly shows that as the ratio of EJ /EC is increased, the total charge dispersion rapidly decreases. Therefore, the energy difference between states with different parities rapidly decreases as well. Indeed, in transmon qubit the logical state of the qubit contains two physical states with even and odd parity while the energy difference between these two states is given by [4],
m m ω˜eo = ωeo cos(πNg), (2.34) where,
r 2m 2m+1 √ 2 2 8E 4 m p m J − 8EJ /EC ωeo = 4 8EC EJ (−1) e , (2.35) π m! EC for which m = 0 for the logical qubit ground state and m = 1 for the excited states. Therefore, a transmon qubit is also less disturbed from single-electron tunneling because this does not bring the qubit out of its computational subspace.
2.4 Qubit-Quasiparticle Interaction
So far, in writing down the qubit Hamiltonian, we have neglected the quasiparticle excitations. This is because, as it follows from Eq. (2.16), the quasiparticle density is negligible at millikelvin temperatures where superconducting qubits are operating. However, a number of experiments firmly confirm that at low temperatures (below around 0.1 Tc for Aluminum) quasiparticles fail to equilibrate with the environment and their density significantly exceeds the expected equilibrium value [18, 19, 20, 21]. These excitations have a detrimental effect on the performance of superconducting devices in a wide range of applications. To name a few, they limit the sensitivity of photon detectors in astronomy [22, 23] and cooling power of micro-refrigerators 14 Chapter 2. Quantum Coherent Superconducting Devices
[24, 25] and cause braiding errors in proposed Majorana-based quantum computation [26, 27, 28, 29]. In superconducting qubits, it has been firmly established both theoretically [30, 31, 32, 33, 34, 35, 36, 37] and experimentally [18, 38, 39, 40, 41, 42, 43, 44, 45] that quasiparticle tunneling causes qubit energy decay and dephasing. In addition, the residual nonequlibrium quasiparticles result in qubit excited state population in excess of thermal equilibrium value [46]. In this section we take quasiparticles into account and discus how their tunneling across the junction results in qubit energy relaxation. The system Hamiltonian in presence of quasiparticles can be divided into three parts
H = Hq + Hqp + Hint, (2.36) where Hq is the qubit Hamiltonian and the second term describes presence of quasiparticles on the left and right superconducting leads,
X X s† s Hqp = Ekγk,σγk,σ. (2.37) s=L,R k,σ
The third term describes quasiparticle tunneling across the junction; from Eq. (2.23) we write,
0 p Hint = Hint + Hint, (2.38)
0 p where the single quasiparticle tunneling, Hint, and pair tunneling, Hint, read up to a global phase factor,
0 ˜ X L R iφ/ˆ 2 R L −iφ/ˆ 2 L† R Hint =t (uk uk0 e − vk0 vk e )γkσ γk0σ + H.c., (2.39) k,k0,σ p ˜X L R iφ/ˆ 2 R L −iφ/ˆ 2 L† R† Hint =t [(uk vk0 e + uk0 vk e )γk↑ γk0↓ k,k0 R L −iφ/ˆ 2 L R iφ/ˆ 2 R L +(vk0 uk e + vk uk0 e )γk0↓γk↑] + (L ↔ R). (2.40)
2.4.1 Energy relaxation induced by quasiparticle tunneling The tunneling Hamiltonian makes possible qubit state transition occurring by exchanging energy with the tunneling quasiparticle. Up to lowest order in tunneling amplitude t˜, the transition from excited state, |1i, to the ground state, |0i, with qubit frequency ω10 is found using Fermi’s golden rule
X 2 Γ10 = 2πhh |h0, {λ}qp|Hint|1, {η}qpi| δ(Eλ,qp − Eη,qp − ω10)ii, (2.41) {λ}qp where {η}qp ({λ}qp) is the initial (final) state of quasiparticles with energy Eη,qp (Eλ,qp). The double angular brackets hh...iidenote averaging over initial quasiparticle states and the summa- tion is over all quasiparticle states. To calculate this rate, we note that the pair tunneling part p of the interaction Hamiltonian, Hint, contains terms creating or annihilating two quasiparticles; this absorbs or releases energy by amount twice the superconducting gap. On the other hand, superconducting qubits are designed such that the qubit frequency is much smaller than twice the gap, ωif 2∆, since this is necessary to avoid breaking Cooper pairs during qubit operation. Therefore, up to the leading order given by Fermi’s golden rule, Chapter 2. Quantum Coherent Superconducting Devices 15 the pair tunneling part does not contribute in the transition rate due to energy conservation. Moreover, we assume low-temperature limit so that the characteristic energy of quasiparticles, δε, (that is proportional to temperature and is measured from the gap) is small compared with superconducting energy gap,√δε ∆. This enables us to approximate the coherence factors, Eqs. (2.8), by uk ' vk0 ' 1/ 2 that in turn simplifies the single-quasiparticle tunneling to, ˆ 0 X φ L† R H = t˜ i sin γ γ 0 + H.c. (2.42) int 2 kσ k σ k,k0,σ
The transition rate then factorizes into terms that separately account for qubit dynamics and quasiparticle kinetics,
φˆ Γ = |h0| sin |1i|2S (ω ) (2.43) 10 2 qp 10 where the quasiparticle current spectral density becomes,
˜2 X X L† R R† L 2 Sqp(ω) =2πt hh |h{λ}qp|γkσ γk0σ + γk0σγkσ|{η}qpi| δ(Eλ,qp − Eη,qp − ω)ii 0 k,k ,σ {λ}qp ˜2 X R† R L L† =4πt hhh{η}qp|γk0σγk0σ|{η}qpih{η}qp|γkσγkσ |{η}qpiδ(Eλ,qp − Eη,qp − ω)ii k,k0,σ 32E Z ∞ = J n(ε)n(ε + ω)f(ε)[1 − f(ε + ω)]dε. (2.44) π∆ ∆
R† R R L L† L Here we used hhh{η}qp|γ γ |{η}qpiii = f(ε ), hhh{η}qp|γ γ |{η}qpiii = 1 − f(ε ) and took L R L R L R Eλ,qp − Eη,qp = Eλ,qp + Eλ,qp − Eη,qp − Eη,qp = ε − ε . The spectral density depends on the quasiparticle distribution function; assuming “cold” quasiparticles meaning their energy (or effective temperature) is small compared with qubit q ∆ frequency, δε ω, we can take 1 − f(ε + ω) ' 1 and n(ε + ω) = 2ω . This simplifies the spectral function and we find; s 8E 2∆ S (ω) = J x (2.45) qp π qp ω where 2 Z ∞ xqp = n(ε)f(ε)dε, (2.46) ∆ ∆ is the density of quasiparticles normalized to the Cooper-pair density. In thermal equilibrium this quantity is given by Eq. (2.16). However, we note that Eq. (2.45) is valid for arbitrary distribution function provided δε is the smallest energy scale of the system. To find the qubit excitation rate, Γ01, one has to calculate Sqp(ω) for ω < 0 that is obtained from Eq. (2.44) by replacing ε → ε − ω, ω → −ω; within our low-temperature assumption, in general we have S(−ω) S(ω) indicating that there is no quasiparticle with energy high enough to excite the qubit. Eq. (2.45) is of central importance in this thesis as it indicates that the qubit decay rate can be decreased by reducing the quasiparticle density near the Josephson junction. 16 Chapter 2. Quantum Coherent Superconducting Devices
Figure 2.3: Panel (a) is reproduced and slightly modified from Ref. [19]. It illustrates experimen- tal data points for the number of quasiparticles in a superconducting resonator and compares experimental findings with theoretical prediction in thermal equilibrium. At low temperatures relevant to the operation of superconducting qubits, the residual quasiparticle density is sig- nificantly higher than theory predictions. Panel (b) is reproduced from Ref. [43] and makes it clear that suppressing the quasiparticle density, that in this case is achieved in a 3D transmon qubit by cooling in magnetic field to generate vortices, can improve the qubit coherence times. It is difficult to control vortices and it is observed that a large number of them can negatively influence qubit performance.
In figure (2.3) we have shown some experimental highlights about quasiparticles and their impact on the qubit decay rate. Panel (a) shows the measured quasiparticle density in a su- perconducting resonator as a function of temperature and reveals that the density saturates when temperature goes below ∼ 160 mK. While a detailed knowledge of the source that gener- ates nonequilibrium quasiparticles is eventually needed to solve quasiparticle-related problems, physicists have been looking for ways to suppress quasiparticle density that promises improving the qubit coherence. One proposal that has been recently realized is to cool down the qubit in a magnetic field that would generate vortices in the device. At the core of a vortex, super- conducting order parameter is suppressed, which makes it possible to trap quasiparticles- we Chapter 2. Quantum Coherent Superconducting Devices 17 describe the trapping mechanism in detail in the next chapter. Panel (b) of figure 2.3 shows measured relaxation times for a 3D transmon qubit as a function of magnetic field. The strength of magnetic field determines the number of generated vortices and, consequently, the level of suppression in the quasiparticle density. The plot makes it clear that up to some point in the magnetic field, vortices could improve the qubit coherence while for magnetic field above ' 200 mG, the qubit performance is negatively affected. This behavior is attributed to the energy dissipation that a large number of vortices can cause [43]. Moreover, it is difficult to control the vortex position. This has motivated us to study another method for suppressing quasipar- ticle density that enables us to control trap size and placement. In the next chapter, we will introduce normal-metal quasiparticle traps and discuss how they work and how they can be optimized by proper trap placement.
Chapter 3 Normal-Metal Quasiparticle Traps
We begin this chapter by explaining how a normal-metal connected to superconducting qubit can act as a sink for quasiparticles. Section 3.1 contains part of our work that has been published under the title Normal-metal quasiparticle traps for superconducting qubits and cited in Ref. [47]. Here we develop a model for the effect of a single small trap on the dynamics of the excess quasiparticles injected in a transmon-type qubit. Section 3.2 containes a paper of the author that has been published with title Optimal configurations for normal-metal traps in transmon qubits and cited in Ref. [68]. Here we build on section 3.1 and discuss how quasiparticle trapping can be optimized. We show proper trap design can increase the slowest decay rate of quasiparticle and at the same time suppress quasiparticle steady-state density and its fluctuations. I co-authored Ref. [47] and contributed by discussing the model and experimental data, comparing simplified analytical results with exact numerics and preparing a number of sug- gested figures for the paper. I contributed to Ref. [68] by doing all of exact and numerical modelings and their corresponding figures to demonstrate enhancing the decay rate of the ex- cess quasiparticle density as well as suppression of the quasiparticle steady-state density due to normal-metal traps, comparing normal-metal traps on the pads with vortex trapping, and the analysis of traps for Xmon qubits.
3.1 Modeling
3.1.1 Introduction Ideal superconducting devices rely on dissipationless tunneling of Cooper pairs across a Joseph- son junction. For example, in a Cooper pair pump [48], the controlled transport of Cooper pairs across two or more junctions can in principle make it possible to relate frequency and current and hence enable metrological applications of such a device [49]. For quantum information purposes, the non-linear relation between the supercurrent and the phase difference across a junction makes the junction an ideal non-linear element to build a qubit [50]. However, in addi- tion to the pairs tunneling, single-particle excitations known as quasiparticles can also tunnel. In the pumps this leads to “counting errors”, limiting the accuracy of the current-frequency re- lation [48, 49]. In qubits, quasiparticles interact with the phase degree of freedom, providing an unwanted channel for the qubit energy relaxation [32,6]. While in many cases it is impossible to prevent the creation of quasiparticles, one may keep them away from the Josephson junctions by trapping. Evacuation of the quasiparticles from the vicinity of the junction provides a way to extend the energy relaxation time (T1) in the steady state, and to restore the steady state after a perturbation, whether caused by qubit operation or some uncontrolled environmental effect. Quasiparticle trapping has been explored for a long time, and various proposal exists on how to implement such a trapping. For example, gap engineering takes advantage of the fact 20 Chapter 3. Normal-Metal Quasiparticle Traps that quasiparticles accumulate in regions of lower gap to steer them into or away from certain parts of the device. Gap engineering was used successfully to limit quasiparticle “poisoning” in a Cooper pair transistor [51], while proved ineffective in a transmon qubit [52]. A vortex in a superconducting film can also act as a well-localized trap, since the gap is completely suppressed at the vortex position. Trapping by vortices has been demonstrated [53, 54, 43, 55], but vortex motion may induce an unwanted dissipation. An island of a normal metal in contact with the superconductor may also serve as a quasiparticles trap [25, 56]. In the limit of weak electron tunneling across the contact, the proximity effect is negligible. The quasiparticles tunneled into the normal metal are trapped there upon losing their energy by phonon emission or inelastic electron-electron scattering. The majority of previous works concentrated on the control of a steady-state quasiparticle population [49, 25, 56]. In contrast, we are interested in the effect of a normal-metal trap on the dynamics of the quasiparticle density. Traps accelerate the evacuation of the excess quasiparticles injected in a qubit in the process of its operation. Our main goal is to determine how the characteristic time of the evacuation depends on the parameters of a small normal- metal island in contact with the superconducting qubit. The characteristic time shortens with the increase of the trap size, saturating at a value dependent on the qubit geometry and the quasiparticle diffusion coefficient. The size at which a trap becomes effective depends on the contact resistance, the energy relaxation rate in the normal-metal island, and the effective temperature of the quasiparticles. We develop a simple model allowing to evaluate the time evolution of the quasiparticle density and find the characteristic evacuation time as a function of the trap parameters. The model is validated by measurements of the qubit T1 relaxation time performed on a series of transmons with normal-metal traps of various sizes. This section is organized as follows: in Sec. 3.1.2 we develop a phenomenological quasiparticle diffusion and trapping model which includes the effect of a normal-metal trap. In Sec. 3.1.3 we study the dynamics of the density during injection and trapping in a simple configuration, and in Sec. 3.1.4 we provide experimental data (obtained by our collaborators in Yale University) supporting our approach.
3.1.2 The diffusion and trapping model
Let us consider a quasiparticle trap made of a normal (N) metal covering part of a super- conducting (S) qubit. The contact between the two superconductor and the normal trap is provided by an insulating (I) layer characterized by a small electron transmission coefficient. In order to relate the quasiparticle tunneling rate to the conductance of the contact, we use the tunneling Hamiltonian formalism applied to a model N-I-S system, see Fig. 3.1,
H = Hqp + HN + HT , (3.1) X † Hqp = εnγnσγnσ , (3.2) nσ X † HN = ξmcmσcmσ , (3.3) mσ
te X † † HT = √ cmσdnσ + dnσcmσ . (3.4) ΩN ΩS m,n,σ Chapter 3. Normal-Metal Quasiparticle Traps 21
1 N dN ⌫N r
0 0 1 2 ✏ 3