Expulsion from an Infinite Cylindrical Superconductor

Miguel C. N. Fiolhais∗ LIP, Department of Physics, University of Coimbra, 3004-516 Coimbra, Portugal

Hanno Ess´en† Department of Mechanics, KTH, 10044 Stockholm, Sweden (Dated: July 15, 2014) Abstract The solutions of the London equations for the magnetic field expulsion from superconductors are presented in this paper for the cylindrical symmetry. The result is analyzed in detail and represented numerically for the case of a uniform external magnetic field in the transverse plane. In particular, several contour plots of the magnetic energy density are depicted for the regions inside and around the superconducting area for a wide range of penetration lengths, showing how the expulsion and penetration of the magnetic field evolve with the ratio between the penetration length and the cylinder radius.

The following paper is published in Physica C: http://www.sciencedirect.com/science/article/pii/S0921453413004541

INTRODUCTION (∇ · A = 0), one obtains:

1 ∇2A − A = 0 , (2) The first phenomenological description of the magnetic λ2 field expulsion from superconductors appeared soon af- ter the discovery of the Meissner-Ochsenfeld effect [1], by or, in terms of the magnetic field, the London brothers in 1935 [2]. Ever since, several dis- 1 cussions have been presented along the years on the na- 2 ∇ B − 2 B = 0 . (3) ture of the magnetic field expulsion configuration in the λ scientific literature [3–17]. These studies comprise both This equation is perhaps the most common form of the static and time-dependent analyses, several different ge- London equation for the magnetic field expulsion inside ometries, and even the case of non-homogenous super- the superconducting region, and, in the cartesian coordi- conductors. However, while the spherical case is studied nate system, corresponds to the well-known Helmholtz’s in [17–19], and the cylinder case in a parallel magnetic equation for a complex wavenumber, for each of the vec- field is presented in [19], to the best of our knowledge tor components. The magnetic field is, therefore, sup- the detailed solution of a superconducting cylinder in a pressed in the interior region through negative exponen- transversal magnetic field cannot be easily found in the tial dependencies, or, in the cylindrical and spherical academic literature. This very same case has been pre- symmetries, through modified Bessel functions of the first viously addressed by Zhilichev [7] by proposing a macro- kind [20]. scopic shell model, but without providing an explicit so- lution for the London condition. To fill this gap, we de- rive the solution of the London equations in this paper, CYLINDRICAL SYMMETRY for the magnetic field expulsion from an infinite homo- geneous superconducting cylinder in a constant external Consider a cylindrical superconductor with radius R magnetic field in the transverse plane. The dependence of in an external constant perpendicular magnetic field, in the strength and configuration of the expelled magnetic such a way the external magnetic field points in the y- field with the London penetration depth is examined in direction and the cylinder axis coincides with the z-axis. detail. In order to obtain the final static vector potential configu- The London equations for the magnetic field flux ex- ration resulting from the magnetic flux expulsion, and as- pulsion out of a superconductor can be expressed, in suming the Coulomb gauge, one must solve the Maxwell arXiv:1407.3393v1 [cond-mat.supr-con] 12 Jul 2014 terms of the magnetic vector potential A, in a single equation for the outer region: equation: ∇ × B = ∇ × (∇ × A) = ∇2A = 0 , (4) 1 and the London equation (2) for the inner region. Due j = − 2 A , (1) λ to the symmetry of the system, the magnetic field has no component in the z-axis direction, where j represents the and λ the London penetration depth. After applying the Maxwell’s equa- 1 ∂A ∂A B = z ρˆ − z φˆ , (5) tion for the static case, and by taking the Coulomb gauge ρ ∂φ ∂ρ

1 and, therefore, the axial component of the vector poten- respectively. As a result, one can establish the following tial Az is the only relevant component needed to compute condition: the magnetic field. For the outer region, this component ∂A ∂A can be determined from equation (4): z = z , (14) ∂ρ+ ∂ρ−  1 ∂2A ∂2A 1 ∂A  z + z + z zˆ = 0 , (6) which, with the use of the relation, ρ2 ∂φ2 ∂ρ2 ρ ∂ρ

∂Iν (x) ν which is simply the Laplace equation for Az. The most = Iν−1(x) − Iν (x) (15) general solution in the outer region is, therefore, ∂x x in the differentiation process, leads to, P+∞ Az = n=0 [An cos (nφ) + Bn sin (nφ)]  n −n  2  × Cnρ + Dnρ . (7) R R I1(R/λ) Bρ>R = B0 1 − 2 + 2λ 2 sin φ ρˆ ρ ρ I0(R/λ) Assuming the magnetic field takes the form of the exter-  2  R R I1(R/λ) ˆ nal field in the limit ρ → +∞, + B0 1 + 2 − 2λ 2 cos φ φ , ρ ρ I0(R/λ)   λ I (ρ/λ) B0 = B0ˆy = B0 φˆcos φ + ρˆsin φ , (8) 1 Bρ

