Arxiv:1407.3393V1 [Cond-Mat.Supr-Con] 12 Jul 2014
Total Page:16
File Type:pdf, Size:1020Kb
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´eny 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 (r · A = 0), one obtains: 1 r2A − 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 r 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: r × B = r × (r × A) = r2A = 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+1 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 ρ ! +1, + B0 1 + 2 − 2λ 2 cos φ φ ; ρ ρ I0(R/λ) λ I (ρ/λ) B0 = B0^y = B0 φ^cos φ + ρ^sin φ ; (8) 1 Bρ<R = 2B0 sin φ ρ^ ρ I0(R/λ) and in order to respect this boundary condition, the vec- I0(ρ/λ) λ I1(ρ/λ) ^ tor potential becomes, + 2B0 − cos φ φ : (16) I0(R/λ) ρ I0(R/λ) −1 Az = const. + A1 C1ρ + D1ρ cos φ ; (9) In summary, the magnetic field vector was easily deter- mined by imposing that the axial component of the vec- where A C = −B . The remaining parameters must be 1 1 0 tor potential must a differentiable function in any point determined by the boundary conditions at the surface of of space, in particular, at the superconductor surface. A the superconductor. detailed analysis of these solutions is presented in the The general solution of the London equation in the next section. interior region, for the axial component of the vector po- tential, is somewhat more complicated: P+1 DISCUSSION OF RESULTS Az = n=0 [En cos (nφ) + Fn sin (nφ)] × [GnIn(ρ/λ) + HnKn(ρ/λ)] ; (10) Several plots of magnetic field energy density are shown in Figure 1, for different values of the London pen- where In(ρ/λ) and Kn(ρ/λ) are the cylindrical modified Bessel functions of the first and second kind, respectively. etration depth. Each of the pictures represents a cross As the modified Bessel functions of the second kind di- section of the cylinder in the x-y plane ranging from mi- verge at ρ = 0, we can exclude them a priori. nus 3 to 3 in steps of 0.01 radius units for both x and y directions. The magnetic energy density was computed To ensure continuity at the superconducting surface, 2 from equations (16) in units of B0 , i.e. the magnetic en- the Az function must take the following form on the inner region, ergy density at an infinite distance. It is clearly visible that the purple circle, corresponding to field expulsion, 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.