## Brazilian Journal of Physics

*versión impresa* ISSN 0103-9733

### Braz. J. Phys. vol.40 no.4 São Paulo dic. 2010

#### http://dx.doi.org/10.1590/S0103-97332010000400015

**Vortex dynamics equation in type-II superconductors in a temperature gradient **

**R. Vega Monroy ^{I,*}; J. Sarmiento Castillo^{I}; D. Puerta Torres^{II}**

^{I}Facultad de Ciencias Básicas. Universidad del Atlántico Km. 7, Via a Pto. Colombia, Barranquilla, Colombia

^{II}Facultad de Ciencias Exactas. Universidad de Cartagena Plaza de la Artillería N. 30-84, Cartagena, Colombia

**ABSTRACT**

In this work we determined a vortex dynamics equation in a temperature gradient in the frame of the time dependent Ginzburg-Landau equation. In this sense, we derived a local solvability condition, which governs the vortex dynamics. Also, we calculated the explicit form for the force coefficients, which are the keys for the understanding of the balance equation due to vortex interactions with the environment.

**Keywords:** Type-II Superconductors; Vortex Equation; Vortex Balance Equation.

**1. INTRODUCTION**

In type-II superconductors, the vortices are the ones in charge for the magnetic properties of these systems since every vortex carries a magnetic flux quantum. In the last years the interest in the vortex motion is associated to many non peculiar properties in HTSC not found in conventional type-II superconductors. In particular, one of the most important effects encountered in HTSC is the Hall anomaly [1,2]. In this sense, it is known that, under the action of the Lorentzs force vortices acquire a movement and therefore losses appear in the superconducting state. The equation of motion, which governs the vortex dynamics in type-II superconductors, has been subject to a great amount of works, which have helped us to understand this phenomenon in these systems. In addition, the interest in the vortex motion is emphasized by the responsibility of this dynamics in a great variety of transport phenomena in type II superconductors. Generally, the vortex dynamics has been considered on the hydrodynamical two fluid model [3,4], where the relative motion between the superfluid and the vortex generates the Magnus force. Next, the normal component reacts to this motion, producing the longitudinal viscous drag force and the transversal Iordanski's force, which are the two components of the medium force. An attempt to describe the vortex dynamics was done by Dorsey [5] in the frame of the time dependent Ginzburg-Landau equation following previous developments of Gorkov and Kopnin [6]. In this work, Dorsey formulates a solubility condition through with a vortex equation can be obtained. Today, there are several attempts to construct a unified theory about the vortex motion. Some approximations use the sophisticated many body formalism [7,8,9]. Other authors apply a simpler theory based on the kinetic Boltzmann equation to study the dynamic behavior of the vortex structures [10], but so far the vortex dynamics is an open question for the solid state physics community. The purpose of the present work is to contribute to a better understanding of the fascinating phenomenon of the vortex motion. In this connection, the goal of this paper is to determine a vortex dynamics equation in a temperature gradient in the frame of the time dependent Ginzburg-Landau equation. Such kind of equation has been introduced in a heuristic way in many works [11] to satisfy experimental data [12,13,14].

The paper is organized as follows: In Section 2, we have obtained the basics equations of the work. In Section 3, some dynamic coefficients are calculated and finally in Section 4 we summarized the main results.

**2. BASIC EQUATIONS**

The present analysis follows the works done by Dorsey [5] in order to obtain the equation, which describes the vortex dynamics. Let us write the dimensionless time dependent Ginzburg-Landau equation for the complex order parameter in the form:

In the above equation γ is the dimensionless relaxation time, k is the Ginzburg-Landau parameter, Φ and are the electric potential and the magnetic vector respectively. The order parameter Ψ in terms of the amplitude *f*(,t) and the fase χ(,t) can be represented as follows Ψ(*,t*) = *f*(,*t*)exp[*i*χ(,*t*)]. The relaxation time has a complex character and can be written

The appearance of the imaginary part in the previous expression is a necessary condition for the gauge invariance conservation of equation (1) [10]. Relaxation processes that entail to dispersion in the vortex dynamics are of two types: the first is associated to the Bardeen- Stephens mechanism of dissipation and the second process is an intrinsic relaxation mechanism that governs the approach of the order parameter to its equilibrium state due to variations in the chemical potential. This mechanism is associated to the change of the order parameter in time because the electrons are forced to pair and de-pair with respect to different potentials due to the vortex motion. In consequence, this processe is the responsible that the relaxation time acquires a complex character as was shown by Kopnin [10] and it is framed in the parameter γ_{2}

