arxiv: v1 [cond-mat.mes-hall] 3 Mar 2016

Transcription

1 Quantum Phase Slip Noise Andrew G. Semeno 1,3 and Andrei D. Zaikin,1 1 I.E.Tamm Department of Theoretical Physics, P.N.Lebede Physical Institute, Moscow, Russia Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT, 761 Karlsruhe, Germany 3 National Research Uniersity Higher School of Economics, 11 Moscow, Russia arxi: [cond-mat.mes-hall] 3 Mar 16 Quantum phase slips (QPS generate oltage fluctuations in superconducting nanowires. Employing Keldysh technique and making use of the phase-charge duality arguments we deelop a theory of QPS-induced oltage noise in such nanowires. We demonstrate that quantum tunneling of the magnetic flux quanta across the wire yields quantum shot noise which obeys Poisson statistics and is characterized by a power law dependence of its spectrum S on the external bias. In long wires S decreases with increasing frequency and anishes beyond a threshold alue of at T. Quantum coherent nature of QPS noise yields non-monotonous dependence of S on T at small. PACS numbers: 73.3.Ra, 74.5.F-, n Keywords: phase slips, shot noise, Can a superconductor generate oltage fluctuations? More specifically, if an external bias is applied to a superconductor could the latter produce shot noise? Posing these questions we, of course, imply that temperature T, characteristic frequencies and/or oltages as well as all other releant energy parameters remain well below the superconducting gap, i.e. the superconductor is either in or sufficiently close to its quantum ground state. At the first sight, positie answers to both these questions can be rejected on fundamental grounds. Indeed, a superconducting state is characterized by zero resistance, i.e. a non-dissipatie current below some critical alue can pass through the system. Hence, neither non-zero aerage oltage nor oltage fluctuations can be expected. These simple considerations although applicable to bulk superconductors become insufficient in the case of ultrathin superconducting wires because of the presence of quantum phase slips (QPS [1 4]. In such wires quantum fluctuations of the superconducting order parameter field = e iϕ play an important role being responsible for temporal local suppression of inside the wire and, hence, for the phase slippage process. Each quantum phase slip eent corresponds to the net phase jump by δϕ = ±π implying positie or negatie oltage pulse δv = ϕ/e (here and below we set = 1 and tunneling of one magnetic flux quantum Φ π/e = δv(t across the wire in the direction perpendicular to its axis. Biasing the wire by an external current I one breaks the symmetry between positie and negatie oltage pulses making the former more likely than the latter. As a result, the net oltage drop V occurs across the wire also implying non-zero resistance R = V/I which may not anish down to lowest T [5, 6], as it was indeed obsered in a number of experiments [7 9]. Hence, in the presence of QPS the current flow becomes dissipatie and according, e.g., to the fluctuation-dissipation theorem (FDT one should also expect oltage fluctuations to occur in the system. While these arguments suggest a positie answer to the first of the aboe questions they do not yet specifically address shot noise. Two key pre-requisits of shot noise are: (i the presence of discrete charge carriers (e.g., electrons in the system and (ii scattering of such carriers at disorder. Although discrete charge carriers Cooper pairs are certainly present in superconducting nanowires, they form a superconducting condensate flowing along the wire without any scattering. For this reason the possibility for shot noise to occur in superconducting nanowires appears by no means obious. In this Letter we will perform a detailed theoretical analysis of QPS-induced oltage fluctuations in ultrathin superconducting wires. In particular, we will demonstrate that quantum tunneling of magnetic flux quanta Φ across the wire causes shot noise which obeys Poisson statistics and shows a non-triial dependence on temperature, frequency and external current. The model and effectie Hamiltonian. The system under consideration is displayed in Fig. 1. It consists of an ultrathin superconducting wire of length L and cross section s and a capacitance C switched in parallel to this wire. The right end of the wire (x = L is grounded as shown in the figure (x is the coordinate along the wire ranging from to L. The oltage V(t at its left end x = fluctuates and such fluctuations can be measured by a detector. The whole system is biased by an external current I = V x /R x. I V(t Rx C Phase Slip QPS wire Vx superconductor FIG. 1: (Color online The system under consideration. The figure also illustrates creation of two plasmons by a QPS. An effectie Hamiltonian for our system can be written in the form Ĥ = Ĥ Ch Iϕ/e+Ĥ wire. (1

2 The first and the second terms in the right-hand side of Eq. (1 account respectiely for the charging energy [1] Ĥ Ch = 1 ( i C (ϕ/e +Q ( and for the potential energy tilt produced by an external current I. The ariable ϕ(t ϕ(,t represents the phase of the superconducting order parameter field (x,t at x =. Here we also set ϕ(l,t. The last term Ĥwire in Eq. (1 describes the superconducting wire. This part of the effectie Hamiltonian can be expressed in terms of both the modulus (x,t and the phase ϕ(x,t of the order parameter field [5, 6, 11]. Here, howeer, we will proceed differently and employ the duality arguments. The duality between the phase and the charge ariables was established and discussed in details in the case of ultrasmall Josephson junctions [1, 1 14]. Later the same duality arguments were extended to short [15] and long [16 18] superconducting wires. According to the results [18] the dual representation for the Hamiltonian of a superconducting nanowire is defined by an effectie sine-gordon model Ĥ wire = ĤTL +ĤQPS. (3 In the absence of quantum phase slips such nanowire can be described as a transmission line with ( Φ Ĥ TL = + ( xχ L kin C w Φ, (4 where L kin = 1/(πσ N s and C w are respectiely the kinetic wire inductance (times length and the geometric wire capacitance (per length, [Φ(x,χ(x ] = iφ δ(x x (5 defines the commutation relation between the canonically conjugate flux (or phase and charge operators, σ N is the normal state Drude conductance of the wire and is the superconducting gap. The term Ĥ QPS = γ QPS cosχ (6 accounts for the effect of quantum phase slips and γ QPS (g ξ /ξexp( ag ξ, a 1, (7 is the QPS tunneling amplitude [6] per unit wire length with g ξ = πσ N s/(e ξ 1 being the dimensionless normal state conductance of the wire segment of length equal to the coherence length ξ. The physical meaning of the quantum field χ(x,t is transparent: It is proportional to the total charge q(x, t that has passed through the point x up to the time moment t, i.e. q(x,t = χ(x,t/φ. Accordingly, the local current I(x, t and the local charge density ρ(x, t are defined as I(x,t = t χ(x,t/φ, ρ(x,t = x χ(x,t/φ, (8 thereby satisfying the continuity equation. The charge Q in Eq. ( equals to Q(t = χ(,t/φ. Keldysh technique and perturbation theory. In order to proceed we will make use of the Keldysh path integral technique. Accordingly, our ariables of interest need to be defined on the forward and backward time branches of the Keldysh contour, i.e. we now hae ϕ F,B (t and χ F,B (x,t. As usually, it is conenient to also introduce the classical and quantum ariables, respectiely ϕ + (t = (ϕ F (t + ϕ B (t/ and ϕ (t = ϕ F (t ϕ B (t (and similarly for the χ-fields. Making use of the Josephson relation between the oltage and the phase one can formally express the expectation alue of the oltage operator across the the superconducting wire in the form where and V(t 1 = 1 e S QPS = γ QPS ϕ+ (t 1 e isqps (9 sin(χ + sin(χ / (1... = D ϕ(td χ(x,t(...e is[φ,χ] (11 implies aeraging with the Keldysh effectie action S corresponding to the Hamiltonian Ĥ = Ĥ Ĥ QPS. Analogously, for the oltage-oltage correlator V(t 1 V(t = 1 {ˆV(t 1, ˆV(t } (where curly brackets denote the anticommutator one has V(t 1 V(t = 1 4e ϕ+ (t 1 ϕ + (t e isqps, (1 Higher oltage correlators are defined similarly. Their analysis, howeer, is beyond the frames of this work. Eqs. (9 and (1 are formally exact expressions which we are now going to ealuate. To this end we will employ the regularperturbationtheory in γ QPS (7 which can be regarded as a small parameter of our theory. In the zero order in γ QPS the problem is described by the quadratic (in both ϕ and χ Hamiltonian Ĥ and all aerages can be handled exactly with the aid of the Green functions G R ab (X,X = i a + (Xb (X, G K ab(x,x = i a + (Xb + (X, (13 where a(x and b(x stand for one of the fields ϕ(t and χ(x,t. As both these fields are real, the adanced and retarded Green functions obey the condition G A ab (ω = G R ba ( ω. With this in mind the Keldysh function GK can be expressed in the form G K ab (ω = 1 coth ( ω (G R ab (ω G R ba ( ω. (14

3 3 Expanding Eqs. (9 and (1 up to the second order in γ QPS and performing all necessary aerages we ealuate the results in terms of the Green functions (13, see Supplemental material for further details. The result of our calculation both for the aerage oltage (9 and for the oltage-oltage correlator (1 can also be expressed in the form of candy diagrams displayed in Fig.. Comparing the result (17 with that found in Ref. 5 we immediately conclude that Γ QPS (I defines the quantum decay rate of the current state due to QPS. In [5] this rate was ealuated from the imaginary part of the free energyγ QPS (I = ImF. Herewe deriehe expression for Γ QPS by means of the real time technique without employing the ImF-method. Making use of the aboe results, ealuating the Green functions (13 (see Supplemental material and keeping in mind the detailed balance condition we obtain P x,x (ω = e ω T Px,x ( ω (19 V = Φ LγQPS ( ς Φ I sinh ( Φ I, ( where = 1/ L kin C w is the plasmon elocity [19], ς(ω = τ λ (πt λ 1Γ( λ ( iω πt Γ λ + iω πt, (1 Γ(λ FIG. : Candy-like diagrams which determine both aerage oltage (9 (upper diagram and oltage noise (1 (six remaining diagrams in the second order in γ QPS. The fields ϕ +, χ + and χ in the propagators (13 are denoted respectiely by way, solid and dashed lines. I-V cure and oltage noise. To begin with, let us briefly re-derie the results [5] for the aerage oltage within the framework of our technique. We obtain V = iγ QPS 4e ( lim ω ωg R ϕχ (x;ω (P x,x ( Φ I P x,x (Φ I, (15 where P x,x (ω = P x,x (ω+ P x,x (ω and P x,x (ω = e iωt e ig(x,x ;t,, (16 G(x,x ;t, = G K χχ (x,x ;t, 1 GK χχ (x,x;t,t 1 GK χχ (x,x ;,+ 1 GR χχ (x,x ;t,. Bearingin mind that lim ω ωg R ϕχ(x;ω = πi, Eq. (15 can be cast to the form V = Φ (Γ QPS (I Γ QPS ( I, (17 where we identify Γ QPS as Γ QPS (I = γ QPS P x,x (Φ I. (18 τ 1/ is the QPS core size in time and Γ(x is the Gamma-function. Here we also introduced the parameter [5] λ = R q /Z w s, where R q = π/e is the superconducting quantum resistance unit and Z w = L kin /C w is the wire impedance. It is satisfactory to obsere that the result (, (1 matches with that found in Ref. 5 by means of a different technique []. Let us now turn to oltage fluctuations. Our perturbatie analysis allows to recoer three different contributions to the noise power spectrum, i.e. S = e it V(tV( = S ( +Sr +Sa. ( The first of these contributions S ( has nothing to do with QPS and just defines equilibrium oltage noise for a transmission line. It reads S ( = i coth ( ( G R 16e ϕϕ ( G R ϕϕ (. (3 The other two terms are due to QPS effects. The term S r is also proportional to coth( and contains the products of two retarded (adanced Green functions: S r = γ QPS coth ( 8e Re [ G R ϕχ(x; (F x,x (G R ϕχ(x ; F x,x (G R ϕχ(x; ]. (4 Here we denoted F x,x ( = P x,x (+Φ I P x,x ( Φ I + P x,x ( +Φ I+ P x,x ( Φ I. (5

4 4 The remaining term S a, in contrast, contains the product of one retarded and one adanced Green functions and scales with the combinations C ± = coth ( ±Φ I coth ( as S a = γ QPS 16e G R ϕχ(x;g R ϕχ(x ; (6 [ ] C ± (P x,x (±Φ I P x,x ( Φ I. ± Eqs. (-(6 together with the expressions for the Green functions (see Supplemental material fully determine the oltage noise power spectrum of a superconducting nanowire in the perturbatie in QPS regime and represent the central result of this work. In the zero bias limit I the term S a anishes, and the equilibrium noise spectrum S = S ( + Sr is determined from FDT, see also [18]. At non-zero bias alues the QPS noise turns non-equilibrium. In the zero frequency limit the terms S ( and Sr tend to zero, and the oltage noise S S is determined solely by S a. Then from Eq. (6 we obtain S = Φ (Γ QPS (I+Γ QPS ( I = Φ coth ( Φ I V, (7 where V is specified in Eqs. (17, (. Combining the result (7 with Eqs. (, (1 we find S { T λ, T Φ I, I λ, T Φ I. (8 At higher temperatures T Φ I (though still T Eq. (8 just describes equilibrium oltage noise S = R of a linear Ohmic resistor R = V /I T λ 3 [5]. In the opposite low temperature limit T Φ I it accounts for QPS-induced shot noise S = Φ V obeying Poisson statistics with an effectie charge equal to the flux quantum Φ. This result sheds light on the physical origin of shot noise in superconducting nanowires: It is produced by coherent tunneling of magnetic flux quanta Φ across the wire. In the dual picture [18] such flux quanta can be iewed as charged quantum particles passing through (and being scattered at an effectie spatially extended tunnel barrier. Note that preiously the result analogous to Eq. (7 was deried for thermally actiated phase slips (TAPS [1]. This similarity appears remarkable gien a crucial physical difference between TAPS and QPS: The former can be regarded as classical (i.e. incoherent and noninteracting objects, whereas the latter are fully coherent [] forming an interacting quantum gas. Another interesting limiting case is that of sufficiently high frequencies and/or long wires /L. In this limit we obtain S ( = λ 8πe ( coth (/E C +(λ/π. (9 This contribution is independent of the wire length L. At low T and /λ E C = e /C we hae S ( 1/, i.e. the wire may generate 1/f oltage noise. Ealuating the QPS terms S r and Sa we obsere that the latter scales linearly with the wire length L whereas the former does not. Hence, the term S r can be safely neglected in the long wire limit. For the remaining QPS term S a we get S a = Lλ γ QPS 4e [ ( ( Φ I ς Φ I ς sinh ( Φ I ς ( ΦI ] + ((/E C +(λ/π sinh ( At T from Eq. (3 we find { S a I λ 1 (I /Φ λ 1, < Φ I/,, > Φ I/.. (3 (31 This result can be interpreted as follows. At T = each QPS eent excites (at least two plasmons [3] (see Fig. 1 with total energy E = Φ I propagating in the opposite directions along the wire. One plasmon (with energy E/ gets dissipated at the grounded end of the wire while another one (also with energy E/ reaches its opposite end causing oltage fluctuations (emits a photon with frequency measured by a detector. Clearly, at T = this process is only possible at < E/ in the agreement with Eq. (31. FIG. 3: (Color online The frequency dependence of the QPS noise spectrum S (3 at λ =.7, large E C and different T in the long wire limit. The inset shows S as a function of T. The result (3 is also illustrated in Fig. 3. At sufficiently small (we still keep /L one obseres a

5 5 non-monotonousdependence of S on T which is a direct consequence of quantum coherent nature of QPS noise. Finally, we point out that the perturbatie in γ QPS approach employed here is fully justified for not too thin wires with λ > λ c [5]. In wires with λ < λ c (characterized by unbound QPS-anti-QPS pairs γ QPS gets effectiely renormalized to higher alues and, hence, the perturbation theory eentually becomes obsolete. Howeer, een in this case our results may still remain applicable at sufficiently high temperature, frequency and/or current alues. In the low energy limit long wires with λ < λ c show an insulating behaior, as follows from the exact solution of the corresponding sine-gordon model [4]. This solution suggests that also oltage fluctuations become large in this limit. In summary, we demonstrated that quantum phase slips generate oltage noise in superconducting nanowires. In the presence of a current bias I quantum tunneling of the magnetic flux quanta Φ across the wire causes Poissonian shot noise with a non-triial power law dependence of its spectrum on both I and frequency. Our predictions can be directly erified in future experiments and need to be obsered while optimizing the operation of QPS qubits [5]. We acknowledge useful discussions with K.Yu. Arutyuno, D.S. Golube and P. Hakonen. [1] K.Yu. Arutyuno, D.S. Golube, and A.D. Zaikin, Phys. Rep. 464, 1 (8. [] A. Bezryadin, J. Phys.: Condens. Matter, 43 (8. [3] A.D. Zaikin, in: Handbook of Nanophysics: Nanotubes and Nanowires (CRC Press, Boca Raton, FL, 1 p [4] A. Bezryadin, Superconductiity in Nanowires (Wiley- VCH, Weinheim, 13. [5] A.D. Zaikin et al., Phys. Re. Lett. 78, 155 (1997. [6] D.S. Golube and A.D. Zaikin, Phys. Re. B 64, 1454 (1. [7] A. Bezryadin, C.N. Lau, and M. Tinkham, Nature 44, 971 (. [8] C.N. Lau et al., Phys. Re. Lett. 87, 173 (1. [9] M. Zgirski et al., Phys. Re. B 77, 5458 (8. [1] G. Schön and A.D. Zaikin, Phys. Rep. 198, 37 (199. [11] A. an Otterlo et al., Eur. Phys. J. B 1, 131 (1999. [1] S.V. Panyuko and A.D. Zaikin, J. Low Temp. Phys. 73, 1 (1988. [13] D.V. Aerin and A.A. Odintso, Phys. Lett. A 14, 51 (1989. [14] A.D. Zaikin, J. Low Temp. Phys. 8, 3 (199. [15] J.E. Mooij and Yu.V. Nazaro, Nat. Phys., 169 (6. [16] A.G. Semeno, PhD thesis (P.N. Lebede Physics Institute, Moscow, 1. [17] A.M. Hriscu and Yu.V. Nazaro, Phys. Re. B 83, (11. [18] A.G. Semeno and A.D. Zaikin, Phys. Re. B 88, 5455 (13. [19] J.E. MooijandG.Schön, Phys.Re.Lett.55, 114(1985. [] For the sake of simplicity here we assume that τ x, where x ξ is the QPS core size in space. [1] D.S. Golube and A.D. Zaikin, Phys. Re. B 78, 1445 (8. [] Experimental eidence for quantum coherent behaior of QPS was demonstrated by O.V. Astafie et al., Nature 484, 355 (1. [3] Due to momentum conseration plasmons can only be created in pairs with the total zero momentum. [4] See, e.g., D. Controzzi, F.H.L. Essler, and A.M. Tselik, Phys. Re. Lett. 86, 68 (1 and references therein. [5] J.E. Mooij and C.J.P.M. Harmans, New J. Phys. 7, 19 (5.

6 Supplemental Material 1 A. Perturbation theory Let us expand the general expressions (9 and (1 up to the second order in γ QPS. It is easy to demonstrate that linear in γ QPS terms anish identically in both expressions after aeragingoer the zero mode contained in the χ-field. In order to ealuate the terms γ QPS it is conenient to make a shift χ +(t Φ It+χ + (t and to decompose the aerages by means of the Wick theorem. As a result we obtain and V(t 1 = γ QPS e V(t 1 V(t = 1 4e ϕ +(t 1 ϕ + (t γ QPS e ϕ + (t 1 χ (x,t cos(φ I(t t +χ + (x,t χ + (x,t cos(χ (x,t/sin(χ (x,t / (S1 ϕ + (t 1 χ + (x,t ϕ + (t χ (x,t sin(φ I(t t+χ + (x,t χ + (x,tcos(χ (x,t/sin(χ (x,t / γ QPS e ϕ + (t 1 χ + (x,t ϕ + (t χ (x,t sin(φ I(t t+χ + (x,t χ + (x,tsin(χ (x,t/cos(χ (x,t / γ QPS e ϕ + (t 1 χ (x,t ϕ + (t χ + (x,t sin(φ I(t t+χ + (x,t χ + (x,tcos(χ (x,t/sin(χ (x,t / γ QPS e ϕ + (t 1 χ (x,t ϕ + (t χ + (x,t sin(φ I(t t +χ + (x,t χ + (x,t cos(χ (x,t/sin(χ (x,t / γ QPS 4e ϕ + (t 1 χ (x,t ϕ + (t χ (x,t cos(φ I(t t +χ + (x,t χ + (x,t cos(χ (x,t/cos(χ (x,t / (S The aerages in Eqs. (S1 and (S are Gaussian and, hence, can be handled in a straightforward manner. After that we immediately arrie at our final results for the I V cure (15, (16 and for the oltage noise spectrum (-(6. Both these results are expressed ia the function P x,x (ω (16 which in turn contains the Green function G(x,x ;t,. B. Analytic structure of the Green functions Let us define a more general Green function G χ (x,x ;σ which depends on the complex time σ and obeys the condition G χ (x,x ;t i = G(x,x ;t,. With the aid of the Kubo-Martin-Schwinger condition one can deduce that the Green function G χ is periodic in the imaginary time direction, i.e. G χ (x,x ;σ = G χ (x,x ;σ i/t. This function is analytic and has branch cuts at Im(σ = N/T for all integer N. The function exp(ig χ (x,x ;σ has the same analytic properties. One can write P x,x (ω = e iωt e igχ(x,x ;t i. (S3 (S4

7 Distorting the integration path and utilizing the property G χ (x,x ;σ = Ḡχ(x,x ; σ together with Eq. (S3 we arrie at the detailed balance condition (19. and C. Green functions The Green functions for the system displayed in Fig. 1 can be ealuated directly with the results: G R χχ (x,x ;ω = ( 4πλ cos G R 1 ϕϕ(ω = ω E C + iω 4e R x ωλ π cot(, (S5 ωl G R χϕ (x;ω = GR ϕχ (x;ω = ( ( ( ω(l x cos ωx iλcos ω E C + iω 4e R x θ(x x +cos ωsin ( ωl + ( ω(l x sin ( ωl ( ω(l x ωλ π cos( ωl cos ( ωx θ(x x ( ω(l x cos 4λ cos sin ( (( ωl ω E C + iω 4e R x sin ( ωl (S6 ( ω(l x ωλ π cos( ωl (S7 The last two expressions take a much simpler form in the long wire limit, in which case all plasmon excitations moing towards the grounded end of the wire eentually disappear and neer pop up again while excitations moing in the opposite direction produce oltage fluctuations measured by a detector. In this limit Eqs. (S6 and (S7 reduce to λe iωx G R ϕχ(x;ω ( (ω +i ω E C + iλ π, G R χχ(x,x ;ω πiλ ω x x ω +i ei. (S8 Here we also set R x as requested in the current bias limit. In order to ealuate the general expressions for the I V cure (15, (16 and for the oltage noise (-(6 it is necessary to compute the integral Υ(ω, = L/ L/ e i (x x P x,x (ω (S9 L/ L/ Separating the left moers and the right moers, making use of the explicit form of the Green function G R χχ(x,x ;ω and introducing the high frequency cutoff ω c 1/τ in order to aoid unphysical diergencies we obtain where Υ(ω, L (z = Performing the integration in Eq. (S11 we find e i x P x, (ω = L ( ω + ( ω, (S1 e izt sinh λ (πtτ sinh λ/ (πt(τ t+isinh λ/ (πt(τ +t i. (ω = λ (πtτ λ πt Γ ( λ iω πt ( Γ λ + iω ω πt e Γ(λ ς(ωe ω. The function ς(ω (1 is directly employedin our results both forthe I V cure( and for the QPS noise spectrum (7, (3. (S11 (S1