arxiv: v2 [cond-mat.supr-con] 5 Apr 2017

Transcription

1 Vortex Dynamics in Type II Superconductors Dachuan Lu (Kuang Yaming Honors School, Nanjing University) (Dated: April 6, 2017) Time dependent Ginzburg-Landau equation is solved for type II superconductors numerically, and the dynamics of entering vortices, geometric defects and pinning effects have been investigated. A superconducting wire with ratchet defects is designed to pump the vortices move in a specific direction and enhance the supercurrent when applying the periodical magnetic field. Some properties of this wire have been investigated numerically and analytically. arxiv: v2 [cond-mat.supr-con] 5 Apr 2017 I. INTRODUCTION High temperature superconductors(htsc) are widely used in industrial production and medical equipments. Almost all the HTSCs are type II superconductors, the vortex pinning effects and vortex dynamic need concerning[1]. Ginzburg and Landau proposed a phenomenological theory for superconducting phase[2] based on Landau second-order phase transitions, and the so-called Ginzburg-Landau equation is widely used to study the vortex in type II superconductors after Abrikosov firstly predicted the vortex lattice in type II superconductors [1][3]. Geometry and pinning has large effects on the configuration of vortex lattice. Time dependent Ginzburg-Landau (TDGL) equation has more fascinating properties than the original stationary one[4], and can be used to investigate the dynamic properties of the vortices. Mu[5] and Du[6] have solved the TDGL equation by finite element method, Pedersen use COMSOL(finite element solver) to solve the TDGL in complicated geometry[7]. Some other methods are widely used to investigate the vortex dynamics. Reichhardt has used molecular dynamics(md) to simulate the microscopic behavior of flux and calculate some macroscopic quantities, such as magnetization hysteresis and critical current density[8] [9] [10]. They modeled the vortex-vortex and vortex-pin force and calculated the total force exerted to each vortex. However, Ginzburg- Landau theory is more suitable for the situation with varying magnetic field, because more vortices will penetrate into the sample while increasing the magnetic field, and molecular dynamics method could hardly calculate the creation and annihilation of the vortices. In Section. II, we normalize the TDGL to the simplest form and In Section. III, we present some numerical results based on the normalized TDGL, the dynamics of entering vortices, geometric defects and pinning effects have been investigated. In Section. IV, we investigate the flux pump which is an application based on our former investigation and can pump the vortices to move in specific direction. In Section. V, we give a summary of this paper and clarify some limitation of our work. II. TIME DEPENDENT GINZBURG-LANDAU EQUATION In 1950, Ginzburg and Landau proposed the famous phenomenological theory[2] to describe the superconductors, in which, complex order parameter is introduced to describe the Cooper pair and normal electrons, a φ 4 like potential is considered which accounts for the phase transition. Ginzburg-Landau theory can precisely describe the magnetic field penetration through the quantized magnetic flux vortices formation in the type II superconductors. To find the dynamics of the vortex in superconductors, we turn to the time dependent Ginzburg Landau (TDGL) equation which is generalized by Schmid[4]. The TDGL equation is read as, ( 2mD t + i q ) Φ Ψ = 1 ( ) 2 2m i qa Ψ + αψ β Ψ 2 Ψ (1) ( ) A σ t + Φ = q 2mi (Ψ Ψ Ψ Ψ ) q2 m Ψ 2 A 1 ( A B a ) (2) µ 0 where Ψ is the order parameter, A and Φ are the vector potential and scalar potential of the electromagnetical field, B a is the applied magnetic field, m and q are the mass and charge of the Cooper pair respectively, and σ is the conductance of the normal electrons. α, β, D are the parameter dependent on the properties of superconductors, in which β is assumed to be a constant while α(t ) = α(0)(1 T T c ) is dependent on the temperature T.

