Determining the Radius of the Magnetic Vortex Core of YBCO123 and Bi2212 ()
1. Introduction
Ginzburg-Landau (GL) theory gives a good account of interactions between magnetic field and electric field in a type II superconductor [1] . The theory agrees well with the BCS theory. However, the BCS theory is not sufficient to explain the high energy pairing in the High-Tc 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 [2] . The plasma of charged particles responds inertially to electric field perturbations by oscillating at the electron plasma frequency. Externally imposed electric fields will induce perturbations in the plasma that are combinations of the time dependence of the externally imposed field and the electron plasma oscillations. The response of charged particles to the electric field force is limited by the inertial force
, and is inversely proportional to the mass of the charged particle. Thus, the lighter electrons will give the primary inertial response to an electric field perturbation in a plasma [2] . Consequently, the electrons give the dominant contribution to the plasma polarization. In the Langmuir waves, plasma electrons oscillate against a stationary ion background in a normal state of a material [3] . The force experienced by an ith particle in motion is given by
(1)
where
is the velocity of the ith particle. The equation of motion of a charged particle in a magnetic field (
) due to an electric field (
) is given by [4] ,
(2)
where
is the velocity of the ith particle, q is the particle charge and
is the Lorentz force on the ith particle. The charge and current densities respectively are obtained by integrating the distribution function over velocity space and summing over species
(3)
And
(4)
where
is the microscopic distribution function whose time derivative is given by
(5)
here,
Putting Equation (2) into Equation (5) yields the Klimontovich equation given as
(6)
This is the plasma kinetic energy. The term
is known as the Coulomb collision operator on the average distribution function
. The new form of Equation (6) becomes
(7)
For a plasma of particles, Coulomb collision effects over short time scales are negligible. Thus, for such plasma processes,
(8)
The solution to this equation is
(9)
where
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
(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:
(11)
Flow velocity:
(12)
Average Kinetic energy (meV):
(13)
Pressure (NM−2):
(14)
Pressure Tensor (NM−2):
(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
(16)
where,
is the velocity function for the particles. To find the density moment,
is used while for the momentum moment,
is used as a velocity function. The momentum equation for the plasma super-particles is
(17)
Each part is tackled separately using the commutation property before combining the solutions:
(18)
(19)
(20)
Substituting the solutions back to Equation (17) we have
(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
(22)
Performing the product yields
(23)
Dividing Equation (23) by
we have
(24)
Equation (24) gives the total flow velocity of the system
(25)
The
indicates the ordering of the various flow components while
(26)
(27)
(28)
In this case,
is the cross flow velocity,
is the parallel flow velocity and
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
(29)
where,
(30)
(31)
represents the diamagnetic flow velocity while
represents the interaction between the fields. This study revolves around the diamagnetic flow velocity,
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
is the diamagnetisation of the super-particles. To work out the various components of the
tensor, a local Cartesian coordinate system with coordinates
,
and
is used such that
is aligned along
at
and
pointing towards the centre of the vortex core and perpendicular to
. The y-axis of the local co-ordinate,
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
,
and
. Equation (31) can be evaluated as follows:
(32)
Note that
is one-dimensional in the
-direction i.e. towards the centre of the cylindrical vortex.
(33)
Recall that
, so that
, where
is the kinetic energy of the super-particles, at absolute zero temperature and Ek is the excitation energy of the quasi particles. Therefore, equation becomes
(34)
The quantity
is known as the pressure-gradient scale length and is proportional to the radius, r of the vortex core, i.e.
(35)
The vector
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
at diamagnetisation B and flux H. The diamagnetisation B is proportional to the
. Therefore,
(36)
where
a constant that depends on the material of the superconductor is, H is the applied external field,
is the lower critical field,
is the upper critical field,
is the permeability of free space and
is the reduced magnetic field [5] . Considering all these cases, Equation (34) becomes
(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
, 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
. It is important to note that the values of
are in meV. Consequently,
(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,
(39)
Hence, using Equations (39) into (38) we get
(40)
From Equation (40), it should be noted that for a three-electron system,
(41)
Substituting 41 into 40, we have
(42)
If we make an assumption that the boundaries of the vortex core are clear when
and
. Therefore, the radius
is experimentally determined at this point and
(43)
At
,
,
and
. Consequently, at
, Equation (42) becomes
(44)
2. Results and Discussion
The average kinetic energy (
) 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
. Equations (36) and (43) has been used to generate the graphs of diamagnetisation,
as a function of the external applied magnetic field
. Figure 1 shows a graph of diamagnetisation,
as a function of the external applied magnetic field
over varied ranges of the applied external magnetic field for both YBCO123 and Bi2212 systems.
In Figure 1(a), the Continuous line represents the theoretical results while the dotted line represents the experimental results. Similar shapes have been obtained through experimental procedures for Type II superconductors. In Figure 1(a), the peak for both the experimental and theoretical occur at the lower critical field. Analysis of Figure 1(b) shows some interesting feature that emanates from the assumption in Equation (43). At
, the theoretical estimation puts the external applied field, at which the vortex core boundaries are clear, at 6T which agrees fully with the experimental procedure that gives 6T. The curve finally ends while it is very close the x-axis where the upper critical field of 120T. Table 1 shows the estimated values of the ratio
from Equation (44) based on experimental values of
.
The results from Table 1 have shown that
(45)
(a)(b)(c)
Figure 1. Graph of diamagnetisation as a function of applied external magnetic field over the range (a)
(b)
(c)
.
Table 1. Estimated values of
.
At
,
. Thus,
can be estimated from the table as
(46)
Equation (46) shows the dependence of the vortex core radius,
, 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.
3. Conclusion
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 (
) of the super-particles.