arxiv: v1 [cond-mat.supr-con] 27 Dec 2018

Transcription

1 Ccution of currents round tringur indenttion by the hodogrph method J. I. Avi 1, B. Venderheyden 2, A. V. Sihnek 3, S. Meinte 1 1 ICTEAM, Université cthoique de Louvin, B-1348, Louvin-L-Neuve, Begium 2 SUPRATECS nd Deprtment of Eectric Engineering nd Computer Science, Université de Liège, B-4 Srt Timn, Begium nd 3 Experiment Physics of Nnostructured Mteris, Q-MAT, CESAM, Université de Liège, B-4 Srt Timn, Begium (Dted: December 31, 218) rxiv: v1 [cond-mt.supr-con] 27 Dec 218 Border indenttions in non-iner conductors, such s superconducting thin fims in the creep regime, ter the distribution of currents nd mgnetic fieds ner nd fr from the indenttion. One of such disturbnces re the discontinuity ines, or d-ines, prboic-ike ine originting from the indenttion where the current density direction chnges brupty. Hodogrph series resuts re obtined for the currents round tringur indenttion nd its corresponding d-ines in conducting stripe of finite width nd in n infinite hf pne, considering two cses: uniform creep exponent nd mixed infinite nd ohmic exponents. The mixed creep exponent cse presents currents distributions resembing the purey ohmic cse, with significnt current disturbnces ony ner the indenttion. For uniform creep exponent, resuts simir to pnr indenttion re obtined, with fr rnged currents fetures nd prboic-ike d-ines with shpes depending on the creep exponent. In prticur, the sme d-ine symptotic behviour is obtined for the tringe indenttion s tht of the pnr defect in the critic stte, resut obtined here just on continuity considertions of the hodogrph expnsions, nd quite simir but not identic to the popur Ben mode. This equivence is due to identic contributions to the Fourier series of the current strem-function in the hodogrph spce obtined from n imges method expnsion. INTRODUCTION Defects hve significnt effects on the eectric properties of type II superconductor. Whie nnoscopic defects serve s pinning centres tht rrest the motion of fux vortices nd void dissiption [1], rger defects crete zones of incresed resistivity tht cn serve s nucetion spots for fux vnches [2, 3]. Border defects such s crcks or voids re ccompnied by ong rnge perturbtions in the current density nd mgnetic fied simir to the Ben mode [4] discontinuity ines, where the current density J turns shrpy. The exct shpe of the ines wi be directy reted to the critic current nd constitutive retionship under extern ppied fieds, such s n extern eectric fied, or n extern mgnetic fied rmp creting induction currents. The resistivity in the fux creep regime [5] cn be modeed s: ρ(j, B, T ) = ρ (J/J c (B, T )) n1, (1) where ρ is the resistivity t J = J c, B the mgnetic fied, T the temperture nd n is the creep exponent. The steep vritions of the resistivity in this regime re so espected to generte thermo-mgnetic instbiities [6, 7]. By tking no mgnetic fied or temperture dependence of J c, expression (1) hs been used by mens of the hodogrph method [8 1] to obtin nytic expressions for the eectric currents round pnr defect, for finite nd infinite creep exponent n. A more bsic mode consists of considering ony the cse n nd to disregrd non-oc effects nd mgnetic fied or temperture depence of the critic current. This eds to the Ben mode [4, 11], in which the mgnetiztion current density of smpe exposed to rmping mgnetic fied is modeed s piece-wise uniform regions, on whose boundries the current cn chnge of direction discontinuousy, trcing the d-ines. From this, prboic d-ines of the form x = y 2 /2h h/2 re obtined for n indenttion of height h on border octed t x =. Here, the ccutions of currents round pnr defect by the hodogrph method [9] re extended to tringur defect. In tht formism the ccuted currents ctuy do not chnge direction in discontinuous wy, however, the d-ine concept is extended to current domin w tht trces ine simir to tht of the Ben s mode, so the d-ine term is so used here. The current strem function is obtined for different creep exponents nd geometries nd finy compred to the Ben mode. STRIPE WITH TRIANGULAR DEFECT A gener sttionry soution method ws deveoped in refs [8 1] for currents in fims with non-iner retionship between the eectric fied E nd the current J of the type E = E (J/ ) n, where E is chrcteristic eectric fied mgnitude nd the corresponding current density nd n 1. In this method the Mxwe equtions for the current strem function ψ re soved in the spce of the current density mgnitude J = ψ nd the nge θ subtended between J nd the x xis, i.e. under the trnsformtion ψ(x, y) ψ(j, θ). It cn be equy ppied for the eectric potenti with φ(x, y) φ(e, θ), where E = φ. In the spce (J, θ), nmey the hodogrph spce, the eectromgnetic equtions become iner nd seprbe when using the power w resistivity, nd cn be soved by eigenfunction expnsions, with gener