where ν and µ are the density of states and the chemical potential respectively. Introducing Ψ and γ in the expression (1) and in addition, introducing the invariant forms for the magnetic vector and the scalar potential

one can obtain two equations, the first for the real part and the other one for the imaginary part as follows:

To form a closed system of equations, we need to derive an equation for the magnetic vector. In this connection, the dimensionless equations for the superconducting and normal current are:

where σ^{(n)} and *b*^{(n)} are the electric and thermoelectric conductivities in the normal state respectively. In this sense, the dimensionless equation for the magnetic vector is:

Now, to find an equation for the potential P, one can use the vector relation ·(*f*^{2}) = 2*f**f*·+*f*^{2}·, so that counting the fact that ·(_{s}+* _{n}*) = 0, i.e.

we obtain from equations (5) and (6) the following relation:

The above equation together with relations (3) and (7) will allow to arrive at the solvability condition for the equation that determines the vortex dynamics in type II superconductors.

On the other hand, locally in the vortex, the order parameter differs from the value in the bulk of superconductor so that a deviation appears in the order parameter and the potentials. In this connection, the quantities *f* and are expanded

where *f*_{0} and *f*_{1} are the amplitudes of the order parameter associated to the equilibrium and non equilibrium states respectively. In addition, in the vortex motion it is possible to consider that vortices move independently in the first approximation of the limit B Hc2, so that one can find an equation of motion for every individual vortex in presence of a temperature gradient. In this limit, assuming that the vortices move uniformly, the quantities* f*, and P are functions solely of -_{L}*t*, where is the electron position vector and _{L} is the vortex velocity. Thus, the temporary derivative can be written in terms of spatial derivative by means of = -_{L}·, so if we introduce this relation in (3), (7) and (9), we obtain the following set of equations for the equilibrium state:

being

the equilibrium current density. For the non equilibrium state we have the equations

where we have considered the fact that *f*_{0} *f*_{1} , _{0}_{1} and the linear approximation of the deviation. In addition we introduced the quantities

The Ginzburg-Landau equations are translational invariant, and therefore *f*_{0}(+) and _{0}(+) are the solutions to the equation (10), being an arbitrary translation vector. If we expand in Taylors series these two quantities, the linear equations (13) and (14), without the nonhomogenous terms, can be written in the form

where we introduced the parameters * f_{d}* ≡ ·

*f*

_{0}and

_{d}≡ ( ·)

_{0}. In order to derive the solvability condition one can integrate the equations (13) and (16) over a cylindrical area. Using the in the plane Gauss theorem and neglecting superficial contributions, the following equations are obtained:

Now we multiply the equation (18) by *f*_{d} and the equation (19) by *f*_{1}. Next both equations are combined to obtain the following expression:

Multiplying the quantity _{1s} by _{d}, and the quantity _{d} by _{1}, then these two expressions are subtracted each other and replacing in the equation (20) to obtain:

In order to simplify, in the k 1 limit, from equations (2) and the expression for _{d}, we have _{d} ≈ - χ_{d} / k and ≈ - χ_{1}/ k . Inserting the above expressions in the left side of the equation (21) we obtain:

Now considering the vector relation ·( α ) = ·α + α· and in addition, ·_{d} = 0, we can apply the in the plane divergence theorem as follows:

On the other hand, inside the vortex, the total current density satisfies the continuity equation:

where we have considered that the total charge density is ρ = ρ_{n}+ ρ* _{s}*, being ρ

_{n}and ρ

*the normal and superconducting charge densities respectively. In this way, from equation (11) we have:*

_{s}So that the expression (24) take the form:

Replacing the expressions (21) and (25) in the equation (23) and assuming that _{n} homogeneous, we have is

The equation (26) is the local solvability condition for the vortex motion. Also it is of the linear order in the vortex velocity. If this condition is not fulfilled, then the non homogenous equations do not have solution. The expression (26) contains an original contribution of this paper and will allow determining the equation that governs the vortex dynamics in type II superconductors in presence of a temperature gradient. The solvability condition obtained in the present work differs from Dorsey's expression [5] by the thermal flow force due to the temperature gradient and also by normal contributions (which give a right relation of the involved dynamic parameters) in the vortex dynamics.

In order to derive the vortex dynamics equation, it is necessary to introduce a coordinate system. In this sense, as much the applied temperature gradient as the transport current are oriented in x direction. The magnetic field is directed along the z-axis. The vortex velocity vector and the displacement vector form angles θ* _{H}* and ø with the x axis respectively. The geometry of the problem is pictured in Figure 1.

From Figure 1 we have:

The differential equation for the scalar potential *P* acquires the following form:

In the above equation _{0} = x_{0}. It is possible to determine the form of the solution to the equation (31) considering the limit r → 0, then χ → θ . Therefore, from the equation (2) we have P ≈ - ≈ _{L}.θ =*ê*_{θ} .

In general way, the solution to the non homogenous differential equation (31) takes the form:

where *p*_{1} is the solution to the homogenous equation (31), *p*_{2} and *p*_{3} are particular solutions to the nonhomogenous equation. From the previous expression, we have that, *p*_{1}, *p*_{2} and *p*_{3} satisfy the differential equations:

As *T* acts on normal electrons for small values of *r*(*r* → 0) then *f*_{0}^{2} → 0 and for high values of r superelectrons do not feel the presence of the temperature gradient, reason why in the equation (35) the second term of the left side can be neglected, having the solution (32) the form:

where the parameter ε = was introduced.

**3. VORTEX DYNAMICS EQUATION**

The solvability condition can be evaluated solving the equation (26). The surface integral can be expressed in terms of the applied transport current at the borders, where

so that the first integral to evaluate is the only one in the left side of the equation (26), where we use the formulae for the transport current * _{t}* and χ

*in cylindrical coordinates shown in the equations (27) and the trigonometric equality sin(θ - ø ) = sinθcosø - cosθsinø,*

_{d}The following term to solve is the first integral of the right side of the equation (26), which with the aid of the relation (36) is written as follows

Developing the product *f*_{0}*f*_{d} and using the expression for , given in equation (30), we have

Solving the second integral in the right side of the equation (26), one must account the relations (12) and (15), obtaining

where we have introduced the projections for _{L} and given by the equations (29) and (30) respectively. Now to calculate the third term of the equation (26), we expand ρ_{n}() in a Taylor's series around the vortex centre

Accounting that, the position vector in cylindrical coordinates has the form = r[sin(θ - θ* _{H}*)

*ê*+cos(θ - θ

_{r}*)*

_{H}*ê*

_{r}_{θ}] and taking the scalar product, the function ρ

*() can be written in the following way*

_{n} Taking the gradient of the function (42) in the limit *r* → 0, we obtain:

where we have considered that ρ* _{n}*() also grows in the -

*ê*direction. Replacing the equation (43) in the solvability condition, the third term of this expression acquires the form:

_{r}Equally, the term associated to _{L} in the solvability condition is obtained as follows:

Finally, the fourth integral in equation (26) is solved using the relations (37) and (29):

The equation of motion is obtained from the solvability condition when we replace the results obtained in equations (38), (40), (41), (44), (45) and (46). In a compact form the equation of motion can be written:

In the above equation, is the quantum of circulation and we have taken account that _{t} = * _{s}*/2 and

_{s}= ρ

_{s}_{s}. In addition we have defined the coefficients

The equation (47) determines the vortex dynamics in type II superconductors in presence of a temperature gradient. Equations of this type have been reported in some works in a phenomenological way [11]. In the present development this equation is obtained in a natural way from the time dependent Ginzburg - Landau equation. The term of the left side of equation (47) represents the Magnus force, which is connected to the momentum transfer between the vortex and the the superfluid. The forces proportional to *D* and *D'* are due to dispersion of quasi particles by the vortices, presenting a momentum transfer between the vortex moving with velocity _{n} and the normal fluid moving with velocity _{n}. The transverse component of this interaction is the Iordanskii force. The forces proportional to *d* and *d'* , according to Kopnin and Kravtsov [15], are due to vortex relaxation processes associated to the vortex interaction impurities which are at rest in the lattice. The transverse component of this interaction is the Kopnin-Kravtsov force. The force proportional to *T* is connected to the vortex thermal flow with a transport entropy *S _{v}*. If a moving vortex transports entropy, it experiences a force in a temperature gradient. The excistence of this force is well known from experimental data [12,13,14].

**4. CALCULATION OF THE COEFFICIENTS d,d' AND S_{v}**

In this section we will determine the explicit form for the coefficients, which are connected to relaxation processes in the vortex dispersion so as the entropy *S _{v}* due to the thermal motion. In general, when the Ginzburg-Landau equation is solved, one take some conjectures about the profile of the order parameter inside the vortex. An approach of the variation of dimensionless order parameter, which we will use, was proposed by Brandt [16]. In the present analysis we will find solutions in the limit

*r*→ 0 where the order parameter has the form

Where ξ_{0} is the coherence length. Replacing the equation (53) in the equation (33), the following equation in cylindrical coordinates is obtained:

where the parameter R = was introduced. The differential equation (54) has a solution in power's series given by the expession *p*_{1}(*r*) = , where *q* = -1. The asymptotic solution to the equation (54) must consider that:

From the previous analysis *p*_{1}(*r*) ≈ 1/*r* in the limit *r* → 0. In the same way for the magnetic vector we have *Q*_{0}(r) ≈ . But on the other hand, *Q*_{0} satisfies the relation *h*_{0} = x_{0}, where *Q*_{0}(*r*) =*h _{0}*(0)

*r*which determines the solution Q

_{0}(r) = , where h

_{0}(0) is the magnetic field in the vortex center. In this sense, accounting the relation (55) in the limit

*r*→ 0, then we have

*a*

_{1}= 0 and therefore

In order to determine *p*_{2}(*r*) in the same way that for *p*_{1}(*r*), we look to the solution of the equation (32), in power series and doing all the previous procedure one obtain

Now we are able to determine coefficients *d , d*' and *S _{v}*. Replacing the expressions (56) and (57) in the equations (49), (50) and (51) in the limit

*r*→ 0:

From the above equations we observe that the parameter *d*, which is associated to the interaction with impurities, has a double nature, the first one is connected to the Bardeen-Stephen mechanism of dissipation, whereas the second contribution is due to variations in the local charge density in the vortex interaction with impurities. The same effect occurs in the obtaining of the coefficient *d' *. Equally it is important to notice that the transport entropy posses a direct relation to relaxation processes due to the condensation and Coopers pair breaking by different pairing potentials via vortex motion as it is shown in equation (60).

**5. SUMMARY**

In the present analysis, we have obtained a vortex dynamics equation in type-II superconductors. Following previous works by Dorsey [5], Gorkov and Kopnin [6] we developed a more general equation for the vortex motion in the frame of the time dependent Ginzburg-Landau Equation. After some assumptions, we derived a local solvability condition, which governs the vortex dynamics. In this sense, we calculated the explicit form for the force coefficients, which are the keys for the understanding of the balance equation due to vortex interactions with the environment. In this connection, we took into account the spatial distribution of the normal electrons in the vortex core. This fact allowed to obtain exact expressions for the coefficients *D* and *D'* , which characterize the interaction between vortices and normal electrons by their dispersion. On the other hand, the approximated expressions for the coefficients *d* and *d'* , which characterize the vortex interaction with impurities were obtained. The nature of these coefficients is connected to the energy dissipation due to local variations of the charge density in the vortex core by the electron pairing. The same source of dissipation is associated to the transport entropy in the temperature gradient.

**Acknowledgements**

The Authors would like to thank the financial support granted by the Universidad del Atlántico.

[1] J.M. Harris, N.P. Ong and Y.F. Yan, Phys. Rev. Lett. **71**, 1455 (1993). [ Links ]

[2] S.J.Hagen et al., Phys. Rev. B **47**, 1064 (1993). [ Links ]

[3] J. Bardeen and M.J. Stephen, Phys. Rev. **140**, A1197 (1965). [ Links ]

[4] P. Nizierese and W.F. Vinen, Philos. Mag. **14**, 667 (1966). [ Links ]

[5] A.T. Dorsey, Phys. Rev. B **46**, 8376 (1992). [ Links ]

[6] L.P. Gorkov and N.B. Kopnin, Sov. Phys. Usp. **18**, 496 (1975). [ Links ]

[7] P. Ao and Thouless, Phys. Rev. Lett. **70**, 2158 (1993) [ Links ]

[8] P. Ao and Q Niu, Phys. Rev. Lett. **76**, 3758 (1996) [ Links ]

[9] G. Blatter, V.G. Geshkenbein and N.B. Kopnin, Phys. Rev B **59**, 14663 (1999) [ Links ]

[10] N.B. Kopnin, Rep. Prog. Phys. **65**, 1633 (2002) [ Links ]

[11] A. Freimuth and M. Zittartz, Phys. Rev. Lett. **84**, 4978 (2000). [ Links ]

[12] R.P. Huebner, Supercond. Sci. Technol. **8**, 189 (1995). [ Links ]

[13] T.T.M. Palstra et al., Phys. Rev Lett. **64**, 3090 (1990). [ Links ]

[14] M. Zeh et al., Phys. Rev. Lett. **64**, 3195 (1990). [ Links ]

[15] N.B. Kopnin and V.E. Kravtsov, JETF Lett. **23**, 578 (1976). [ Links ]

[16] E.H. Brandt, J. Low Temp. Phys. **24**, 709 (1977). [ Links ]

(Received on 27 August, 2010)

* Electronic address: rvega@uniatlantico.edu.co