Abstract

Cubic superconducting sample was simulated using time-dependent Ginzburg-Landau model under oscillating magnetic field with and without additional background static magnetic field. Vortex dynamics including entrance and exit from the sample was simulated. Magnetization and carrier concentration densities of the sample were studied as a function of external magnetic field variations. Anomalies in carrier concentration density were observed at certain values of the magnetic field which were correlated with the entrance and exit processes of vortices. Area swept by superconductor magnetization with magnetic field was observed to have a hysteresis-like behavior where area representing energy dissipated per cycle. This energy accumulation was suggested to cause instability in superconductor over the number of cycles and may result in thermal quenching. Temporal distribution of energy components showed consistency with the pattern observed for carrier concentration and magnetization under oscillating magnetic field. Rapid phase changes with magnetic oscillations resulted in oscillations in energy components, and irregular peaks and ripples in superconducting energy represent the situation of exit and entry of vortices. While the rise in interaction energy with cycles is referred to vortex relaxation time in a cycle, this energy is expected to accumulate and take other forms (e.g., heat) and is predicted to cause thermal quenching. In the presence of background static magnetic field, this energy dissipation was calculated to increase significantly while superconductor is subjected to oscillating magnetic field.

1. Introduction

In general, the behavior of superconductor in applied electric and magnetic field is attributed to the vortex dynamics in it, which generally results in energy dissipation (e.g., heat) and resulting thermal quenching of superconductivity. According to Bardeen-Cooper-Schrieffer (BCS) theory, spin singlet state (opposite spin electron pair), bound by phonon interaction, constitutes the ground state of condensate [1]. The proposed macroscopic quantum theory by Ginzburg and Landau (named Ginzburg-Landau (GL) theory) used Higgs mechanism of spontaneous symmetry breaking and quartic potential [25] to generate a local mass term of vector potential. Development of homogeneous ferromagnetic spins destroys the superconductivity if it exceeds critical magnetic field of the superconductor. This model successfully describes the Meissner-Ochsenfeld effect [4, 68]. The prediction of penetrating strong magnetic fields in type-II superconductors by Abrikosov [9] gave further credibility to GL model which was later verified experimentally [1012]. Over the past 70 years, not only macroscopic properties of superconductors (such as categorizing superconductors in type-I and type-II, and description of vortex state in type-II superconductors [9]) but also mesoscopic superconducting samples were successfully described using this theory. GL model is probably the most accurate phenomenological model to describe macroscopic properties of superconductors [13, 14]. Gorkov [15] proved that close to critical temperatures microscopic BCS theory also reduces to GL theory. Various numerical methods, including finite difference [1618], finite element [1921], and spectral method [22], have been developed for the solution of this model. In this report, we used the semi-implicit finite difference method with staggered grid scheme.

Generally, high-temperature superconductor (HTS) based magnetic Levitation (maglev) vehicle systems consist of HTS on a permanent magnet (PM) guideway [2325]. In general, uniform magnetic field is considered for PM guideways but practically this is not the case; there are situations of magnetic field variations (e.g., cracks, magnetic contacts, and structural and magnetic defects) resulting in a variable magnetic field. The frequency of this variation depends on the speed of the HTS on PM guideway (i.e., speed of the vehicle). The magnetic field in such a case is not exactly oscillating (which was reported earlier [18]), but it is (practically) a constant magnetic field having small jitters (or oscillations) in its amplitude appearing due to imperfections in PM guideway. In this work, we studied such a magnetic field and compared the results with results reported earlier for perfectly oscillating magnetic fields [18]. To investigate this type of energy loss, we studied the behavior of type-II superconductor exposed to an oscillating magnetic field simulated using time-dependent Ginzburg-Landau (TDGL) equations [15, 26, 27] near critical temperature with and without an additional static background magnetic field. The dynamics of Abrikosov vortices in oscillating magnetic field is studied. A hysteresis-like behavior of sample magnetization was observed under oscillating magnetic fields, while background static magnetic field plays a vital role in thermal quenching of superconductor. The behavior of different energy components is also discussed in detail.

2. Ginzburg-Landau Theory