2 2 The LHS(left hand side) of Eq. 1 and Eq. 2 are time dependent term, and the RHS are original Ginzburg- Landau theory. To describe the superconductor in vacuum, we choose the following appropriate boundary conditions, ( i Ψ qaψ) n = 0 A = B a ( (3) A t + Φ) n = 0 The first and the last equation in Eq. 3 implies that the supercurrent and the normal current is zero at the boundary, and the second equation illustrates that the magnetic induction is equal to the applied magnetic field at the boundary. by introducing London penetration depth λ, Ginzburg- Landau coherence length ξ and Ginzburg-Landau parameter κ = λ/ξ. The corresponding quantities are shown as follows, (x, y, z) (λx, λy, λz ), t ξ2 D t A = qξ A, Φ = αdκ 2 2µ 0 b Φ, (4) Ψ = α β Ψ, σ = 1 µ 0Dκ σ 2 A. Normalization For simplicity in simulation, the parameters in Eq. 1 and Eq. 2 can be arranged into dimensionless quantities By applying the normalization procedure, the TDGL equations can be read as, ( ) ( ) 2 i t + iκφ Ψ = κ + A Ψ + Ψ Ψ 2 Ψ (5) ( ) A σ t + Φ = 1 2iκ (Ψ Ψ Ψ Ψ ) Ψ 2 A ( A B a ) (6) and the boundary conditions are, ( i Ψ + AΨ) n = 0, A = B a, ( A t + Φ) n = 0, B. Gauge Invariance (7) Because the free energy is invariant under the U(1) gauge transformation, Ψ = Ψe iκχ, A = A + χ, Φ = Φ χ t (8) where χ is an arbitrary function, we can fix the gauge by varying the χ and make the scalar potential Φ to be zero, namely, χ t = Φ. The terms containing Φ in the Ginzburg-Landau equation(eq. 5, Eq. 6) and the boundary conditions (Eq. 7) can be dropped. The U(1) gauge invariance implies that the particle number is conserved in this system, for the Hamiltonian is commuted with the particle number operator ˆQ = Ψ n Ψ n. This property indicates that there is n no dissipation in this system. III. NUMERICAL SIMULATION In this section, we investigate the numerical solution of dimensionless TDGL(Eq. 5, Eq. 6) by using the finite element method which is implemented by COMSOL [11]. For simplicity but with universality, we consider 2D geometry because the superconducting layers of many HTSCs have strong 2D properties. We set the superconducting layer in xy plane and designate (x, y, t) as the spatial and temporal coordinate. The magnetic field is set along the z direction. In numerical implementation, the order parameter is divided into real part and imaginary part, Ψ = RΨ+iIΨ, the vector potential is decomposed into scalar functions in different directions, namely A = A x i + A y j. In the following section, we will first investigate some universal properties of the vortex dynamics, and then turn to the flux pump, which is the application of the vortex dynamics and intends to enhance the supercurrent.

3 3 A. Disk Sample with SO(2) symmetry The disk remains unchanged under the rotation of any angle, and when we apply the magnetic field along the z- axis as described in the Section. III, the vortices will form from the boundary and get into the sample, consistent with the some previous works[12][13]. Because of the system tends to the minimum free energy, the vortices will be rearranged to reduce the interaction energy, but the original symmetry is broken. The evolution process of the vortices is shown in Fig. 1. We also calculate the energy density of this process by using the following equations, H sup = 1 κ 2 Ψ 2 Ψ Ψ 4 H mag = (B a A) 2 H int = i κ A (( Ψ) Ψ Ψ ( Ψ )) + A 2 Ψ 2 (9) where H sup is the superconducting energy, H mag is the magnetic energy and H int is the interaction energy. The total energy is, H tot = H sup + H mag + H int (10) FIG. 1. The evolution of vortices penetrating into the disk sample. we set κ = 4, B a = 1.5, σ = 1 and the diameter is d = 4, they are all dimensionless parameter. The density plot illustrates the magnitude of order parameter, the higher magnitude, the lighter color, (a) the magnetic field penetrates into the sample uniformly and destroy the Meissner state, (b) the magnetic field tends to exist in form of flux vortices, (c) the vortices get into the sample, (d) the vortices are rearranged to a low energy configuration. The vortices tends to repel each other, and the repulsive force is Φ ( 2πµ 0 λ K x 3 1 λ), which decays with the distance x dramatically, as Fig. 2 illustrates. The first ( order modified bessel function of the second kind, K r 1 λ) λ r ( as r 0 and K r ) 1 λ πλ 2r e r/λ as r. The corresponding short range interaction is large, while the long range interaction decays fast, so we could set a cutoff in numerical simulation. FIG. 3. Different parts of energy change along with time. The blue line is the superconducting energy, increases dramatically when the vortices are entering the sample. The red line is the interaction energy and has a maximum at the point when the vortices are entering the sample. The grey line is the magnetic energy and decreases along with the magnetic field penetrating into the sample and forming vortices. The yellow line is the total energy and always decreases to find the lowest energy state. 2 1 Φ 0 f/ 2 π μ 0 λ B. Disk Sample with Defect x/λ FIG. 2. Repulsive force of the adjacent vortices From the previous section, we can gain some general idea of how vortices form and evolve, in this section, we will investigate the sample with geometrical defects which will break the symmetry and make the motion of vortices anisotropy.