2 2 soution: ψ(j, θ) I with ( J = A (B D θ) ) (1n) C θ J m C m sin(mθ φm) J m D m sin(mθ ωm), (2) τ ± m = 1 2 (1 n ± (n 1) 2 4nm 2 ). (3) The boundry conditions for ψ in the hodogrph spce re simpe ony in the cse of smpes with stright borders, in which θ is constnt. Once the soution is found in the hodogrph pne, it cn be trnsformed bck to the x, y spce, i.e. (J, θ) (x, y). By using the compex vribe z = x iy, this trnsformtion is given by: z = F Ceiθ J n eiθ J {ib (n 1) D (1 i(n 1)θ} n e iθ D m 1 τm [m cos(mθ ω m) 1 e iθ iτm sin(mθ ω ] m) C m 1 τ m [m cos(mθ φ m) 1 iτ m sin(mθ φ m) ], (4) The technique hs been ppied to geometries ike corners, bridges, current eds nd pnr defects. In Ref. [9] pnr defect of width 2 perpendicur to the current in n infinite superconductor ws considered. In tht work, the d-ine ws identified s the core of current domin w strting ner the defect. In the imit n it ws found there tht, symptoticy, the d-ine trces curve with x y 2 /1.94, quite simir to the d-ine in the Ben resut tht goes s x = y 2 /2. STRIPE OF FINITE WIDTH Equtions (2) nd (4) re soved for strip crrying tot current I, of width nd with tringur constriction of bse 2 nd height, s depicted in Fig.1(). The soution for the strem function ψ(j, θ) is seprted in 3 hodogrph regions, s shown in Fig.1(b), nmey ψ 1 in region (1) for J, with the mgnitude of J fr from the defect, ψ 2 in region (2) for J, where is the mgnitude of the current density in the border directy opposite to the defect, nd ψ 3 in region (3) for J. The soutions fufiing the boundry conditions described in Fig.1(b) re: ψ 1 I = 1 J sin(θ), (5) () A 2 B C J E D y x (b) J C θ ψ= D ψ=i ψ= E A B ψ=i (3) (2) (1) ψ=i /2 3/4 FIG. 1. () Geometry of superconducting strip with tringur constriction, indicting with cpits the points with direct equivent in hodogrph spce. In prticur the current density is zero in the corner B nd reches t (, ) nd fr from the constriction. The nge subtended by segment CB is 3/4. The strem function ψ is chosen to be on the right nd I on the eft boundries of the strip. Ony the upper prt y of the smpe is soved, the remining is obtined by symmetry. (b) Hodogrph pne representtion of soution for ψ ong with the boundry conditions foowing exmpe of pnr constriction in ref. [1]. ψ 3 I ψ 2 I 4θ = 2 = 1 The, D(2), C(2) [ ( J θ ) τ C ψ=i ψ=i ψ=i B ] C (2) J sin(θ), (6) D (3) J sin(( 2)θ ), (7) 2 nd D(3) 2 re obtined by enforcing continuity of ψ nd ψ/ J t J = nd J = ong with Fourier nysis for θ [/2, 3/4]. Appying continuity t : 8 3/4 /2 It gets: ψ 1() sin(θ)dθ = 8 3/4 /2 ψ 2() sin(θ)dθ (8) = 2 m D(2) rτ C (2) (9) Where r = /. From the derivtive continuity t, it is obtined: = τ τ rτ C (2) (1)

3 3 From continuity t, by using: 8 3/4 sin(( 2)θ 3 /2 2 ) sin(θ)dθ = 8 It gets: C (2) = 2 8 m D(2) rτ () 2 ( 2) 2 (11) D (3) 2 () 2 ( 2) 2 (12) And from continuity of the derivtive t : τ D(2) rτ τ C (2) = 8 Combining (9) nd (1): = 2 m τ τ τ D (3) 2 τ 2 () 2 ( 2) 2 (13) (14) And combining (12), (13) nd (14) the iner system of equtions for D (3) 2 is obtined: τ () 2 (rτ 1) = D (3) 2 (τ τ 2 ) () 2 ( 2) 2 (15) The spce coordintes (x, y) ssocited to ψ i (J, θ) cn be obtined from the soution (4) for the compex vribe z = x iy in the different regions: z 1 = 1 i eiθ J 1 τ 1 [ cos(θ) iτ sin(θ) ], (16) Where the integrtion constnts hve been chosen to gurntee spce continuity of ψ. The prmeters re dependent on D (3) 2, nd those depend t the sme time on r, tht is found numericy by enforcing spce continuity of the soutions. In the imit n the coefficients τ m m 2 nd τm. This grety simpify the strem function ccution. In prticur the hodogrph region (3) nd r = ( )/ due to current conservtion. Then C (2) = 2/m from (12) nd = 2(r()2 1)/m from (9). The current stremines, or the sme the contour or equipotenti ines of the strem function, cn be seen for different n vues in FIG. 2. Those noticby chnge direction round the tringe constriction much fr from the constriction thn in the ohmic cse. The non-oc mgnetic fied H d given by Biot-Svrt s w hs we behved oc minimum round the sme re of current direction chnge, nd s in [2], this minimum ws used here s the d-ines, tht re dispyed in FIG. 3. This pproch is prctic s in mgneto-optic experiments it is the mgnetic fied tht is mesured rther thn the current distribution. The fied H d ws ccuted by Fourier trnsform method [12], nmey, H d = F 1 [kf(ψ)/2], where F is the 2D Fourier trnsform, nd k the wve-vectors ssocited to the domin discretiztion. As it cn be seen in FIG. 3, the d-ines re wys nrrower thn the one corresponding to the Ben cse. This is expected s in stripe it is impossibe to keep constnt current density with constrictions present. Then, if J, > J c, fux fow regime must be considered ner the indenttion, where for finite n the current density J diverges. z 2 z 3 ( ) = i e iθ r e iθ = ( ) r e iθ 4 J eiθ ( J τ 1 iτ sin(θ) ] C (2) τ 1 ( J ) τ 1 [ cos(θ) ) τ 1 [ cos(θ) iτ sin(θ) ], (17) D (3) 2 J 2 1 τ 2 1 [ ( 2) cos(( 2)θ 3 2 ) iτ 2 sin(( 2)θ 3 2 ) ], (18) Fux Fow The previous hodogrph ccution ssumes tht the power w resistivity is vid even for currents bove the current, s in refs. [8 1]. If we consider the fux fow resistnce, then better representtion of the resistivity is given by: ρ = { ρ (J/J c ) n1, J J c ρ, J > J c, (19) Assuming tht = J c for this cse, here it is combined the creep regime for regions with J <, nd ohmic regions for J > using n = 1. Lets obtin first the retionship for the series expnsions coefficients for region (1) with rbitrry creep exponent ñ with corresponding τ m = 1 2 (1ñ (ñ 1) 2 4ñm 2 ), ong with region (2) nd region (3) both with creep exponent n nd τ ± m = 1 2 (1 n ± (n 1) 2 4nm 2 ). The continuity retionships (9), (12) nd (13) sty unmodified. Then by enforcing z 1 = z 2 it gets:

4 4 () (b) (c) (d) FIG. 2. Current stremines in bck nd mgnetic fied contours in coor in smpe with tringur defect of height.5 nd bse.1 ) Ohmic cse n = 1, r = b) n = 51, r=.9593 c) Critic stte or n, r =.95. d) Hybrid cse with n for J < nd n = 1 for J, r = τ τ τ 1 = D(2) τ 1 C (2) rτ τ τ 1 8 () 2 1 (2) Tht reduces to (1) when ñ = n. And by combining eqns. (9), (12), (13) nd (2) n eqution cn be obtined to ccute D (3) 2. Here is tken the prticur cse with ñ nd n = 1. With tht choice, the retionships for the coefficients re finy: = 2 m D(2) r C (2) (21) = r C (2) (22) C (2) = 2 m D(2) r 8 C (2) = D(2) r 8 D (3) 2 () 2 ( 2) 2 (23) D (3) 2 ( 2) () 2 ( 2) 2 (24) Combing the bove equtions (22), (23) nd (24) it is obtined: r 8m 1 = D (3) ( 2 (1 r 8m ) ( 2)(1 r 8m ) ), () 2 ( 2) 2 (25) tht ows ccuting D (3) 2 nd then D(2), C(2) nd. The r prmeter is gin found numericy by enforcing continuity of the different soutions, nd the

