^{1}

^{2}

A plasma has been defined as a quasi-neutral gas of charged particles showing collective behaviour. Plasmas can support waves depending on the local condi-tions, the presence of external electric and magnetic fields. A characteristic property of plasmas is their ability to transfer momentum and energy via col-lective motion. An example in this case, is the Langmuir waves where plasma electrons oscillate against a st ationary ion background. In a superconductor, two plasma electrons arise one that is made up of normal electrons and the other that is made up of super-electrons. In this study, we consider a system of super-electrons forming a super-particle. The motion of the plasma su-per-particles around a magnetic vortex core has been studied in the YBCO123 and Bi2212 systems.

Ginzburg-Landau (GL) theory gives a good account of interactions between magnetic field and electric field in a type II superconductor [_{c} superconductors. While the lower and the upper critical fields of a superconductor have been determined using the GL theory, finding the radius of the magnetic vortex core theoretically has been a challenge. Therefore, there is need for a new approach in studying the motion of charged particles in an electromagnetic field, especially in high-superconductor systems. In this study, we consider assemblage of charged particles, each moving independently in prescribed electromagnetic fields which constitute a plasma [^{th} particle in motion is given by

F i = m d v i d t (1)

where v i = d x i d t is the velocity of the i^{th} particle. The equation of motion of a charged particle in a magnetic field ( B ) due to an electric field ( E ) is given by [

d v i d t = q m ( E + v × B ) (2)

where v i = d x i d t is the velocity of the i^{th} particle, q is the particle charge and F i = q ( E + v i × B ) is the Lorentz force on the i^{th} particle. The charge and current densities respectively are obtained by integrating the distribution function over velocity space and summing over species

ρ ( x , t ) = ∑ i q i ∫ d 3 v ⋅ f i ( x , v , t ) (3)

And

J ( x , t ) = ∑ i q i ∫ d 3 v ⋅ v f i ( x , v , t ) (4)

where f = f ( x , v , t ) is the microscopic distribution function whose time derivative is given by

d f d t = ( ∂ f ∂ t + ∂ x i ∂ t ∂ f ∂ x + ∂ v i ∂ t ∂ f ∂ v ) ∑ i = 1 N δ ( x − x i ( t ) ) δ ( v − v i ( t ) ) (5)

here, v = | v | Putting Equation (2) into Equation (5) yields the Klimontovich equation given as

d f d t = ∂ f ∂ t + v ⋅ ∂ f ∂ x + q m ( E + v × B ) ⋅ ∂ f ∂ v = C ( f ) (6)

This is the plasma kinetic energy. The term C ( f ) is known as the Coulomb collision operator on the average distribution function f = f ( x , v , t ) . The new form of Equation (6) becomes

∂ f ∂ t + ∂ ∂ x ⋅ v f + ∂ ∂ v ⋅ [ q m ( E + v × B ) ] f = C ( f ) (7)

For a plasma of particles, Coulomb collision effects over short time scales are negligible. Thus, for such plasma processes,

∂ f ∂ t + ∂ ∂ x ⋅ v f + ∂ ∂ v ⋅ [ q m ( E + v × B ) ] f = 0 (8)

The solution to this equation is

f = f ( c 1 , c 2 , c 3 , ⋯ , c N ) (9)

where c i is the constant of motion. The distribution is kinetically stable when the energy ϵ is the constant of motion and the equilibrium constant of distribution depends on it. Fluid moments are obtained by integrating low order powers of the product of velocity and the distribution function over velocity space in the laboratory frame, that is

∫ d 3 v v i ⋅ f : i = 0 , 1 , 2 , (10)

Theoretical Formulation

For a material in a superconducting state, two particles arise: normal particles that consist of normal electrons and the super-particles made up of super-electrons. The motion of the plasma super-particles around a magnetic vortex core contains the magnetic flux within the vortex core.

Some of the velocity moments of the distribution function f are listed below

Density of states:

N = ∫ d 3 v f (11)

Flow velocity:

V = 1 N ∫ d 3 v v ⋅ f (12)

Average Kinetic energy (meV):

E = N ( ∫ d 3 v m v r 2 3 f ) = m v r 2 2 (13)

Pressure (NM^{−2}):

p = ∫ d 3 v m v r 2 3 f = N E (14)

Pressure Tensor (NM^{−2}):

P = ∫ ( d 3 v ) m v r v r f = p Ι + π = p + π (15)

A system of interacting super-particles consists of a number (n) of electrons each carrying a charge e. The total charge of the super-particle will be denoted by . Physical properties of plasma such as energy can be generated from the velocity moment of the plasma kinetic equation of the super-particles that is expressed as

∫ d 3 v g ( v ) { ∂ f ∂ t + ∂ ∂ x ⋅ v f + ∂ ∂ v ⋅ [ Q m ( E + v × B ) ] f } = 0 (16)

where, g ( v ) is the velocity function for the particles. To find the density moment, g = 1 is used while for the momentum moment, g = m v is used as a velocity function. The momentum equation for the plasma super-particles is

∫ d 3 v m v { ∂ f ∂ t + ∂ ∂ x ⋅ v f + ∂ ∂ v ⋅ [ Q m ( E + v × B ) ] f } = 0 (17)

Each part is tackled separately using the commutation property before combining the solutions:

∫ d 3 v m v ∂ f ∂ t = ∂ ∂ t ∫ d 3 v f m v = m N d v d t (18)

∫ d 3 v m v ∂ ∂ x ⋅ v f = ∂ ∂ x ∫ d 3 v f m v v = ∇ ⋅ ( p + π ) (19)

∫ d 3 v m v ∂ ∂ v ⋅ [ q m ( E + v × B ) ] f = d d v ∫ d 3 v f [ m v Q m ( E + v × B ) ] = − [ N q ( E + v × B ) ] (20)

Substituting the solutions back to Equation (17) we have

m N ∂ v ∂ t = N Q ( E + v × B ) − ∇ ⋅ p − ∇ ⋅ π (21)

To determine the perpendicular guiding center (radius) drifts, the perpendicular flow responses are obtained by taking the cross product of the momentum equation with the magnetic field B. Thus, Equation (21) becomes

m N ∂ ∂ t ( B × v ) = N Q ( B × E + B × v × B ) − B × ∇ ⋅ p (22)

Performing the product yields

m N ∂ ∂ t ( B × v ) − N Q ( B × E ) + N Q ( B ⋅ V ) B + B × ∇ ⋅ p = 0 (23)

Dividing Equation (23) by N q B 2 we have

m Q B 2 ∂ ∂ t ( B × v ) + 1 B 2 ( B × E ) + 1 B 2 ( B ⋅ V ) B + 1 N Q B 2 ( B × ∇ ⋅ p ) = 0 (24)

Equation (24) gives the total flow velocity of the system

V = V ‖ + ϵ V ∧ + ϵ 2 V ⊥ (25)

The ϵ indicates the ordering of the various flow components while

V ‖ = V ‖ b ^ = 1 B 2 ( B ⋅ V ) B (26)

V ∧ = 1 B 2 ( E × B ) + 1 N Q B 2 ( B × ∇ ⋅ p ) (27)

V ⊥ = m q B 2 ∂ ∂ t ( B × v ) (28)

In this case, V ∧ is the cross flow velocity, V ‖ is the parallel flow velocity and V ⊥ is the perpendicular flow velocity. The directions, here, are relative to the direction of the magnetic flux, which is in the z-direction. Of importance, here, is the is the cross flow velocity given as

V ∧ = 1 B 2 ( E × B ) + 1 N Q B 2 ( B × ∇ p ) (29)

where,

V d = 1 N Q B 2 ( B × ∇ p ) (30)

V E = 1 B 2 ( E × B ) (31)

V d represents the diamagnetic flow velocity while V E represents the interaction between the fields. This study revolves around the diamagnetic flow velocity, V d that describes the motion of super-electrons while exerting pressure p on the magnetic flux. This pressure contains the flux along cylindrical vortices. The quantity B is the diamagnetisation of the super-particles. To work out the various components of the ∇ B tensor, a local Cartesian coordinate system with coordinates e ^ x , e ^ y and e ^ z is used such that e ^ z is aligned along b ^ at b ^ and e ^ x pointing towards the centre of the vortex core and perpendicular to b ^ . The y-axis of the local co-ordinate, e ^ y = b ^ × e ^ x becomes tangential to the surface of the vortex core (interface between the flux and the plasma-su- per-particles’ flow). Thus, the unit vectors characterizing this local coordinate system will be e ^ x , e ^ z = b ^ and e ^ y = b ^ × e ^ x . Equation (31) can be evaluated as follows:

V d = 1 N Q B 2 ( B × ∇ p ) (32)

Note that ∇ p is one-dimensional in the e ^ x -direction i.e. towards the centre of the cylindrical vortex.

1 N Q B 2 ( B × ∇ p ) = 1 N Q B | e ^ x e ^ y e ^ z 0 0 B d p d x 0 0 | = 1 N Q B ( d p d x ) e ^ y (33)

Recall that p ≅ N E , so that E p = 1 , where E is the kinetic energy of the super-particles, at absolute zero temperature and E_{k} is the excitation energy of the quasi particles. Therefore, equation becomes

V d = E ( Q e ) B ( 1 p d p d x ) e ^ y (34)

The quantity 1 L p = − 1 p d p d x is known as the pressure-gradient scale length and is proportional to the radius, r of the vortex core, i.e.

r = 1 L p = − 1 p d p d x (35)

The vector e ^ y shows that plasma superfluid moves tangentially to the interface as the flux concentrates about the centre of the vortex. These interactions form a cylindrical vortex core of radius r = 1 L p at diamagnetisation B and flux H. The diamagnetisation B is proportional to the μ 0 ( H c 2 + H 0 H ) e − ( H c 1 H ) . Therefore,

B = η μ 0 ( H c 2 + H 0 H ) e − ( H c 1 H ) (36)

where η a constant that depends on the material of the superconductor is, H is the applied external field, H c 1 is the lower critical field, H c 2 is the upper critical field, μ 0 is the permeability of free space and H H c 2 = b is the reduced magnetic field [

V d = − 1 η μ 0 ( H c 2 + H 0 H ) ( E ( Q e ) r ) e ( H c 1 H ) (37)

When the external magnetic field H is increased, the resulting magnetic flux (ф) mounts more pressure on the superelectrons and weakens the diamagnetisation. This results into an increase in the radius, r, of the core. At H = H c 2 , the supe-particles are so weak that the magnetic flux permeates the whole spectrum of the super-particles causing the superconductivity of the material to breakdown and r → ∞ . It is important to note that the values of E are in meV. Consequently,

V d = 1 η ( H c 2 + H 0 H ) ( E μ 0 ( Q e ) r ) e − ( H c 1 H ) (38)

At absolute zero temperature, the thermal velocity of the super electrons is zero and therefore, Equation (38) relies on the acoustic velocity. Hence, for a super-electron at absolute zero temperature,

V d = E m p (39)

Hence, using Equations (39) into (38) we get

Q = − 1 η ( H c 2 + H 0 H ) ( E ( E m p ) μ 0 r e − ( H c 1 H ) ) e (40)

From Equation (40), it should be noted that for a three-electron system,

μ 0 E 3 m e = 12.57 × 10 − 7 E 27.3 × 10 − 31 = ( 7.608 × 10 8 ) E (41)

Substituting 41 into 40, we have

Q = − 1 η ( H c 2 + H 0 H ) ( E r e ( H c 1 H ) × 1.314 × 10 − 9 ) e (42)

If we make an assumption that the boundaries of the vortex core are clear when B = B 0 and H = H 0 . Therefore, the radius r = r 0 is experimentally determined at this point and

B = B 0 = η μ 0 ( H c 2 + H 0 H 0 ) e − ( H c 1 H 0 ) (43)

At B = B 0 = μ 0 e − ( H c 1 H 0 ) , r = r 0 , H = H 0 and η = ( H 0 H c 2 + H 0 ) . Consequently, at B = B 0 , Equation (42) becomes

Q e = − ( 1 r 0 ) ( 1.728 E × 10 − 18 ) 1 2 e ( H c 1 H 0 ) (44)

The average kinetic energy ( E ) of super-particle consisting of three super-electrons in a YBCO123 system, by calculation, is found to be 24 meV while that of Bi2212 is 25.1 meV. Msass of the model m p = 3 m e . Equations (36) and (43) has been used to generate the graphs of diamagnetisation, B as a function of the external applied magnetic field H .

In

The results from

( ( 1 r 0 ) ( 1.728 E × 10 − 18 ) 1 2 e ( H c 1 H 0 ) ) ≈ 3 (45)

Superconductor | r 0 ( Å ) | ( 1 r 0 ) ( 1.728 E × 10 − 18 ) 1 2 e ( H c 1 H 0 ) |
---|---|---|

YBCO123 | 20 [ | 3.2 |

Bi2212 | 22 [ | 3.1 |

At H = H 0 ≫ H c 1 , e − ( H c 1 H 0 ) ≅ 1 . Thus, r 0 can be estimated from the table as

r 0 = 1 3 ( 1.728 E × 10 − 18 ) 1 2 Å (46)

Equation (46) shows the dependence of the vortex core radius, r 0 , on the kinetic energy and the diamagnetisation of the super-electrons. Increase in the kinetic energy increases the centrifugal force on the super-particles and a consequent increase in the radius of the vortex core.

It has been shown through the working that the radius of the magnetic vortex core, when the boundary is clear, increases with the kinetic energy ( E ) of the super-particles.

The authors declare no conflicts of interest regarding the publication of this paper.

Mukubwa, A. and Masinde, F. (2019) Determining the Radius of the Magnetic Vortex Core of YBCO123 and Bi2212. Open Access Library Journal, 6: e5661. https://doi.org/10.4236/oalib.1105661