4 where t = T/T c is reduced temperature, and ɛ(r) = 1, 0 models the pinning center. As Fig. 5 illustrated, vortex will be attracted by the defects which are often approximately modeled by an attractive parabolic wells[8]. 4 FIG. 4. Comparison of vortex dynamics in the samples with two types of geometrical defects. As illustrated in the Fig. 4, we find that it is easier for vortices to form at the defects with the obtuse angle, while the defects with acute angle tend to remain Meissner state. This phenomenon can be explained by the repulsive force between two vortices. For the acute angle situation, the vortices forming from two edges will repel each other, and the total force has non-zero component along the bisector of the acute angle, so the vortices tend to avoid the defects with acute angle. Meanwhile, the vortices tend to be attracted by the surface for the Bean-Livingston image force[14], as also observed in the experiment[15]. Due to the boundary condition of the superconductor, there will be an image anti-vortex at the opposite position of the interior vortex. And the interior vortex will be attracted by the image anti-vortex, but the interaction is short range, the vortex will only be attracted to the boundary when it is very close to the boundary, and will not feel the boundary when the distance from the boundary is larger than the magnetic field penetration depth λ. C. Vortex Pinning When the superconductors have at least one of defects, disorders and impurities, the magnetic field will get into the sample in form of vortices and pinned. By varying the parameters κ, α, β and changing the geometry of the samples, we can simulate different pinning states and investigate their properties[16]. One model to describe the δt pinning is achieved by modifying α, namely the coefficient of the ψ term in dimensionless TDGL, to be dependent on temperature. ( t + iκφ) Ψ = ( i κ + A) 2 Ψ + ɛ (r) (1 t) Ψ Ψ 2 (11) Ψ FIG. 5. When increasing the applied magnetic field, the vortex interacts with the defect. The defect obviously breaks the symmetry of the disk, and due to its attraction, the magnetic field will form a vortex from the boundary at the position which is the closet to the defect. And the vortex will site at the defect, until the magnetic field increased higher than the pinning attraction potential. IV. FLUX PUMP High temperature superconductors are widely used to carry large current. However, the movable vortices in type II superconductors may cause large dissipation effects, in a result, reducing the critical current. The vortex dynamics deserve attention, because most of the HTSC are type II superconductors, and when the applied magnetic field is higher than H c1, the superconductors will change into mixed state, in which movable vortices appears. In this section, we report using periodical magnetic field to enhance the supercurrent in a type-ii superconducting wire with ratchet defects, which may correspond to vortex-rectification effect[17]. A. Geometry and Dynamics of the Ratchet Superconducting Wire We designed the superconducting wire with strong unidirectional properties. As Fig. 6 illustrates, the infinite long superconducting wire is built by the blocks with the size of 2λ 2λ, the top and bottom of each block are applied the periodical boundary condition. The geometrical defects and the two sides of the blocks are applied the boundary conditions of TDGL, as Eq. 7 suggests.