5 5 y n=19 S n=31 S n=11 S n= S Ben H1 S H2 HP n=19 HP n=11 HP n= HP x FIG. 3. Minimum in nonoc mgnetic fied for different creep exponents in stripe (S) nd infinite hf pne (HP). Hybrid cse (H1 nd H2) correspond to n for J < nd n = 1 for J. Vues of r prmeter used re.9716,.9648,.95475,.95 nd for n=19 S, n=31 S, n=11 S, n = S nd H1 S respectivey. Ben corresponds to Ben mode prbo. current stremines nd mgnetic fied cn be seen in FIG.2(d). The d-ine for this cse cn so be found in FIG.3. Appying continuity conditions for ψ nd the terms invoving the sine functions of z it gets: τ 2 2 = 8 τ2 1D(2) 2 = 8 8 Combing eqns. (3) nd (31): () 2 ( 2) 2 (3) ( 1)( 3) () 2 ( 2) 2 τ τ 1 (31) τ 2 () 2 ( 2) 2 τ2 1 τ τ 1 ( 1)( 3) When n this expression simpifies to: 1 () 2 ( 2) 2 () 2 = 1 ( 1)( 3) Tht cn be better inverted s: = (32) (33) DEFECT IN INFINITE HALF PLANE A strem function for the infinite hf pne x with tringur defect tht stisfies the boundry conditions is given by: ψ 1 = For J <. And for J : ψ 2 = J sin(θ), (26) J τ sin(( 2)θ ), (27) 2 Notice the different units of distnce of the C nd D coefficients now. The z functions re given by: z 1 = 1 i e iθ J τ 1 τ 1 [ cos(θ) ] iτ sin(θ), (28) = (()2 1) () 2 ( 2) 2 ( 1)( 3) = 32 k (34) Buy using the orthogonity properties of the expnsion of expression of the form sin(θ) m = sin((2)θ3/2)b k, nd the sum on the right is ccuted expicity. The currents stremines nd mgnetic fied cn be seen in FIG.4(). Fux fow Agin the cse with creep exponent ñ for J < nd n for J is investigted. Continuity conditions now produce: 2 = 8 () 2 ( 2) 2 (35) e iθ z 2 = (1 ) 2 τ 2 1 ( J ) τ 2 1 [ ( 2) cos(( 2)θ 3 2 ) iτ 2 sin(( 2)θ 3 2 ) ], (29) τ 2 τ2 1D(2) 2 = 8 8 ( 1)( 3) () 2 ( 2) 2 τ τ 1 (36)

