Attractive vortex interaction and the intermediate-mixed state of superconductors Ernst Helmut Brandt1 and Mukunda P. Das2 1Max-Planck-Institut f¨ur Metallforschung, D-70506 Stuttgart, Germany 2Deptartment of Theoretical Physics, The Australian National University, Canberra, ACT 0200, Australia (Dated: July 8, 2010) The magnetic vortices in superconductors usually repel each other. Several cases are discussed when the vortex interaction has an attractive tail and thus a minimum, leading to vortex clusters and chains. Decoration pictures then typically look like in the intermediate state of type-I super- conductors, showing lamellae or islands of Meissner state or surrounded by Meissner state, but with the normal regions filled with Abrikosov vortices that are typical for type-II superconductors in the mixed state. Such intermediate-mixed state was observed and investigated in detail in pure Nb, TaN and other materials 40 years ago; last year it was possibly also observed in MgB2, where it was called “a totally new state” and ascribed to the existence of two superconducting electron bands, one of type-I and one of type-II. The complicated electronic structure of MgB2 and its consequences for su- perconductivity and vortices are discussed. It is shown that for the real superconductor MgB2 which possesses a single transition temperature, the assumption of two independent order parameters with separate penetration depths and separate coherence lengths is unphysical. PACS numbers: 2 1. VORTEX LATTICE FROM GL THEORY dense triangular lattice with induction B = (2/√3)Φ0/a (a = vortex spacing). When Ba and B reach the upper 2 critical field Bc2 = µ0Hc2 =Φ0/(2πξ ), the order param- The most successful phenomenological theory of su- eter vanishes and the superconductor becomes normal perconductors was conceived in 1950 by Ginzburg and conducting. Since B(Ba) increases monotonically (see Landau (GL) [1]. When written in reduced units (length Fig. 1 (top)), one may say that the pressure of the vortex unit λ, magnetic field unit √2Hc, energy density unit 2 lattice is positive and the vortices repel each other, held µ0H , where λ is the magnetic penetration depth and c together by the applied magnetic field Ba, which acts as Hc the thermodynamic critical field) the GL theory con- an external pressure. From thermodynamics follows that tains only one parameter, the GL parameter κ = λ/ξ (ξ in equilibrium (in absence of pinning and surface bar- is the superconducting coherence length). The kind of riers) one has B = µ0H = µ0∂F/∂B where F (B,κ) √ a a solutions of GL theory is quite different when κ< 1/ 2 is the GL free energy density. The numerical GL solu- or κ > 1/√2. Physically, this is due to the fact that tion Ba(B) of the vortex lattice [3] for 1/√2 κ 10 the energy of the wall between a normal and a super- ≤ ≤ may be fitted well by the expression (with b = B/Bc2, conducting domain is positive when κ < 0.71 (type-I) h = Ba/Bc2, m = b h = M/Bc2, hc1 = Bc1/Bc2): and negative when κ> 0.71 (type-II). This implies that − 3 a type-II superconductor is unstable against the sponta- c1b neous formation of many small superconducting and nor- h(b) = hc1 + 2 , 1+ c2b + c3b mal domains. A thesis student of Lev Landau in Moscow, ln κ + α(κ) 1 1 + ln 2 Alexei Abrikosov in 1953 discovered a periodic solution h 1 = , α(κ)= + , c 2κ2 2 √ of the GL theory which describes the occurrence of a reg- 2κ 2+2 3 2 − arXiv:1007.1107v1 [cond-mat.supr-con] 7 Jul 2010 c1 = (1 h 1) /(h 1 p) , ular lattice of vortex lines (published in 1957 [2]). Each − c c − Abrikosov vortex (or fluxon, flux line) carries one quan- c2 = (1 3hc1 +2p)/(hc1 p) , −15 2 − − tum of magnetic flux Φ0 = h/2e = 2.07 10 Tm , 2 · c3 = 1+(1 hc1)(1 2hc1 + p)/(hc1 p) , and the supercurrent circulates around a singular line on − − − − p = [(2κ2 1)β + 1 ] 1 (1) which the complex GL function ψ(x,y,z) (or supercon- − A ducting order parameter ψ 2) is zero. For finding this with β = 1.15960 (1.18034) for the triangular (square) vortex lattice solution, in| which| Landau initially did not A vortex lattice. At κ = 1/√2 the B(B ) is a vertical believe [2], Abrikosov obtained the Nobel prize in physics a straight line at B = B 1 = B 2 like in type-I supercon- 50 years later in 2003. a c c ductors. This means the vortices do not interact with Within the original GL theory in bulk superconduc- each other and all vortex arrangements have the same tors, with increasing applied magnetic field Ba = µ0Ha free energy. when κ =1/√2. the superconductor is first in the Meissner state which The vanishing vortex interaction for κ = 1/√2 can has induction B = 0. When Ba reaches the lower critical also be seen from Kramer’s [4] asymptotic interaction field B 1 Abrikosov vortices start to penetrate. With fur- between two vortices at large distance r λ, V (r) = c ≫ ther increasing B the vortex lines form a more and more d1(κ)K0(r/λ) d2(κ)K0(√2r/ξ) where the constants are a − 2 1 0.9 Ginzburg−Landau Theory 0.8 Type−II Superconductors 0.7 0.6 c2 0.5 B / 0.4 0.3 B = B a 0.2 0.76 κ = 0.1 23 1.5 1.2 κ = 1 0.85 0.707 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 B / B a c2 1 BCS−Gor’kov Theory 0.9 κ = 0.85, T/T = 0.7 c 0.8 FIG. 2: Intermediate state of a type-I superconductor. The B = B 0.7 a normal domains are dark. Coexistence of triangular lattice of flux tubes and of laminar domains. Tantalum disk of tickness 0.6 c2 d = 33 µm, diameter D = 5 mm, T = 1.2 K , Ba = 34 mT. 0.5 (Courtesy U. Essmann). Optical microscope. B / B GL−Theory 0.4 0 0.3 0.2 0.1 0 0.75 0.8 0.85 0.9 0.95 1 B / B a c2 FIG. 1: Top: Ideal reversible magnetization curves of a type- II superconductor with no vortex pinning and no demagne- tizing effects (long cylinder or slab in a parallel magnetic field Ba) from Ginzburg-Landau Theory. Plotted is the induction B versus Ba for several GL parameters κ = 1/√2 ... 3. Bot- tom: The same magnetization curve B(Ba) from the micro- scopic BCS-Gor’kov theory for GL parameter κ 0.85 and ≈ temperature T/Tc 0.7, schematic from [28]. The dashed ≈ part for B < B0 is unphysical, the real B(Ba) contains the depicted vertical line from B =0 to B = B0 obtained from a Maxwell construction that divides the area under the dashed curve into equal halves. d1 = d2 at κ = 1/√2. A similar vortex potential [5] was derived such that it reproduces the nonlocal elastic FIG. 3: Laminar domains in the intermediate state of pure properties of the triangular vortex lattice [6] in the entire superconductors at T = 1.2 K. a) Type-I superconductor, Tantalum disc, d = 33 µm, D = 4 mm, Ba = 68 mT, angle ranges of B and κ and for all wavelengths of the elastic ◦ strain. For parallel vortices this approximate interaction between magnetic field and disc 15 . Optical microscope. b) Type-II superconductor, Niobium disc, d = 40 µm, D = 4 potential reads mm, Ba = 74 mT. Optical microscope. c) As b) but electron 2 microscope. From [11]. Φ0 ′ ′ V (r)= ′2 [ K0(r/λ ) K0(r/ξ ) ] (2) 2πλ µ0 − with the lengths λ′ = λ/ ψ 2 1/2 λ/√1 b, ξ′ = equal to half the self energy of a vortex with flux 2Φ0, | | ≈ − since the two terms diverging as ln r compensate. ξ/p2(1 b) and K0(x) a modified Bessel function with ± the limits− ln x (x 1) and (π/2x)1/2 exp( x) (x 1). For κ > 1/√2 the potential (2) is repulsive. At κ = Physically,− the first≪ (repulsive) term in V (r−) Eq. (2)≫ de- 1/√2 it vanishes at all distances, which means that all scribes the magnetic interaction and the second (attrac- vortex arrangements have the same energy. This is an tive) term originates from the overlap of the vortex cores exact result of the GL theory [7]. The particular case 2 that reduces the condensation energy 1 ψ . Note κ =1/√2 and its surroundings (κ 0.71, T Tc) were that this V (r) has a finite maximum V (0),∝ approximately− | | investigated in detail in [8, 9]. ≈ ≈ 3 FIG. 5: Decoration of a square disk 5 5 1 mm3 of high- × × purity polycrystalline Nb at T = 1.2 K and Ba = 1100 Gauss. As Ba is increased, magnetic flux penetrates from the edges in form of fingers which are composed of vortex lattice shown enlarged in the right picture. The rectangular cross section of the disk causes an edge barrier (section 8.5.). As soon as this edge barrier is overcome, single flux lines or droplets of vortex lattice (lower right) pull apart from these fingers and jump to the center, filling the disk with flux from the center. (Courtesy U. Essmann) FIG. 4: Flux-line lattice observed at the surface of type-II superconductors in an electron microscope after decoration with Fe microcrystallites (“magnetic smoke”).