5 5 FIG. 8. Motion of vortex in the superconducting wire with extra periodical magnetic field, parameters are B 0 = 1.3, B c = 0.8, ω = 0.12.(a) t = 57, (b) t = 70, (c) t = 85, (d) t = 100. FIG. 6. Geometry of superconducting wire. We argue this geometry could easily make the vortices form and move in an efficient way and can be simply generalized to other geometry. As discussed in Section. III B, we find that the defects with obtuse angle could easily form vortices, so the vortex will easily form at the angle A(Fig. 7) when we turn on the magnetic field B 0 along z direction. For the size of each block is comparable with the size of vortex, there will only form one vortex in each block. B. Vortex Transportation and Supercurrent By integrating the current density along the top of each block, we can find the supercurrent in the superconducting wire, J s = 1 2iκ (Ψ Ψ Ψ Ψ ) Ψ 2 A, (12) I s = top J s dl. (13) By comparing with supercurrent in the superconducing wire without the ratchet-like defects, we find that the ratchet-like defects do enhance the unidirectional properties of the superconducting wire(fig. 9), the average supercurrent of the wire without ratchet-like defects is almost zero, but of the one with defects has non-zero value. The corresponding phenomenon of the superconducting wire without the ratchet-like defects is that the vortex didn t move from one block to another, the vortex only vibrates along with the periodical magnetic field B p. FIG. 7. Geometry of superconducting wire. To make the vortices move continuously, we add a periodical magnetic field along z direction with wave like form, B p = B c sin( 2π L y ωt)k, where B c, ω is the magnitude and frequency of the periodical magnetic field, L is the length along y direction of each block. The total applied magnetic field is B a = B 0 k + B p k. Fig. 8 shows how the vortex move in the superconducting wire. FIG. 9. Supercurrent of wire with ratchet and without ratchet, the +, are the ratchet direction, the ratchet in Fig. 6 is +. The blue, red and grey line are Ratchet(+), Ratchet(-) and No Ratchet respectively.

6 6 By tuning the frequency of B p to match the motion of vortex, we can make the vortex move in the highest speed, correspondingly, the current in the superconducting wire becomes largest. By using the parametric sweep, we can find the optimal frequency of B p (Fig. 10) before, ( the interaction decays exponentially in long range but K r 1 λ) λ r at small r. Eq. 15 can be approximated as, f B = 1 2 B c 2πµ 0 λ 2 Eq. 14 can be simplified as, x (t) = n nearest s kt + x (t) r r B,n r r B,n 2 (16) (17) where k = ωl/2π and s = 1 2 B c 2πµ 0η λ. By solving Eq. 17, 2 we can find the average velocity is, v e k(l 2x 0 ) s k 2 (l x 0 ) s kx 0 e k(l 2x 0 ) s (s + k ( l + x 0 )) (18) FIG. 10. ω 0.12 is the optimal frequency for this superconducting wire. where x 0 is the initial position of the moving vortex. The supercurrent density can be characterized by the average velocity according to, C. Analytical Analysis of the Optimal Frequency To find the optimal frequency of B p, we turn to investigate the dynamics of single vortex. The overdamped equation for the single vortex is, ηv = f B + f vv (14) where η BH c2 /c 2 ρ n (ρ n is the resistivity of normalstate) is the damping coefficient due to the dissipation processes inside and around the vortex cores[18], v is the velocity of the vortex, f B is the force exerted by the magnetic field, and f vv is the force between two vortices. But for the periodical boundary condition, the vortexvortex force can be neglected. The periodical magnetic field will make the vortices move along the hypotenuse of the ratchets when the magnitude B c is below the Bean-Livingston surface barrier H s. The force f B is actually the type of vortex-vortex force[8], f B = 1 2 B c 2πµ 0 λ 3 K 1 ( r r B,n ) r r B,n λ r r B,n n (15) where r B,n = (ωlt/2π + nl)j + ωl tan θt/2πi, n Z is the position vector of the n th vortex which is generated by the periodical magnetic field B p, r is the position vector of the moving vortex. As shown in Fig. 2 and discussed j = 1 2i (Ψ ( p m Ψ) + ( p m Ψ) Ψ) = Re(Ψ ( p m Ψ)) ρv (19) And when k satisfies k 2s L x 0, the velocity reaches the maximum. V. SUMMARY In this study, time-dependent Ginzburg-Landau equation is solved in different geometrical regions numerically by finite element method. The distribution of the order parameter is vividly depicted and the magnetic vortices can be found. Some derived quantities, such as magnetization and current density, are also calculated. By modifying the parameters in the TDGL equation, the pinning effect is investigated. In the last part, we investigate a special superconducting wire with ratchet-like defects, in which the unidirectional motion of the vortices can enhance the supercurrent. However, some limitations in this study need to be fixed in the future. The TDGL equations(eq. 1, Eq. 2) used in this study are only suitable for the s-wave superconductors such as niobium, but most of the type II superconductors, such as BSCCO and YBCO, are d-wave or even more complex pairing symmetry, thus, Eq. 1, Eq. 2 are no longer suitable. Some researchers have generalized the Ginzburg-Landau equation for d-wave and p-wave superconductors[19][20], these can be adapted for our future study. [1] G. Blatter, M. V. Feigel man, V. B. Geshkenbein, A. I. Larkin, and V. M. Vinokur, Reviews of Modern Physics 66, 1125 (1994). [2] L. D. Landau and V. Ginzburg, Zh. Eksp. Teor. Fiz. 20,