6 6 And the eqution for : τ 2 () 2 ( 2) 2 τ 2 1 τ τ 1 = (37) ( 1)( 3) And gin for the hybrid cse, very high ñ is considered for which the imit ñ is tken in the bove expression, nd n = 1, which gives: 2 () 2 ( 2) = 1 ( 1)( 3) (38) The current stremines nd mgnetic fieds nd d-ine cn be found in FIG.4(b) nd FIG.3 respectivey. DISCUSSION The tringur defect here studied shres sever properties with the pnr defect nyzed in [9]. In prticur, the current density diverges on the tip of the tringe for finite n due to the existence of the hodogrph region (3), opposed to the critic stte n in which the current density remins finite in tht spot. For the infinite hf pne it cn be verified tht J = t the tringe tip, by evuting z 1 (θ = 3/4 δ, J = ) from eqution (28) with the id of expression (34) for : z 1 = 1 i e i( 3 4 δ) 8 ie i( 3 4 δ) (1)m cos(δ) (1)m sin(δ), (39) And by tking the imit δ it is obtined z 1 = 1, probing tht J = there. In the cse of the critic stte in the stripe, it cn be equy ccuted tht J = t the tip of the tringe, the sme s on the opposite side of the stripe. In fct J = on the ine connecting both points, since J = on both sides of the constriction, nd the current is pssing through zone where the mximum J is, nd therefore r = ( )/. In the hybrid cse for the stripe, most of the current significnty bove psses ner the tringe tip, nd the current density J is very cose to in most of the constriction. Additiony, is so very cose to, s it cn be seen from the r = vue. In the cse of the infinite hf pne with hybrid creep exponents, it is seen tht perturbtions on the currents re ony importnt ner the indenttion, simir to the purey ohmic cse, virtuy disppering t distnces rger thn 1, s it cn be seen in the current stremines in FIG.2(d) nd their d-ines in FIG.4(b). The d- ines for the stripe nd infinite pne so most identic in tht cse, just chnging in their horizont extension. The d-ine in the critic stte for the infinite hf pne trces curve quite cose to the Ben prbo, but identic ctuy to the resut obtined in [9] of x y 2 /1.94 for the pnr defect. In the ppendix 1 it is shown tht both the pnr defect nd tringur indenttion hve ctuy the sme symptotic d-ine. This is done by repcing the expression (34) for, obtined here just from continuity conditions, opposed to the ccution for the pnr defect [9] where is obtined from expnsion of current vortices nd nti-vortices. The sme procedure is so done in the ppendix 1 to obtin pproximte expressions for the nd D(2) 2, tht hve the sme function form nd imit for high n s the pnr defect. A key point of the derivtion is tht the behviour of the currents is miny determined by vortex soution to London eqution in trnsformed hodogrph pne of coordintes η = (1/ n) og(e/e ), θ. Tht vortex is centred in η =, θ = /2, nd is the min contributor to the Fourier expnsions of ψ in both the pnr defect nd in the tringe indenttion. CONCLUSION The hodogrph series method is used to sove the probem of currents circuting in stright stripe of finite width nd n infinite hf pne with tringur border indenttion, showing the non-oc mgnetic fieds produced by those currents. Both cses present simir properties s the pnr defect redy discussed in [9, 1], with diverging currents in the defect tip for finite n, nd bounded current for the critic stte. An hybrid cse of fux-fow nd creep ws so considered, showing the presence of d-ines, but of short rnge. For the tringur defect in the infinite hf pne, the sme symptotic d-ine is obtined s tht of the pnr defect in the critic stte, which is expined here by the simir contribution of vortex currents in modified hodogrph spce η, θ. ACKNOWLEDGEMENTS Grnt ARC 13/18-8 for Concerted Reserch Actions, finnced by the Woni-Brusses Federtion. J. I. A. so cknowedges FNRS support reference 216/V 3/5/13 (W) -IB/LS-6967.