Primary variables in GL model are order parameter, ψ, and magnetic vector potential, A, occupying a superconductor sample in a three-dimensional region Ω having boundary Γ. In nondimensional state, these primary variables are linked to physical quantities, given as follows:Density of superconducting charge carriers, Induced magnetic field, Current density, .

Ginzburg-Landau free energy is taken from Gorkov and Eliashberg [28] given as follows:where is charge and is effective mass of superconducting charge carriers, is external magnetic field, and are reduced Plank's constant (i.e., h/2π), and , while α and β are phenomenological parameters depending on environmental factors (e.g., temperature).

We introduce London penetration depth () and coherence length () as length scale parameters. We introduce a length scale and nondimensionalized physical quantities, special and time coordinates () and t, order parameter ψ, magnetic vector potential , and magnetic field according to the following transformations:where γ is relaxation parameter. For simplicity, we choose length scale in terms of coherence length (i.e. l = ξ). Ginzburg-Landau parameter κ is defined as . κ defines the type of superconductor, for type-I and for type-II superconductors. Minimizer of Ginzburg-Landau energy functional (1) satisfies Euler-Lagrange equations, generally known as Ginzburg-Landau (GL) equations:σ is electrical conductivity defined as , while is normal conductance. External magnetic field is assumed to be uniformly distributed over the sample surface. GL equations were solved using natural boundary conditions, i.e., supercurrent across the boundaries is zero and magnetic field at the boundaries is given as follows:where is unit normal vector on the boundary Γ of the region Ω.

TDGL model is gauge invariant under transformation:where gauge χ is a function of space and time. For the present study, we choose zero electric potential gauge () known as Coulomb Gauge. So (3) and (4) reduce toEnergy density Htot=Htot(x,y,z,t) in superconductor is distributed into three parts: superconducting energy density Hsup=Hsup(x,y,z,t), interaction energy density Hint=Hint(x,y,z,t), and magnetic energy density Hmag=Hmag(x,y,z,t) given as[18, 29] where length scales are taken in terms of London penetration depth (i.e. l=λ). Total energy (Htot) density is the sum of these three energy densities, Htot = Hsup + Hint + Hmag; total energy is the integral of total energy density over the region Ω, .

Gauge invariant discretization [17, 3033] is the most popular and widely used approach to solve TDGL equations, which is first-order accurate in time and second-order accurate in space, and other finite element [19, 20], finite difference [34, 35], and spectral method [22] have also been developed. Link variable schemes [16, 17, 36] are the most commonly used schemes to solve coupled GL equations. In the present work, in order to investigate carrier concentration, magnetization, and vortex dynamics in superconductor, we used coupled nonlinear TDGL equations (8) and (9) solved by finite difference scheme using link variables, described by Winiecki and Adams [16], implemented in three-dimensional code we have developed [18]. Equations (8) and (9) along with boundary condition equations (5) and (6) form the basis of the present work. The present work was performed close to critical temperature allowing the assumption of thermal suppression of surface barrier [3739]. The frequency of oscillating magnetic field was set to 0.00278/t in nondimensionalized coordinates for the present work.

3. Results and Discussion

Three-dimensional, cubic superconductor domain of size , periodic along z-axis (along the direction of applied magnetic field), was discretized with grid size 0.5 ξ in each direction, and magnetic field was applied in z-axis direction with Ginzburg-Landau parameter κ=4. Initially, a constant magnetic field was applied, and on reaching equilibrium state, it was allowed to oscillate as cosine wave. Cubic superconductor under two different types of magnetic field schemes, (A) and (B) (where θ=ωt=2πft), was simulated to investigate the effect of static background magnetic field on the superconductor in the presence of an oscillating magnetic field. The system was initialized with the completely superconducting sample () and no magnetic vector potential (A=0). The system was relaxed to minimize Gibbs free energy. During this energy minimization process, vortices entered the sample and relaxed their positions to minimum energy. Magnetic field oscillations were introduced as the sample reached minimum energy state, resulting in a new series of entrance and leaving of vortices, corresponding to different positions of magnetic field wave. Figure 1 shows a variation of average carrier concentration of the sample with magnetic field. The pattern of carrier concentration is symmetric for positive and negative half cycles but nonsymmetric for increasing and decreasing parts of magnetic field oscillation, shown in Figure 1 (for magnetic field of type (B)), indicating a lag in vortex dynamics behind the oscillating magnetic field. For the situation of magnetic field of type (A), the system is initially subjected to , and after sufficient time, on reaching equilibrium state, a small oscillation with amplitude 0.15 in magnetic field was introduced.

