arxiv: v1 [cond-mat.mes-hall] 9 Jan 2019

Transcription

1 Fortschritte der Physik, January 1, 19 Quantum phase slips and voltage fluctuations in superconducting nanowires FQMT-15 Andrew G. Semenov 1,3 and Andrei D. Zaikin,1, arxiv: v1 [cond-mat.mes-hall] 9 Jan 19 We argue that quantum phase slips QPS may generate non-equilibrium voltage fluctuations in superconducting nanowires. In the low frequency limit we evaluate all cumulants of the voltage operator which obey Poisson statistics and show a power law dependence on the external bias. We specifically address quantum shot noise which power spectrum S Ω may depend nonmonotonously on temperature. In the long wire limit S Ω decreases with increasing frequency Ω and vanishes beyond a threshold value of Ω at T. Our predictions can be directly tested in future experiments with superconducting nanowires. 1 Introduction To conduct electric current without any resistance is the most fundamental property of any bulk superconducting material. Usually, the behavior of such superconductors is well described within the framework of the standard mean field theory. The situation changes dramatically as soon as a superconducting structure e.g., forming a narrow wire becomes sufficiently thin. In such structures thermal and/or quantum fluctuations start playing an important role being responsible for temporal local suppression of the superconducting order parameter = e i ϕ inside the wire and, hence, for the phase slippage process. This process gives rise to novel physical phenomena which cannot adequately be described with the aid of the mean field theory. In the low temperature limit thermal fluctuations are irrelevant and the system behavior is dominated by quantum phase slips QPS [1 ]. Each QPS accounts for the net phase jump by δϕ = ±π implying a voltage pulse δv = ϕ/e and tunneling of one magnetic flux quantum Φ π/e = δv t d t across the wire in the direction perpendicular to its axis. Formally such QPS events can be considered as quantum particles interacting logarithmically between each other in space-time [5]. Accordingly, the ground state of ultrathin superconducting wires can be described in terms of a d gas of interacting quantum phase slips with effective fugacity proportional to the QPS tunneling amplitude per unit wire length [6] γ QP S g ξ /ξexp ag ξ, a 1. 1 Here g ξ = πσ N s/e ξ 1 is the dimensionless normal state conductance of the wire segment of length equal to the coherence length ξ, is the mean field order parameter value, σ N and s are respectively the wire Drude conductance and cross section. At T long superconducting wires suffer a quantum phase transition [5] governed by the dimensionless parameter λ s to be defined later. In ultrathin wires with λ < superconductivity is completely destroyed by quantum fluctuations, and such wires can even go insulating at T =. In relatively thicker wires with λ > quantum fluctuations are not so pronounced, the wire resistance R decreases with T and one finds [5] R = d ˆV γqp S T λ 3, T Φ I, d I γqp S I λ 3, T Φ I. Here and below ˆV denotes the expectation value of the voltage operator across the wire. According to this result the wire non-linear resistance does not vanish down to lowest temperatures, as it was indeed observed in many experiments [7 1]. The physics behind this result is transparent: An external current I flowing along the wire breaks the symmetry between positive and negative voltage pulses making the former more likely than the latter. As a result, the net voltage drop ˆV occurs across the wire also implying non-zero resistance. Corresponding author andrei.zaikin@kit.edu 1 I.E.Tamm Department of Theoretical Physics, P.N.Lebedev Physics Institute, Moscow, Russia Institute of Nanotechnology, Karlsruhe Institute of Technology KIT, 761 Karlsruhe, Germany 3 National Research University Higher School of Economics, 11 Moscow, Russia Copyright line will be provided by the publisher 1