7 7 () (b) FIG. 4. Current stremines in bck nd mgnetic fied contours in coor in hf infinite pne with tringur defect of height.5 nd bse.1 ) Critic stte or n. b) Hybrid cse with n for J < nd n = 1 for J. Mgnetic fieds re ccuted with currents in re twice the size shown. Appendix 1: Asymptotics of currents for hf infinite pne As per [8 1], eq. (2) is the soution of the eqution ppering from the hodogrph trnsformtion: J 2 n 2 ψ J 2 J ψ J 2 ψ θ 2 = (4) Whie the trnsform expressions for spce coordintes (4) re obtined from the equtions: z Eσ E = ψ ψ eiθ J 2 iρ (41) θ E ) z ( θ = eiθ E2 ψ J E iρ ψ (42) θ Where σ = J/ E. The foowing chnge of vribe cn be done for (4): η = (1/ n) og(e/e ) = n og( J ) (43) ψ = e βη h (44) Where β = n1 2, nd eq.(4) trnsforms into London n eqution: 2 h η 2 2 h θ 2 β2 h = (45) Considering η nd θ s crtesin coordintes, chnge of vribe to por coordintes r, χ is performed, with η = r cos χ, θ = θ sin χ. A vortex soution is obtined in these por coordintes: h = 2h e βη r cos χ 2 (46) The soution for ψ cn be constructed by dding different vortices of the sme strength on the ine η = nd ternting signs such tht ψ = t θ = 4, 3 4,... This is done by dding terms of the form (46) with sign (1) m1 centred t points θ m = m 2, nd converting bck to the η, θ coordintes: ψ = h (1) m1 rm η e β(rmη) (47) r m m Where r m = η 2 (θ m/2). The coefficients cn be obtined then by ccuting the Fourier terms: = 8 3/4 /2 ψ(, θ) sin(θ)dθ (48) The min contributor to tht integr is the term with m = 1 of (47), so the other m terms re discrded. For n 1, ψ(, θ) is very peked round θ = /2 so the upper imit cn be extended to infinite, nd the foowing expressions re obtined: = h = h 8 2 β2 () 2 β β2 () 2 (49) β2 ( 2) 2 β β2 ( 2) 2 (5) Where the sme principe ws used to obtin (5). The prmeter h is obtined by compring the behviour of (49) for high n: h 32m 2, (51) n 3/4 to the vue redy ccuted here in (34) from continuity conditions for n nd equ to 32m/,

8 8 then: h = n3/4 2 (52) The d-ine cn be obtined from the symptotic behviour of the current, tht cn be obtined by dding up the series (28) when repcing (34), for n : z 1 = 1 i 8 () 2 1 J eiθ () 2 1 [ cos(θ) i() 2 sin(θ) ], (53) This expression cn be written s function of the series I invoving cosine functions: z 1 = 1 i 8 ( eiθ I i I ) (54) θ I = = () 2 () 2 1 J cos(θ) () 2 1 (55) ( ) () J 1 cos(θ) () 2 1 (56) Then, for τ = og( /J) >. For z, J nd τ 1, nd (J/ ) ()2 1 = e τ(()2 1) e τ()2. I ( ) 1 1 () 2 e τ()2 cos(θ) (57) 1 With b = 8/. By repcing I nd I, nming e 1 = exp(γ 2 /4τ) /8 τ: x = be 1 (sin γ γ cos γ) (64) 2τ y = be 1 (cos γ γ sin γ) (65) 2τ Negecting the terms not proportion to e 1. Then, eiminting e 1, nd soving for τ: τ = γ 2 (x sin γ y cos γ) (x cos γ y sin γ) (66) Combining (64) nd (65) nd repcing bck (66), nd expnding for γ (1 ξγ) (1ξγ) b 2 x2 = (ξ γ) 3 γ eγ (ξγ) (67) In which ξ = y/x. As discussed in [9], the d-ine corresponds to the ine where the spce derivtives of γ(x, y) diverge. This cn be obtined from the condition x/ γ =. Then (67) is pproximted for ξ 1, derived with respect to γ whie keeping ξ constnt nd equed to zero, from which it is obtined γ = ξ/2. Repcing this bck in (67) nd retining the owest order in ξ, the expression for the d-ine is obtined x = y 2 /(8/ 2e). REFERENCES I e τ()2 cos(θ) e τ()2 () 2 1 cos(θ) (58) The exponenti term in the second sum cn be omitted for τ 1, nd repcing the nge γ = θ /2: I e τ()2 cos(γ) cos(γ) () 2 1 (59) The first sum cn be estimted with the Euer-Mcurin formu, giving exp(γ 2 /4τ) /8 τ 1/2 nd the second dds to 1/2 sin(2γ)/4. Then: I 8 2 τ eγ /4τ sin 2γ (6) 4 I = I γ θ eγ 16τ 3/2 Then, for y, x 1, from (54): 2 /4τ cos 2γ (61) 2 x = b (I sin γ I cos γ) (62) y = b (I cos γ I sin γ) (63) [1] E. H. Brndt, Rep. Prog. Phys. 58, 1465 (1995). [2] J. Brisbois, O.-A. Admi, J. I. Avi, M. Mott, W. A. Ortiz, N. D. Nguyen, P. Vnderbemden, B. Vnderheyden, R. B. G. Krmer, A. V. Sihnek, Phys. Rev. B 93, (216). [3] R. G. Mints, E. H. Brndt, Phys. Rev. B 54, (1996). [4] C. P. Ben, Rev. Mod. Phys. 36, 31 (1964 ). [5] V. M. Vinokur, M. V. Feige mn, V. B. Geshkenbein, Phys. Rev. Lett. 67, 915 (1991). [6] A. L. Rkhmnov, D. V. Shntsev, Y. M. Gperin, T. H. Johnsen, Phys. Rev. B 7, (24). [7] D. V. Denisov, A. L. Rkhmnov, D. V. Shntsev, Y. M. Gperin, T. H. Johnsen, Phys. Rev. B 73, (26). [8] A. Gurevich, J. McDond, Phys. Rev. Lett. 81, 2546 (1998). [9] A. Gurevich, M. Friesen, Phys. Rev. B 62, 44 (2). [1] M. Friesen, A. Gurevich, Phys. Rev. B 63, (21). [11] A. M. Cmpbe, J. E. Evetts, Critic Currents in Superconductors, Monogrphs on Physics. Tyor & Frncis Ltd., London (1972). [12] J. I. Vestgården, P. Mikheenko, Y. M. Gperin, T. H. Johnsen, New J. Phys. 15, 931 (213 ).