During decaying half of magnetic field cycle (i.e., ), carrier concentration does not increase as rapidly as decrease in magnetic field, till a certain value of magnetic field, where it is difficult for the sample to keep the vortices bound any further within the sample and a set of loosely bound vortices leave the sample, causing a sudden increase in carrier concentration shown by steps between positions (b) and (c) in Figure 1, carrier concentration peaks at minimum magnetic field, i.e., . During increasing half of magnetic field cycle (i.e., ), penetration of magnetic field in large quantity starts from the sample edges till vortices are formed and they repel magnetic field penetrations from the edges. Relaxation of vortices causes an increase in carrier concentration indicated by sections (e) and (f) in Figure 1, and formation and relaxation of the first additional set of vortices are indicated by this small peak; the same process repeats itself for every subsequent set of vortices entering the sample indicated by (g) in Figure 1.

Energy dissipation () in magnetic oscillations is the area covered by hysteresis loop (M-H curve) of sample magnetization with magnetic field, which takes the following form:where is sample magnetization (), is volume bounded by the superconductor, and is the applied external magnetic field. This type of behavior of magnetic field trapped inside the superconductor lagging behind the magnetic field oscillations indicates the energy dissipation during the cycle. The area of loop, attributed to the energy dissipated during a cycle, was observed to increase with frequency of magnetic oscillations, due to limited velocities of vortices. Figure 2 shows variation of average sample magnetization () with external magnetic field oscillations at normalized frequency f=0.00278/t. It was observed that during decreasing and increasing half cycles of magnetic field oscillations, sample magnetization behaved in a similar way (comparing with plots behavior) to that of carrier concentration in Figure 1. During decreasing part of magnetic field sample magnetization had a rapid drop between (b) and (c) in Figure 2 at the same magnetic field intensity as observed between (b) and (c) in Figure 1. This behavior of a sudden drop in magnetization, between (b) and (c) in Figure 2, and sudden rise in carrier concentration, between (b) and (c) in Figure 1, represent the same process and are attributed to the process of a set of vortices leaving the sample. This trapped energy can be anticipated to transform in other forms of energy, mostly heat, which accumulates with a number of cycles, resulting in thermal quenching of the superconductor. It can also be inferred that, at faster vehicle speeds, resulting in higher frequencies, this energy accumulation may result in a rapid quenching of superconductor. This energy loss during a cycle was observed to be ~10 times greater in the presence of additional static background magnetic field (type A) compared with that in its absence (type B), which can be seen by the area swept by both curves in Figure 2.

Temporal distribution of three energy components and total energy for the sample under the magnetic field of type (A) is shown in Figure 3. Initially, (at t=0) sample was in a nonequilibrium state with no magnetic field penetration and carrier concentration at maximum. Nonequilibrium state of magnetic field distribution across the boundary is the reason of large initial value of Hmag, but as the field penetration starts, Hmag rapidly drops to equilibrium values, because the number of vortices will remain the same before introduction of field oscillations and there would only be rearrangement of vortices in the sample, and Hmag do not change during the rearrangement of vortices. Initially, due to uniform penetration of vortices from the boundaries, Hsup increases rapidly because average carrier concentration decreases and the second term in (10) dominates. As the vortices are formed and start relaxing, Hsup starts decreasing forming an initial hump-like structure, indicating nonvortex to vortex transition. Hint is responsible for the interaction of magnetic and superconducting energies. The second term in (11) is responsible for the formation of vortices while the first term controls vortex dynamics inside the superconductor.