2 A. G. Semenov and A. D. Zaikin: Quantum phase slips and voltage fluctuations in superconducting nanowires FQMT-15 Can one also expect to observe non-vanishing voltage fluctuations in superconducting nanowires? The presence of QPS-induced equilibrium voltage fluctuations in such nanowires can be predicted already on the basis of the result combined with the fluctuation-dissipation theorem. The issue of non-equilibrium voltage fluctuations e.g., shot noise is more tricky. This issue requires a detailed theoretical analysis which is the main goal of the present paper. Below we will demonstrate that non-equilibrium voltage fluctuations in ultrathin superconducting wires are caused by the process of quantum tunneling of magnetic flux quanta Φ which can be described by Poisson statistics. In particular, we will investigate QPS-induced shot noise of the voltage in such wires and predict a highly non-trivial dependence of the noise power spectrum on temperature, frequency and external current. The system setup and the Hamiltonian In what follows we will consider an experimental setup depicted in Figure 1. It consists of a thin superconducting wire of length L and cross section s attached to a voltage source V x by means of a resistor R x. A capacitor C is switched in parallel to the superconducting wire. The right end of the wire is grounded as shown in the figure. Voltage fluctuations at its left end can be measured by a detector. At low enough energies the superconducting wire can be described by an effective Lagrangian [5, 6, 11] 1 Cw e d x ϕx, t 1 ϕx, t, 3 L kin x is the coordinate along the wire ranging from to L, C w denotes geometric wire capacitance per unit length and L kin = 1/πσ N s is the kinetic wire inductance times length. Below we will employ the well known property enabling one to describe any electric circuit either in terms of the phases node variables or in terms of the transferred charges loop variables, see, e.g., [1] for further details. It is remarkable that under certain conditions these two approaches can turn effectively dual to each other. The duality of this kind was established and discussed in details, e.g., for ultrasmall Josephson junctions [13 16], as well in the case of short [17] and long [18] superconducting wires. According to the results [18] the dual representation for the Hamiltonian of a superconducting nanowire is defined by an effective sine- Gordon model Ĥ wire = Ĥ T L + Ĥ QP S. I Vt Rx C Superconducting wire Figure 1 online color at: The system under consideration. Here Ĥ T L = L d x ˆΦ x L kin + ˆχx C w Φ defines the wire Hamiltonian in the absence of quantum phase slips. It describes an effective transmission line in terms of canonically conjugate flux or phase and charge operators obeying the commutation relation [ ˆΦx, ˆχx ]= i Φ δx x. 6 The term L Ĥ QP S = γ QP S d x cos ˆχx 7 accounts for QPS effects in our nanowire. Note that the quantum field χx, t is proportional to the total charge qx, t that has passed through the point x up to the time moment t, i.e. qx, t = χx, t /Φ. Hence, the local current I x, t and the local charge density ρx, t are defined by the equations I x, t = χx, t /Φ, ρx, t = χx, t /Φ. 8 In order to construct the total Hamiltonian Ĥ for our system it is also necessary to include the charging energy of a capacitor C as well as the Hamiltonian of the external circuit. Below we will assume that the external resistor R x is very large, i.e. the wire is biased by a constant current I = V x /R x which does not fluctuate. Then for the total Hamiltonian we have Ĥ = Ĥ wire + ˆQ ˆχ/Φ C Vx 5 I ˆϕ e. 9 Here the second and the third terms in the right-hand side represent respectively the charging energy and the Copyright line will be provided by the publisher