7 (1950). [3] A. A. Abrikosov, in Dokl. Akad. Nauk SSSR, Vol. 86 (1952) pp [4] A. Schmid, Physik der kondensierten Materie 5, 302 (1966). [5] M. Mu, SIAM Journal on Scientific Computing 18, 1028 (1997). [6] Q. Du, Journal of mathematical physics 46, (2005). [7] T. S. Alstrøm, M. P. Sørensen, N. F. Pedersen, and S. Madsen, Acta applicandae mathematicae 115, 63 (2011). [8] C. Reichhardt, C. Olson, J. Groth, S. Field, and F. Nori, Physical Review B 53, R8898 (1996). [9] R. Richardson, O. Pla, and F. Nori, Physical review letters 72, 1268 (1994). [10] C. Reichhardt, C. Olson, and F. Nori, Physical Review B 57, 7937 (1998). [11] W. B. Zimmerman, Multiphysics Modeling With Finite Element Methods (series on Stability, vibration and control of systems, serie) (World Scientific Publishing Co., Inc., 2006). [12] V. Schweigert and F. Peeters, Physical review letters 83, 2409 (1999). [13] B. Baelus, F. Peeters, and V. Schweigert, Physical Review B 63, (2001). [14] P. S. Deo, V. Schweigert, and F. Peeters, Physical Review B 59, 6039 (1999). [15] A. Geim, S. Dubonos, I. Grigorieva, K. Novoselov, F. Peeters, and V. Schweigert, Nature 407, 55 (2000). [16] A. E. Koshelev, I. A. Sadovskyy, C. L. Phillips, and A. Glatz, Physical Review B 93, (2016). [17] J. Van de Vondel, C. de Souza Silva, B. Zhu, M. Morelle, and V. Moshchalkov, Physical review letters 94, (2005). [18] J. Bardeen and M. Stephen, Physical Review 140, A1197 (1965). [19] J.-H. Xu, Y. Ren, and C.-S. Ting, Physical Review B 53, R2991 (1996). [20] J.-X. Zhu, C. Ting, J. Shen, and Z. Wang, Physical Review B 56, (1997). [21] H. K. Onnes, Commun. Phys. Lab. Univ. Leiden 12, 1 (1911). [22] A. K. Geim, S. V. Dubonos, I. V. Grigorieva, and K. S. Novoselov,, 55 (2000). [23] C. C. de Souza Silva, J. Van de Vondel, M. Morelle, and V. V. Moshchalkov, Nature 440, 651 (2006).