Az = E0G0I0(ρ/λ) + E1G1I1(ρ/λ) cos φ , (11) shrinks to an ellipsoid as the London depth increases and the magnetic field penetrates more and more. The mag- where, netic field energy density reaches its maximum value at the cylinder’s border in the horizontal axis due to the E1G1I1(R/λ) = −B0R + A1D1/R . (12) concentration of field lines in this region, and the mini- mal region is located at the top and bottom of the cylin- Since the electric current is distributed in volume, der. As the penetration depth increases, the maximal j = − 1 A, there are no pure surface currents at the cylin- λ2 regions tend to fade while entering the cylinder to form drical surface, allowing some penetration of the magnetic the ellipsoidal shaped region. field. Therefore, the following relation can be used as a Furthermore, one can see that no dramatic changes boundary condition: are perceptible until the penetration depth becomes one k = ˆn × B+ − B− = 0 , (13) order of magnitude lower than the radius of the cylin- der, when the ellipsoid starts to form. It is also worth where k is the surface current, ˆn is the unitary vector noting that when the London depth equals the radius of normal to the surface and, B+ and B− are the magnetic the cylinder the magnetic field expulsion is almost imper- field at the cylinder surface in the outer and inner limits, ceptible, or in other words, the superconductor becomes

2 transparent to the magnetic field. This is also visible in Figure 2, where the dependence of the magnetic energy density along the x-axis is shown at y = 0 for different penetration depths. In particular, a single jump on the order of magnitude, from λ = 0.1 R to λ = R, is enough to have a significant phenomenological impact, i.e. the difference between an almost fully expelled field and a considerable field penetration. Finally, it is also worth discussing the range of ap- plicability of these solutions in both type-I and type-II superconductors. In type-I superconductors, the London penetration length is of the same order of magnitude or smaller than the coherence length [21, 22]. Therefore, the most accurate description is not given by the London equations alone, but by the Pippard’s model instead [23]. However, in type-II superconductors, the London condi- tion indeed provides a very good description of the elec- tromagnetic field inside the superconductor in the Meiss- ner state, i.e. below the lower critical field. Above the lower critical field, the superconductor enters a mixed state, where the magnetic field is allowed to penetrate the superconducting region through quantized vortices [24], and the London solutions are no longer an accurate de- scription.

CONCLUSIONS

We presented a detailed solution of the London equa- FIG. 1: Contour plots of the magnetic energy density for λ = 0 (top left), λ = 0.03 R (top right), λ = 0.1 R (center tion for the magnetic field expulsion from an infinite left), λ = 0.2 R (center right), λ = 0.4 R (bottom left), and cylindrical superconductor, which can be easily obtained λ = R (bottom right). by using physical and differentiability arguments at the boundary conditions. The field expulsion is analyzed with the use of contour plots of the magnetic energy den-

2 0 4.5 sity to show how it evolves with the value of the penetra- λ = 0 / B 2 4 λ = 0.1 R tion length. In this case, the magnetic field penetrates B λ = 0.4 R λ = R through the formation of an elliptic region which becomes 3.5 more pronounced as the penetration depth increases. 3 2.5

2

1.5 ∗ Electronic address: miguel.fi[email protected] 1 † Electronic address: [email protected] 0.5 [1] Meissner W. and Ochsenfeld R., Naturwiss, 21, 787 0 •3 •2 •1 0 1 2 3 (1933) x / R [2] London F. and London H., Proc. Roy. Soc., A149 (866), 71 (1935) [3] Cullwick E. G., Proc. IEE, Part C, Monographs, 103, 441 (1956) FIG. 2: Dependence of the magnetic energy density with the [4] P. G. de Gennes, of Metals and Alloys, x-coordinate at y = 0 for different penetration lengths. Perseus books, Reading (1999) [5] Pfleiderer J., Nature, 210, 614 (1966) [6] Karlsson P. W., Arch. f. Elektrotech., 67, 29 (1984) [7] Zhilichev Y. N., IEEE Trans. Appl. Supercond., Vol. 7, [9] Ejrnaes M. and Nappi C., Supercond. Sci. Technol., 16, 3874 (1997) 1224 (2003) [8] Prozorov R., Giannetta R. W., Carrington A. and [10] Bad´ıa-Maj´osA., Am. J. Phys., 74, 1136 (2006) Araujo-Moreira, F. M., Phys. Rev. B, 62, 115 (2000) [11] Kudinov E. K., Phys. Solid State, 48, 1147 (2006)

3 [12] Fiolhais M. C. N., Ess´enH., ProvidˆenciaC. and Nord- [19] Batygin V. V. and Toptygin I. N., Problems in Electro- mark A. B., Prog. Electromagn. Res. B, 27, 187 (2011) dynamics, 2nd ed., Academic, London (1978) [13] Ess´enH., Europhysics Letters, 94, 47003-1 (2011) [20] M. Abramowitz and I. A. Stegun, Handbook of Math- [14] Ess´enH. and Fiolhais M. C. N., Am. J. Phys., 80, 164 ematical Functions with Formulas, Graphs, and Mathe- (2012) matical Tables, Dover, New York (1965) [15] Fiolhais M. C. N. and Ess´enH., Int. J. Theor. Phys., 52, [21] Ginzburg V. L. and Landau L. D., Zh. Eksp. Teor. Fiz., 1701 (2013) 20, 1064 (1950) [16] Fiolhais M. C. N. and Kleinert H., Phys. Lett. A, 377, [22] Fiolhais M. C. N. and Kleinert H., Phys. Lett. A, 277, 2195 (2013) 2195 (2013) [17] Matute E. A., Am. J. Phys., 67, 786 (1999) [23] Pippard A. B., Proc. R. Soc. Lond. A, 203, 210 (1950) [18] Reitz J. R., Milford F. J., and Christy R. W., Foun- [24] Abrikosov A. A., Zh. Eksp. Teor. Fiz., 32, 1442 (1957) dations of Electromagnetic Theory, 4th ed., Addison- Wesley, Reading (1993)

4