3 Fortschritte der Physik, January 1, 19 potential energy tilt produced by the current I. The operator ˆQ accounts for the charge across the capacitor and the phase operator ˆϕ corresponds to the variable ϕt ϕ, t which represents the phase of the superconducting order parameter x, t at x =. Setting ϕl, t, we conclude that the latter operator is related to the fluctuating voltage across the wire V t as ˆV t = ˆϕ/e. 1 V t 1 V t V t 3 = 1 8{ ˆV t 1 T ˆV t ˆV t 3 + T ˆV t ˆV t 3 ˆV t 1 + ˆV t T ˆV t 1 ˆV t 3 + T ˆV t 1 ˆV t 3 ˆV t + ˆV t 3 T ˆV t 1 ˆV t + T ˆV t 1 ˆV t ˆV t 3 + T ˆV t 1 ˆV t ˆV t 3 + T ˆV t 1 ˆV t ˆV t 3 }, 15 FQMT-15 3 Keldysh perturbation theory The task at hand is to investigate fluctuations of the voltage V t in the presence of quantum phase slips. We will proceed with the aid of the Keldysh path integral technique. As usually, we define our variables of interest on the forward and backward time parts of the Keldysh contour, ϕ F,B t and χ F,B x, t, and introduce the classical and quantum variables, respectively ϕ + t = ϕ F t + ϕ B t / and ϕ t =ϕ F t ϕ B t and similarly for the χ-fields. Making use of Eq. 1, let us express the general correlator of voltages in the form V t 1 V t...v t n = 1 e n ϕ + t 1 ϕ + t... ϕ + t n e i S QPS, 11 S QP S = γ QP S d t and... = D ϕt D χx, t...e i S [ϕ,χ] L d x sinχ + sinχ / 1 13 indicates averaging with the effective action S [ϕ,χ] corresponding to the Hamiltonian Ĥ = Ĥ Ĥ QP S. It is important to emphasize that Eq. 11 defines the symmetrized voltage correlators. E.g., for n= one has V t 1 V t = 1 ˆV t 1 ˆV t + ˆV t ˆV t 1, 1 while for n= 3 one can verify that [] T andt are, respectively, the forward and backward time ordering operators. Equation 11 is a formally exact expression which we will now evaluate perturbatively in the tunneling amplitude γ QP S 1. In the zero order in γ QP S the problem is described by the quadratic in both ϕ and χ action S. In that case it is necessary to employ the averages ϕ + t = ϕ t = χ x, t =, χ + x, t = Φ I t, 16 as well as pair averages the Green functions: G K ab X, X = i a + X b + X + i a + X b + X, G R ab X, X = i a + X b X, 17 ax and bx stand for one of the fields ϕt and χx, t. As both these fields are real, the advanced and retarded Green functions satisfy the condition G A ab ω= ω. Due to linearity the Keldysh function GK can G R ba then be expressed in the form G K ab ω= 1 coth ω G R ab ω GR ba ω 18 even in presence of the external bias V x. For the same reason the retarded Green function can be evaluated either from the full quantum mechanical treatment or just from simple electrotechnical arguments. For the function G R ϕϕ we obtain Gϕϕ R ω= 1 ω E C + i ω ωλ e R x π cot, 19 ωl v E C = e /C, v = 1/ L kin C w is the plasmon velocity [19] and the parameter λ already introduced above is defined as λ = R Q /Z w with R Q = π/e being the "superconducting" quantum resistance unit and Z w = Lkin /C w being the wire impedance. The corresponding expressions for G R ϕχ and GR χχ turn out somewhat lengthy. In order to simplify them it is useful to bear in mind that due to momentum conservation plasmons in our system can only be created in pairs with Copyright line will be provided by the publisher 3

4 A. G. Semenov and A. D. Zaikin: Quantum phase slips and voltage fluctuations in superconducting nanowires FQMT-15 Average voltage To begin with, let us evaluate the expectation value of the voltage operator. Our perturbation yields i γqp S ˆV = 8e d x d x lim ω ωgr ϕχ x;ω P x,x Φ I P x,x Φ I, P x,x ω=p x,x ω+ P x,x ω and P x,x ω= d t e i ωt e igx,x ;t, i G x,x;, i G x,x ;,, 3 G x, x ; t,= G K χχ x, x ; t,+ 1 GR χχ x, x ; t,. Figure Candy-like diagrams which determine both average voltage V upper diagram and voltage-voltage corrrelator V V six remaining diagrams in the second order in γ QP S. The fields ϕ +, χ + and χ in the propagators 17 are denoted respectively by wavy, solid and dashed lines. the total zero momentum. All plasmons moving towards the grounded end of the wire eventually disappear and never pop up again while excitations moving in the opposite direction produce voltage fluctuations measured by a detector. Then in the interesting for us long wire limit the general expressions for Gϕχ R and GR χχ reduce to Gϕχ R x;ω ω+i λe i ωx v ω E C + i λ π, Gχχ R x, x ;ω πi λ ω x x ω+i ei v. 1 In Eqs. and 1 we also set R x as requested in the current bias limit. Expanding Eq. 11 up to the second order in γ QP S and performing all necessary averages, we evaluate the results in terms of the Green functions 17. Pictorially these results can also be represented in the form of the socalled candy diagrams [1]. These diagrams for the first and the second moments of the voltage operator are displayed in Figure. Making use of the identity lim ω ωgϕχ R x;ω = πi, we can reduce Eq. to a simple form ˆV =Φ ΓQP S I Γ QP S I, 5 Γ QP S I = d x d x P x,x Φ I. 6 Employing the above results and making use of the generalized detailed balance condition P x,x ω=e ω T P x,x ω, 7 we obtain vlγqp S Φ V = ς Φ I sinh Φ I, 8 Γ λ ςω= τ λ i ω πt Γ λ + i ω πt πt λ 1, 9 Γλ τ 1/ is the QPS core size in time and Γx is the Euler Gamma-function. The result 8, 9 yields Eq. in the corresponding limits and matches with the analogous expression derived in [5] by means of a different technique. 5 Relation to ImF -method Comparing Eq. 5 with the corresponding result for the average voltage [5] we immediately conclude that the Copyright line will be provided by the publisher