Figure 3 (region II) represents all the three component energy densities and total energy density with the magnetic field of type (A). As system reaches equilibrium, oscillations in magnetic field are introduced as cosine wave, in the range 0.6κ ≥ 0.3κ. In region II as decreases, the balance between external magnetic field and vortex magnetization is disturbed, resulting in an increase in intervortex spacing and exit from the sample, and all the energies start oscillating with the oscillations in , and the process of repeated entrance and exit of vortices starts. Peak shape analysis of Hsup and Hint shows small kinks in energies, indicating the moments when outermost (or least bound) vortices leave the sample; after that, a sudden dip was observed in Hsup and Hint because the remaining vortices which are now free from force exerted by the vortices just left and relaxed themselves; this kink is also visible in total energy pattern, similar to the one observed at positions (b) and (c) in Figures 1 and 2. Interestingly, energy patterns of Hsup and Hint are not symmetric for increasing and decreasing half cycles of , i.e., entrance and exit of vortices. Peak shape analysis shows a quick entrance and slow exit of vortices during increasing and decreasing half cycles of , resulting in hysteresis-type behavior and energy accumulation which have already been discussed. An irregularity in Hint is observed after a number of cycles, due to energy accumulation. This effect was only observed in Hint due to its dependence on vortex dynamics, having limited velocities, resulting in Hint reaching equilibrium value with a delay compared to Hmag and Hsup. Repeated entrance and leaving of vortices over a cycle do not give sufficient time to the sample to completely attain equilibrium state, resulting in accumulation of energy in the sample in the form of interaction energy. This accumulated energy represented the area enclosed by magnetic field and sample magnetization shown in Figure 2, which may appear in the form of heat resulting in thermal quenching of superconductor in practical applications. This effect is significantly high (~10 times) for the samples under magnetic field of type (A) compared to the samples under magnetic field of type (B). Figure 4 shows the comparison of Hint of the systems with and without steady magnetic field (of types (A) and (B)) under the same oscillating component of magnetic field. It was observed (by the comparison of interaction energies in Figure 4) that, in the presence of steady background magnetic field, the energy dissipation (or the abnormalities in Hint) is ~10 times higher than that observed for the case without additional external magnetic field. Therefore, a static background magnetic field plays a significant role in addition to oscillating magnetic field in thermal dissipation and subsequently quenching of superconductor moving on PM guideway in a Maglev vehicle system; the rate of such energy dissipation/accumulation is directly dependent on the speed of the vehicle.

4. Conclusions

Time-dependent Ginzburg-Landau model was used to investigate type-II superconductor under the oscillating magnetic field, with additional static background magnetic field. Finite difference scheme was used for cubic superconducting sample. Carrier concentration, magnetization, and energy in the presence and absence of static background magnetic field were studied. Nonuniform behavior in carrier concentration with the oscillating magnetic field was observed and was related to vortex dynamics on point-to-point basis. Lagging of average sample magnetization behind magnetic field was observed during magnetic field oscillations. Area of hysteresis loop was measured and was attributed to dissipated energy in the sample, accumulated over time. The same effect caused anomalies in interaction energy density of the sample under magnetic field oscillations. Ripples in superconducting energy density were observed due to entrance and exit of vortices during magnetic field oscillations, which is the result of the same effect of small peaks (and sudden rise) observed in average carrier concentration with increasing (and decreasing) external magnetic field magnitude. Comparison between energy patterns in the presence and absence of background steady magnetic field indicated that more energy is dissipated in the presence of background magnetic field than that in its absence, suggesting that, in practical devices (e.g., Maglev vehicles) when magnetic field is oscillating (magnetic guideways with defects), its necessary to shield the superconductor from background magnetic fields to avoid energy dissipation and thermal quenching of the superconductor.

Data Availability

All the results published in this article can be reproduced by using equations and parameters described herein.

Conflicts of Interest

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

Acknowledgments

This work was supported by National Science Foundation of China, Grant nos. 11174030 and 11504020. Author Hasnain Mehdi Jafri acknowledges the Higher Education Commission of Pakistan for the Ph.D. scholarship.

Supplementary Materials

Figure S1: vortex structure of the sample at (a) Bz=0.6κ (oscillating magnetic field), (b) Bz=0.3κ (oscillating magnetic field), (c) Bz= 0.15κ (oscillating magnetic field), and (d) Bz= 0.15κ (nonoscillating magnetic field) for cubic sample of size 20ξ×20ξ×20ξ, with parameters κ=4 and σ=1. (Supplementary Materials)