Magnetic Field 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 electric current 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ρ