5 Fortschritte der Physik, January 1, 19 Imt s Ret QPS pair [5] which yields γqp S δf 1/T d x d x dτe S pai r, 3 FQMT-15 Contour deformation -i S pai r = Φ Iτ+V x,τ; x,, 35 Figure 3 online color at: Integration contour. quantity Γ QP S I 6 can be interpreted as a quantum decay rate of the current state due to QPS. In [5] this rate was evaluated from the imaginary part of the free energy by means of the so-called ImF -method, see, e.g., []. It is instructive to establish a detailed relation between the latter approach and the Keldysh technique employed here. To this end let us define the generalized Green function G χ x, x ;σ which depends on the complex time σ and obeys the condition G χ x, x ; t i =G x, x ; t, at t >. It reads G χ x, x ;σ= i T d t cothπt t σ dω π e i ωt Gχχ R x, x ;ω Gχχ R x, x ; ω 3 This function is analytic, has branch cuts at Imσ = N/T for all integer N and is periodic in the imaginary time, i.e. G χ x, x ;σ=g χ x, x ;σ i/t. 31 Note that Eq. 31 just follows from the Kubo-Martin- Schwinger condition. On the imaginary axis the function G χ coincides with the Matsubara Green function G χ x, x ; i τ= ig M χχ x, x ;τ. 3 Evaluating the quantum decay rate Γ by means of the pioneered by Langer ImF -method one employs a general formula Γ= ImF, 33 F is the system free energy. In order to establish the QPS contribution to Γ it is necessary to evaluate the corresponding correction to the free energy δf. In the leading order in γ QP S it suffices to consider just one QPS-anti- τ is the imaginary time interval between QPS and anti-qps events and V x; x ;τ, describes the interaction between these two events occuring respectively at the points x and x. This interaction term can be expressed via the Matsubara Green function as V x; x ;τ,= G M χχ x, x ;τ 1 G M χχ x, x; 1 G M χχ x, x ;. 36 Note that the integral over τ in Eq. 3 is formally divergent at low temperatures. As a result, the free energy acquires an imaginary part ImF derived by means of a proper analytic continuation of δf. Evaluating the integral 3 by the steepest descent method we first determine a stationary point τ s from the stationary condition for the action Φ I = τ G M χχ x, x ;τ s. 37 A closer inspection demonstrates that this stationary point delivers a maximum to the action rather than a minimum, thereby indicating an instability with respect to QPS-mediated decay to lower energy states. In this case the correct recipe is to deform the integration contour along the steepest descent path. This procedure is illustrated in Figure 3. The initial integration path is a vertical line going from to i β note that the Matsubara technique operates with imaginary times. The deformed contour is directed along the real-time axis after passing through the point τ s. With this in mind we obtain γqp S δf τ s d x d x dτe Φ Iτ V x;x ;τ, d x d x i dτe Φ I τ s +i τ V x;x ;τ s +i τ,. 38 Copyright line will be provided by the publisher 5

6 A. G. Semenov and A. D. Zaikin: Quantum phase slips and voltage fluctuations in superconducting nanowires FQMT-15 Taking the imaginary part of this expression, we get ImF = γqp S + d x d x dτe Φ I τ s +i τ V x;x ;τ s +i τ, d x d x dτe Φ I τ s i τ V x;x ;τ s i τ,. 39 Further expressing Eq. 39 as a single integral along the contour passing through the point τ s in the direction perpendicular to real τ axis, we find ImF = d x d x dτe Φ I τ s +i τ V x;x ;τ s +i τ,. Combined with Eqs. 36 and 3, this expression can also be rewritten in the form ImF = d x d x d t e Φ I τ s +i t e ig χx,x ;t i τ s, i G χx,x;, i G χx,x ;,. 1 Employing the relation P x,x ω= d t e i ωt e ig χx,x ;t i i G χx,x;, i G χx,x ;,, we arrive at the final result ImF = d x d x P x,x Φ I. This result together with Eq. 33 demonstrates that Eq. 5 indeed defines the decay rate of the current states in a superconducting nanowire due to QPS, thereby proving the equivalence of the ImF -approach employed in [5] and the real time Keldysh technique combined with duality arguments elaborated here. The latter technique appears much more convenient for the analysis of voltage fluctuations to be developed below. 6 Voltage fluctuations power spectrum, i.e. S Ω = d t e i Ωt V t V =S Ω + Sr Ω + Sa Ω. 3 The first of these contributions S Ω = i Ω coth Ω 16e Gϕϕ R Ω GR ϕϕ Ω defines equilibrium voltage noise for a transmission line and has nothing to do with QPS. The remaining two contributions are due to QPS effects. The term S r Ω contains the products of two retarded advanced Green functions: γ SΩ r = QP S Ω coth Ω 16e d x [ d x Re Gϕχ R x;ω F x,x ΩG R ϕχ x ;Ω F x,x G R ϕχ x;ω ], 5 F x,x Ω= P x,x Ω+Φ I P x,x Ω Φ I + P x,x Ω+Φ I + P x,x Ω Φ I. 6 The term S a Ω, in contrast, contains the product of retarded and advanced Green functions and reads S a Ω = Ω 3e d x d x Gϕχ R x;ωgr ϕχ x ; Ω 7 [ Q ± P x,x Ω±Φ I P x,x Ω Φ I ]. ± Here we denoted Ω±Φ I Q ± = coth Ω coth. 8 Eqs. 3-8 together with the expressions for the Green functions 19-1 fully determine the voltage noise power spectrum of a superconducting nanowire in the perturbative in QPS regime. At non-zero bias values the QPS noise turns nonequilibrium. In the zero frequency limit Ω the terms S Ω and Sr Ω tend to zero, and the voltage noise S Ω S is determined solely by S a Ω. Then from Eq. 7 we obtain S = Φ ΓQP S I +Γ QP S I Φ I = Φ coth V, 9 Let us investigate the second moment of the voltage operator, i.e. the voltage noise. Our perturbative analysis allows to recover three different contributions to the noise V is defined in Eqs. 5, 8. In the low temperature limit T Φ I Eq. 9 accounts for QPS-induced shot noise S = Φ V obeying Poisson statistics with an 6 Copyright line will be provided by the publisher

7 Fortschritte der Physik, January 1, 19 effective charge equal to the flux quantum Φ. In other words, shot noise in superconducting nanowires is produced by quantum tunneling of the magnetic flux across the wire. In the dual picture tunneling flux quanta Φ can be viewed as charged quantum particles passing through and being scattered at an effective "tunnel barrier" which role is played by the nanowire. Let us also note that shot noise can also be generated in the high temperature regime dominated by thermally activated phase slips TAPS [3]. In the latter case this shot noise is again described by Eq. 9 with the TAPS decay rate substituted instead of the QPS one, i.e. Γ QP S I Γ T AP S I. Our analysis also allows to recover higher correlators of the voltage operator 11. Let us define the voltage cumulants C n = i n lim 1 t n z log t t e i z d t 1 V t 1 z= 5 In fact, within the accuracy of our analysis the terms C k γqp S with k < n generated in the right-hand side of Eq. 5 can be safely dropped and C n just coincides with the Fourier transformed correlators 11, i.e. C = S etc. As before, proceeding perturbatively in γ QP S and making use of Eqs. 8, 9, at T we find π vlγ C n = Φ n 1 QP S τλ Φn V = λ Γ Φ I λ. 51 λ Let us now consider voltage fluctuations at non-zero frequencies. Below we will restrict our analysis to voltage noise and stick to the limit of sufficiently high frequencies and/or long wires v/l Ω. In this case the term S Ω turns out to be independent of the wire length L. Evaluating the QPS terms SΩ r and Sa Ω, we observe that the latter scales linearly with the wire length L as the former does not. Hence, the terms S Ω and Sr Ω can be safely neglected in the long wire limit. For the remaining QPS contribution S a Ω we get S a Ω = vlλ 8e [ Φ I ς Ω Φ I ς + Ω Φ I ς Φ I ] sinh Ω/EC + λ/π sinh Ω At T from Eq. 5 we find S a Ω I λ 1 I Ω/Φ λ 1, Ω<Φ I /,, Ω>Φ I / Figure online color at: The temperature dependence of the QPS noise spectrum S Ω 5 at λ = 3., large E C and different Ω in the long wire limit. The inset shows S Ω as a function of Ω at different temperatures. This result is explained as follows. At T = each QPS event can in general excite N plasmons N = 1,... with total energy E = Φ I and total zero momentum. N plasmons carrying total energy E/ propagate towards the grounded end of the wire and eventually get dissipated there, while the remaining N plasmons also with total energy E/ propagate in the opposite direction reaching the opposite wire end and causing voltage fluctuations emit a photon with frequency Ω measured by a detector. At T = this process is only possible at Ω < E/ in the agreement with Eq. 53. The result 5 is also illustrated in Figure. At sufficiently small Ω we still keep Ω v/l one observes a non-monotonous dependence of S Ω on T which is a direct consequence of quantum coherent nature of QPS noise. Note, that at non-zero T the expression 5 does not coincide with the zero frequency result 9 even in the limit Ω. The point here is that before taking the zero frequency limit in Eq. 5 one should formally set L. Then one gets vlt γ S a Ω I = QP S Φ ς Φ I ς Φ I sinh Φ I. 5 Comparing this expression with Eq. 9 combined with Eqs. 8, 9 obtained in the true zero frequency limit meaning that the limit Ω was taken prior to sending the wire length L to infinity one can establish the iden- FQMT-15 Copyright line will be provided by the publisher 7

8 A. G. Semenov and A. D. Zaikin: Quantum phase slips and voltage fluctuations in superconducting nanowires FQMT-15 tity S I,T S a Ω I,T = RI,T. 55 It follows immediately that both expressions 5 and 9 coincide only at T =, while at any non-zero T the noise power S I,T 9 exceeds one in Eq. 5 and in contrast to the latter increases monotonously with temperature. Finally, we point out that the perturbative in γ QP S approach employed here is applicable for not too thin wires with λ >, i.e. in the "superconducting" phase. In thinner wires with λ < characterized by unbound QPS-anti- QPS pairs a non-superconducting phase the QPS amplitude γ QP S gets effectively renormalized to higher values [5] and, hence, the perturbation theory becomes obsolete in the low energy limit. However, even in this case our results may still remain applicable at sufficiently high temperature, frequency and/or current values. This work was supported by the Russian Science Foundation under grant No [1] D. V. Averin and A. A. Odintsov, Phys. Lett. A 1, [15] A. D. Zaikin, J. Low Temp. Phys. 8, [16] G. Schön and A. D. Zaikin, Phys. Rep. 198, [17] J. E. Mooij and Yu. V. Nazarov, Nat. Phys., [18] A. G. Semenov and A. D. Zaikin, Phys. Rev. B 88, [19] J. E. Mooij and G. Schön, Phys. Rev. Lett. 55, [] A. V. Galaktionov, D. S. Golubev, and A. D. Zaikin, Phys. Rev. B 68, [1] A. G. Semenov and A. D. Zaikin, to be published 16. [] U. Weiss, Quantum Dissipative Systems World Scientific, Singapore, 8, p. 7. [3] D.S. Golubev and A.D. Zaikin, Phys. Rev. B 78, Key words. Quantum phase slips, shot noise. References [1] K. Yu. Arutyunov, D. S. Golubev, and A. D. Zaikin, Phys. Rep. 6, 1 8. [] A. Bezryadin, J. Phys.: Condens. Matter, 3 8. [3] A. D. Zaikin, in: Handbook of Nanophysics: Nanotubes and Nanowires CRC Press, Boca Raton, FL, 1 p. -1. [] A. Bezryadin, Superconductivity in Nanowires Wiley- VCH, Weinheim, 13. [5] A. D. Zaikin, D. S. Golubev, A. van Otterlo, and G. T. Zimanyi, Phys. Rev. Lett. 78, [6] D.S. Golubev and A. D. Zaikin, Phys. Rev. B 6, [7] A. Bezryadin, C. N. Lau, and M. Tinkham, Nature, 971. [8] C. N. Lau, N. Markovic, M. Bockrath, A. Bezryadin, and M. Tinkham, Phys. Rev. Lett. 87, [9] M. Zgirski, K. P. Riikonen, V. Touboltsev, and K. Y. Arutyunov, Phys. Rev. B 77, [1] X. D. A. Baumans, D. Cerbu, O-A. Adami, V. S. Zharinov, N. Verellen, G. Papari, J. E. Scheerder, G. Zhang, V. V. Moshchalkov, A. V. Silhanek, J. Van de Vondel, Nat. Commun. 7, [11] A. van Otterlo, D. S. Golubev, A. D. Zaikin, and G. Blatter, Eur. Phys. J. B 1, [1] A. M. Zagoskin, Quantum Engineering Cambridge University Press, New York, 11. [13] S. V. Panyukov and A. D. Zaikin, J. Low Temp. Phys. 73, Copyright line will be provided by the publisher