Arxiv:1007.1107V1 [Cond-Mat.Supr-Con] 7 Jul 2010 Adu(L 1.We Rte Nrdcduis(Length Units Reduced and in Ginzburg Written Unit When by [1]

Home , Abrikosov vortex

arXiv:1007.1107v1 [cond-mat.supr-con] 7 Jul 2010 adu(L 1.We rte nrdcduis(length units reduced and in Ginzburg written unit When by [1]. 1950 (GL) in Landau conceived was perconductors h uecnutri rti h esnrsaewhich state Meissner the in induction first has is superconductor the field hrincreasing ther u fmgei u Φ quan- flux one Each magnetic carries line) [2]). of flux 1957 tum fluxon, in (or (published reg- vortex lines a Abrikosov vortex of occurrence of the lattice solution describes ular which periodic theory a GL discovered the of Moscow, 1953 in in Landau Lev Abrikosov of Alexei student thesis nor- A and domains. superconducting mal sponta- small the many of against formation unstable neous is superconductor type-II a os ihicesn ple antcfield magnetic applied increasing with tors, 2003. in physics in later prize years not Nobel 50 the did obtained initially Abrikosov Landau [2], which believe in solution, lattice vortex utn re parameter order ducting hc h ope Lfunction on GL line complex singular the a which around circulates supercurrent the and ouin fG hoyi ut ieetwhen of different kind quite is The theory GL length). of coherence solutions superconducting the is only tains H h nryo h albtenanra n super- a and normal when when negative a positive and between is wall domain the conducting of energy the or µ 0 c h otscesu hnmnlgclter fsu- of theory phenomenological successful most The ihnteoiia Lter nbl superconduc- bulk in theory GL original the Within H trcievre neato n h intermediate-mixed the and interaction vortex Attractive > κ h hroyai rtclfil)teG hoycon- theory GL the field) critical thermodynamic the .VRE ATC RMG THEORY GL FROM LATTICE VORTEX 1. c 2 B λ where , antcfil unit field magnetic , c 1 1 2 bioo otcssatt eert.Wt fur- With penetrate. to start vortices Abrikosov / etrmn fTertclPyis h utainNation Australian The Physics, Theoretical of Deptartment n parameter one √ ASnumbers: PACS length coherence separate and depths assumption penetration the temperature, separate transition single a possesses ecnutvt n otcsaedsusd ti hw tha shown is It discussed. are vortices and perconductivity attlynwsae n srbdt h xsec ftosup two struc of electronic existence complicated the The to type-II. of ascribed one and and type-I state” new totally “a n te aeil 0yasao atya twspsil al possibly was it year last ago; years 40 an observed materials ty was other are state and that intermediate-mixed Such vortices Abrikosov state. with o mixed filled in state Meissner regions like of normal a look islands the thus typically or and lamellae then showing tail pictures conductors, attractive Decoration an chains. has and interaction vortex the when .Pyial,ti sdet h atthat fact the to due is this Physically, 2. B λ h antcvrie nsprodcosuulyrpleac repel usually superconductors in vortices magnetic The B stemgei eerto et and depth penetration magnetic the is .When 0. = a h otxlnsfr oeadmore and more a form lines vortex the > κ 0 0 1 h Lparameter GL the , . a-lnkIsiu ¨rMtlfrcug -00 Stutt D-70506 für Metallforschung, Max-Planck-Institut | 1(yeI) hsipisthat implies This (type-II). 71 ψ = | B 2 √ h/ szr.Frfidn this finding For zero. is ) a 2 H ece h oe critical lower the reaches 2 ψ e c rs emtBrandt Helmut Ernst ( nrydniyunit density energy , 2 = ,y z y, x, < κ . 07 o supercon- (or ) 0 10 . B 1(type-I) 71 κ a < κ − = 15 = Dtd uy8 2010) 8, July (Dated: λ/ξ µ 1 Tm / 0 √ H ( 2 2 a ξ , 1 n uud .Das P. Mukunda and with otxltie At lattice. vortex es raglrltiewt induction with lattice triangular dense ( h lob enfo rmrs[]aypoi interaction distance asymptotic large [4] at Kramer’s vortices two from between seen be also a efitdwl yteepeso (with expression the by well fitted be may rtclfield critical tagtln at line straight ahohradalvre ragmnshv h same the with have interact arrangements when not vortex energy. do all free vortices and the other means each This ductors. steG reeeg est.TenmrclG solu- GL numerical The density. energy free tion GL the is trvnse n h uecnutrbcmsnormal becomes Since superconductor the conducting. and vanishes eter atc spstv n h otcsrplec te,held field other, magnetic each applied repel the vortex vortices by the the of together and pressure the positive that is say lattice may one (top)), 1 Fig. neulbim(nasneo inn n ufc bar- surface and pinning has of one absence riers) (in that equilibrium follows thermodynamics in From pressure. external an d a 1 = ( h aihn otxitrcinfor interaction vortex vanishing The otxsaig.When spacing). vortex = κ B ) h lUiest,Cner,AT00,Australia 0200, ACT Canberra, University, al B ftoidpnetodrprmtr with parameters order independent two of β K h o h elsprodco MgB superconductor real the for t ueo MgB of ture ( a uruddb esnrsae u with but state, Meissner by surrounded r A a sunphysical. is s oosre nMgB in observed so c c c b c /B 0 netgtdi eali ueN,TaN Nb, pure in detail in investigated d ia o yeI uecnutr nthe in superconductors type-II for pical p 3 2 1 1 ( h nemdaesaeo yeIsuper- type-I of state intermediate the = ) ( B 1 = r/λ te.Svrlcssaediscussed are cases Several other. h iiu,laigt otxclusters vortex to leading minimum, (2 [ = (1 + 1 = (1 = (1 = = c ftevre atc 3 o 1 for [3] lattice vortex the of ) rodcigeeto ad,oeof one bands, electron erconducting 2 , B . ) 56 1104 o h raglr(square) triangular the for (1.18034) 15960 h m ln − c c 2 B 1 B − − κ d = κ = a 2 + 2 2 a 2 + at Germany gart, 3 h ( κ n t osqecsfrsu- for consequences its and κ b − µ = κ h = c + 1 − 2 B − α κ 1 0 ) 1 = c tt fsuperconductors of state 2 1) K ) H 1 ( µ ( B h 3 h κ B 1 = 2 + β 0 / 0 c c c c ) / 2 H 1 = ( ( c 1 a A 2 2 α , √ √ )(1 h 1 b hr twscalled was it where , nrae oooial (see monotonically increases ) Φ = a = p b ] 1 + c M/B / + 2. 2 3 1 ) √ r/ξ = − / B ( − c ( B κ 0 the 2 3 c h 2 / − = ) µ 2 p b a hr h osat are constants the where ) h c c (2 2 1 ) 0 1 2 iei yeIsupercon- type-I in like c 2 and ∂F/∂B , 1 , πξ − , + h 1 2 p c B 2 2 + 1 p ,teodrparam- order the ), ) B B ( which ) , B = B / 2 r ec h upper the reach ( a (2 = κ a κ B h n2 ln + 1 / hc csas acts which , where savertical a is ) − ≫ c √ c 1 = 1 1 √ b 2 − /B / λ 2+2 √ = ≤ p , / c ) 3)Φ 2 F √ V 2 B/B κ ): , ( can 2 ( , ,κ B, ≤ r 0 = ) /a (1) c 10 2 2 ) , 2

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 ﬂux 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 eﬀects (long cylinder or slab in a parallel magnetic ﬁeld 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 microscopic 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”). Top: High- purity Nb disks 1 mm thick, 4 mm diameter, of different crys- mm, at T =1.2 K and Ba = 34 mT, taken in 1969, from tallographic orientations [110] and [011], at T = 1.2 K and the review [11]. The picture shows the coexistence of Ba = 800 Gauss (Bc1 = 1400 Gauss). Due to demagnetiz- triangular lattice of flux tubes and of laminar domains. ing effects and the small κ 0.70, round islands of Meissner ≈ As recently shown by Prozorov [12, 13], the equilibrium phase are surrounded by a regular vortex lattice (“interme- topology of the intermediate state are flux tubes (in ab- diate mixed state”). Bottom: High-purity Nb foil 0.16 mm sence of bulk pinning and geometric barrier). Interest- thick at T = 1.2 K and Ba = 173 Gauss. Round islands of vortex lattice embedded in a Meissner phase. (Courtesy U. ingly, pin-free disks and strips of constant thickness show Essmann) hysteretic magnetization curves due to edge barriers for penetrating flux tubes and exiting lamellae. This topo- logical hysteresis of type-I superconductors was observed and calculated analytically already in 1974 [14], where it 2. TYPE-I SUPERCONDUCTORS was shown that ellipsoidal samples exhibit no magnetic hysteresis. Magnetic hysteresis loops in pin-free type-I Superconductors with κ < 1/√2 are called as type superconductor disk and strips look very similar as the I. For Ba < Bc = µ0Hc (thermodynamic critical field) irreversible magnetization loops in pin-free type-II super- these should exhibit ideal Meissner effect with B = 0 conductors with edge barrier [15], in particular, they go (except for a surface layer of thickness λ) and B = B through the origin (M = B B = 0 at B = 0, no a − a a for Ba >Bc. But in samples of finite size the demagneti- flux is trapped) and become reversible at large Ba (at zation factor N > 0 allows the observation of an “inter- Ba >Bc1/2 for type-II). mediate state” containing superconducting and normal The structure of the flux tubes or lamellae in type-I domains with induction B = 0 and B = Bc, respectively, superconductors and their widening (mushrooming) and in the range (1 N)Bc

[22]. In the type-II superconductors Vanadium and Nio- bium (see Fig. 3 of [11]) and PbIn alloys they showed more or less regular triangular vortex lattice with struc- tural defects, e.g., dislocations and vacancies. These pictures were consistent with the prediction from GL theory that the vortices repel each other. However, later decoration of very pure Nb showed a domain structure similar to that observed in type-I superconductors but with the normal domains replaced by domains of vortex lattice. The Figs. 3a and 3b compare the laminar domains in the type-I superconductor Tantalum and in the type- II superconductor Nb. In an optical microscope these domains look very similar, but the electron microscope (Fig. 3c) reveals that in Nb the domains which carry the magnetic flux consist of vortex lattice with constant lattice spacing a0. In Fig. 4 (top) one sees islands of Meissner state embedded in a vortex lattice, while Fig. 4 (bottom) shows Meissner state with islands of vortex lattice. In Fig. 5 FIG. 6: Ideal magnetization curves of long cylinders or slabs tongues of vortex lattice penetrate from the specimen in a parallel magnetic field Ba (demagnetization factor N = 0, edge into the Meissner state of high-purity polycrystalline left) and of a sphere (N = 1/3, right). Top row: Type-I Nb. superconductor with Meissner state for Ba < (1 N)Bc, in- − This observation of an upper limit a0 of the vortex termediate state for (1 N)Bc Ba Bc, and normal state − ≤ ≤ spacing in the intermediate state of low-κ type-II super- for Ba > Bc. The dashed line indicates the enhancement of the penetration field when the wall energy of the normal conductors, and its explanation by an attractive inter- and superconducting domains in the intermediate state is ac- action, suggests that the magnetization curves in the counted for. Middle row: Type-II/1 superconductors with ideal (pin-free) case should exhibit a jump of height 2 Ginzburg-Landau parameter κ close to 1/√2, for N = 0 ex- B0 = (2/√3)Φ0/a0 at the field of first vortex penetra- hibit a jump of height B0 in their induction and magnetization. This means, as soon as it is energetically favorable tion at Ba = Bc1 due to an attractive interaction between flux for vortices to penetrate, they immediately jump to their lines. For N > 0, this jump is stretched over a finite range of equilibrium distance a0, which for N = 0 means a uni- Ba, allowing one to observe an intermediate mixed state with domains of Meissner phase and vortex-lattice phase. Bottom form induction B0, see Fig. 1 (bottom). Such jumps of row: Type-II/2 superconductors with κ 1/√2. Meissner B(Ba) from 0 to B0 were indeed observed in clean Nb ≫ phase for Ba < (1 N)Bc1, mixed state (vortex lattice) for and in lead alloys with increasing content of Thallium − (1 N)Bc1 Ba Bc2, and normal state for Ba > Bc2. (which changes κ), and excellent agreement of B0 with − ≤ ≤ the directly measured a0 was observed [23–25]. Jumps B0, and thus vortex attraction, were also observed in the for κ = 0.5 the shear modulus c66 of this vortex lattice TaN system with varying nitrogen content and various is positive at all inductions 0 superconductivity (type-I or type-II). We pends heavily on the geometry of the Fermi surfaces. also discussed the behavior of pure Nb at temperatures Then comes the electron-phonon interaction that medi- below the GL regime. In the conventional single gap ates the pairing. In view of the complicated band struc- superconductor within GL theory there are a well defined tures and their Fermi surfaces, the novel superconductors coherence length ξ and penetration depth λ. In terms of are treated beyond the standard BCS model. We are them one defines κ = λ/ξ. As stated earlier its unique reminded that single-sheeted Fermi surfaces only come value, if less or bigger than 0.71, defines precisely type- from the simple metals like alkalies considered in the jel- I or type-II superconductors, in which the surface free lium model and those are not superconductors under am- energy between normal and Meissner domains is positive bient conditions. Therefore, if the multi-sheeted Fermi or negative, respectively. This definition tells about the surfaces are connected in the entire Brilliouin zone, then equilibrium states of a superconductor. Here with this an anisotropic GL theory is useful. Otherwise for uncon- brief background we shall discuss the interesting physics nected Fermi surfaces one can use either (i) a multi-band of the novel superconductor MgB2 and the vortex matter BCS theory [52] or (ii) a generalization of the GL theory therein. Needless to say that it has become a news item invoking many (or many components of) order param- in the past one year with the claim that a new form of eters [53]. This two-component order-parameter theory matter is discovered [49, 50]. Detailed discussions follow was used to explain successfully two specific-heat transi- in Sec. 6. tions in the high Tc superconductor YBCO (see Choy et al. in [53]). As said before, in reality the bands are formed due 5. VORTICES IN MGB2 to ionic and electronic self-consistent potential. In principle we are dealing with a charged (Coulomb) system, In 2001 by a surprising discovery Nagamatsu et al. [57] hence the bands are not entirely independent, and further reported Tc of an intermetallic MgB2 around 39 K. Soon due to the electron-phonon interaction (whether small after (see [58] database for the first three years) it was or large) the bands are coupled. For this reason multi- realized that MgB2 offers a host of novel properties that (inter- and intra-)band interaction terms explicitly ap- has been explored during the past 9 years. Magnesium pear in the microscopic Hamiltonian [52]. In the case diboride has a crystallographic structure consisting of of multi-order parameter GL theory, one has the super- honeycomb boron layers separated by magnesium layers, conducting free energy composed of individual parts of which have also honeycomb structure. Magnesium atoms the GL free energy and the additional free energy due are ionized in this structure donating their s-electrons to to coupling between the order parameters as mentioned the conduction band. There are two 3-dimensional π above. This choice of coupling is non-unique. If ψ1 and bands arising from the boron pz orbitals. One of them ψ2 are two order parameters the coupling or interaction is electron-like, whereas the other is hole-like. From the term can be one or more terms from the following list: px,y orbitals two σ bands occur. They are 2-dimensional 2 2 ∗ + + + confined to the boron planes. Bonds within the boron γ1 ψ1 ψ2 , γ2(Π ψ Πψ2 + hc), γ3(ψ ψ2 + ψ1ψ ) | | | | 1 1 2 layers are strongly covalent, whereas between the layers with Π = +2πiA/Φ0. The first term of this list is they are metallic. From this a complicated Fermi surfaces a systematic∇ expansion of GL, the second term means geometry emerges. Two cylinders around the Γ-A-Γ lines coupling through the gradient of order parameter and are two σ bands, and two webbed tunnels are attributed vector potential A, and the third term is the internal to π bands, see Fig. 7 (from Mazin and Antropov [59]). Josephson coupling term. This last term is supposed to This overall picture is consistent with all the electronic provide a “minimal coupling” [54–56]. Instead of the phe- structure calculations [59]. nomenological GL theory, we start with the microscopic It is further found [60] that the optical bond-stretching Bogolyubov pairing Hamiltonian and use the mean-field E2g phonons couple strongly to the holes at the top approach; we derive the GL functional with a coupling of σ bands, whereas the 3-dimensional π electrons term that will be Josephson-like. An important point are weakly coupled to the phonons. For this reason here is that the coefficients α of ψ 2, and the coefficient of different coupling strengths of σ and π bands, a i | i| γ3 of the Josephson term, do depend both on the intra- two-gap superconductor merges. The E2g phonons band and inter-band pairing. Therefore, as a matter of strongly coupled with σ band, produce high Tc in principle the order parameters can never be independent. MgB2. A further coupling with the π bands enhances 7

In the past nine years a variety of experiments have shown that MgB2 is a two band superconductor. The main experiments are speciﬁc heat (Choi et al. in [59, 62], point contact and tunnelling spectroscopy [63, 64], Ra- man spectroscopy [65], penetration depth measurements [66], angle-resolved photo-emission spectroscopy [67], small-angle neutron-scattering [68], and muon-spin re- laxation studies [69] etc. The value of the gap on the σ bands ∆σ ranges from 6.4 to 7.2 meV, while on the π bands ∆π varies from 1.2 to 3.7 meV. Almost all parameters such as coherence lengths (ξσ and ξπ) [66], upper critical ﬁelds (Hc2(σ) and Hc2(π)) [70, 71], and penetration depths (λσ and λπ) [66] are anisotropic. Besides this the anisotropic factor is also strongly temperature dependent. These experiments convincingly suggest that MgB2 is a two band superconductor, where phonons mediate pairing in the respective bands. Further, it is investigated that inter-band coupling helps to increase Tc. In the absence of inter-band coupling, it may give two isolated superconductors with two lower Tc.

6. A NEW PHASE OF MGB2 !

In a recent paper Babaev and Speight [72] studied phenomenologically a two-band superconductor with two order parameters, starting with a free energy that in- cludes a Josephson-like interaction γ3 between the two FIG. 7: Band structure and Fermi surface of MgB2 from bands. Choosing γ3 0 in the decoupled band limit, → Mazin and Antropov [59]. Green and blue cylinders (hole- they considered κi (i = 1, 2) in two different regimes like) are from the σ bands, and the blue (hole-like) and red to produce type-I (κ1 < 0.71) and type-II (κ1 < 0.71) (electron-like) tubular networks are from the π bands. materials. Their prediction leads to, what they call a “semi-Meissner state”, particularly in the case of very low magnetic fields. Babaev and Speight claim that this T . Table 1 presents values of intra- and inter-band c state is thermodynamically stable. Instead of homoge- (electron-phonon) coupling-constant matrix-elements neous distribution, the vortices form aperiodic clusters or calculated by various groups. The numbers show that vortexless Meissner domains, arising out of short range the electron-phonon coupling λ (a dimensionless param- repulsion and long range attraction. Below we shall argue eter [61] not to be confused with the penetration depth that in the particular case of MgB2, this is an unrealistic λ) for the σ band is very strong and for the π band construction. it is weak. The inter-band coupling is still weaker but In what is claimed as a remarkable manifestation of appreciable. These numbers are useful ingredients to be a new type of superconductivity in MgB2, Moshchalkov incorporated in the two-order parameter calculations. and co-workers [49, 50] reported last year a novel “type- In a first principles calculation using Eliashberg theory, 1.5 superconductor”, contrary to type-I and type-II su- Geerk et al. [60] emphasize the importance of inter-band perconductors. It appears these authors are influenced by pairing and low energy optical modes arising out of the idea of “semi-Meissner state” of Babaev and Speight. boron atoms vibrating perpendicular to the boron planes. Apparently both names carry similar meaning. Let us analyse carefully the results of Moshchalkov et al. λ λ λ λ σσ ππ σπ πσ Similar to Babaev and Speight, Moshchalkov et al. use 1.017 0.448 0.213 0.155 Golubov et al. [60] a two-band GL functional with an inter-band coupling term. They constructed the GL parameters κ = λ /ξ 0.96 0.29 0.23 0.17 Liu et al. [60] i i i (i stands for σ and π bands of MgB2) from the mea- 0.78 0.21 0.15 0.11 Choi et al. [59] sured values of energy gaps, Fermi velocities and plasma frequencies of respective bands. The values are so cho- Table 1. Matrix elements of electron-phonon interactions sen that κσ > 0.71 (=3.68) and κπ < 0.71 (= 0.66). 8

Thus, MgB2 has properties of both type-I and type-II mixed state to occur is that the vortex interaction has simultaneously, unlike in the case of Nb where κ is close an attractive tail, or that the theoretical magnetization to 0.71 as discussed above. This construction is in ab- curve is S-shaped, see Fig. 1 (bottom). But in addition, sence of the coupling term, so the bands are fully in- the demagnetization factor of the decorated specimen has dependent/noninteracting. We have stated earlier that to be N > 0, else only Meissner state or a uniform vortex energy gaps ∆σ range from 6.4 to 7.2 meV, while the ∆π state will occur; the detailed appearance further depends varies from 1.2 to 3.7 meV. Moshchalkov et al. chose the on pinning, on the edge barrier, and on the magnetic value of ∆σ = 7.1 meV and ∆π = 2.2 meV from some history. This intermediate mixed state has been known experiments and computed ξi from the one-band BCS for long time from experiments [11, 22–26, 34, 75] and 2 relation ξi = hvF /(2π ∆i). Then they used calculated computations [28, 29, 31–33] for pure Nb, TaN, PbIn, plasma frequencies of respective bands, from which they and PbTl. In this respect, the decoration and scanning obtained the London penetration depths λσ(0) = 48.8 SQUID results of [49, 50] are not so new. nm and λπ(0) = 33.6 nm and finally κi = λi/ξi. As discussed above, several origins of vortex attrac- Here we note a few points. tion are conceivable, depending on the material, its pu- A. (1) The numerical estimates for ξi are obtained by rity, magnetic history, and temperature. Computer sim- using the one-band BCS formula. (2) One band calcu- ulations using such a vortex interaction (short-range re- lated plasma frequency is used to compute λi. One does pulsive, long-range attractive) yield vortex arrangements not know how these are relevant to interacting two-band that are similar to the observed ones [49, 50]. The details MgB2. of the resulting clusters or chains also depend on the type B. For critical fields Moshchalkov et al. gave their mea- and strength of vortex pinning that is always present and sured values for Hc1(0) = 241 mT, Hc2(0) = 5.1 T for the should be included in such simulations, in particular at σ band. From these data they also obtained λσ and ξσ very low induction B. from the one band formula. The value of Hc1(0) = 241 mT measured for their sample is much larger than the values reported by others, which are in the range of 25 7. SUMMARY AND CONCLUSIONS - 48 mT, see Fig. 25 of Ref. [71]. Another point of con- cern is the reported value of thermodynamic critical field Type-I superconductors with demagnetization factor Hc(0) = 230 mT [73], whereas Moshchalkov et al. have N > 0 in the field range (1 N)B

[32] U. Klein, J. Low Temp. Phys. 69, 1 (1987). [56] Internal Josephson coupling can be obtained from the [33] H. W. Weber, M. Botlo, F. M. Sauerzopf, H. Wiesinger, two-band interacting microscopic model, where it arises and U. Klein, Jap. J. Appl. Phys. 26, 917 (1987). from inter-band interaction. Its presence leads to the pos- [34] M. Laver et al., Phys. Rev. B 79, 014518 (2009); sibility of an internal collective mode, known as Leggett’s see also J. Shelten, H. Ullmaier, and W. Schmatz, mode [54]. Attempts are made to observe this mode in 48 phys. stat. sol.(b) , 619 (1971). the Raman spectroscopy for MgB2 without clear success [35] E. H. Brandt, Phys. Lett. 39A, 227 (1972). (see S. G. Sharapov et al. and G. Blumberg et al. in [36] A. B. Pippard, Proc. Roy. Soc. (London) A216, 547 Ref. [54]). (1953). [57] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, [37] E. H. Brandt, Phys. Lett. 39A, 193 (1972). and J. Akimitsu, Nature 410, 63 (2001). [38] M. C. Leung, J. Low Temp. Phys. 12, 215 (1973). [58] http://www.iitap.iastate.edu/htcu/mgb2preprints.html. [39] A. M. Grishin, A. Yu. Martinovich, and S. V. Yam- [59] J. Kortus et al., Phys. Rev. Lett. 86, 4656 (2001); J. M. pol’skiy, Zh. Exp. Teor. Fiz. 97, 1930 (1990) [Sov. Phys. An and W. E. Pickett, Phys. Rev. Lett. 86, 4366 (2001); JETP 70, 1089 (1990). H. Harima, Physica C 381, 18 (2002); H. J. Choi et al., [40] A. I. Buzdin and A. Yu. Simonov, Physica C 168, 421 Nature 418, 758 (2002); I. I. Mazin and V. P. Antropov, (1990). Physica C 385, 49 (2003). [41] V. G. Kogan, N. Nakagawa, S. Thiemann, Phys. Rev. B [60] K. P. Bohnen et al., Phys. Rev. Lett. 86, 5771 (2001); A. 42, 2631 (1990). Liu et al., Phys. Rev. Lett. 87, 087005 (2001); T. Yildirim [42] M. V. Feigelman, V. B. Geshkenbein, and A. I. Larkin, et al., Phys. Rev. Lett. 87, 037001 (2001); Y. Kong et al. Physica C 167, 177 (1990). Phys. Rev. B 64, 020501 (2001); A. A. Golubov et al., J. [43] J. R. Clem, Phys. Rev. B 43, 7837 (1991). Phys.: Condens. Mat. 14, 1353 (2002); J. Geerk et al., [44] D. J. C. Jackson and M. P. Das, Solid State Comm. 90, Phys. Rev. Lett. 94, 227005 (2005). 479 (1994). [61] E. H. Brandt and N. B. Brandt, Zh. Eksp. Teor. Fiz. 60, [45] D. J. C. Jackson and M. P. Das, Supercond. Sci. Technol. 1141 (1971) [Sov. Phys. JETP 33, 618 (1972)]. 9, 713 (1996). [62] F. Bouquet et al., Phys. Rev. Lett. 87, 047001 (2001); H. [46] E. H. Brandt, R. G. Mints, and I. B. Snapiro, Phys. Rev. D. Yang et al., Phys. Rev. Lett. 87, 167003 (2001). Lett. 76, 827 (1996). [63] P. Szabo et al., Phys. Rev. Lett. 87, 137005 (2001); R. [47] G. Blatter and V. B. Geshkenbein, Phys. Rev. Lett. 77, S. Gonnelli et al., Phys. Rev. Lett. 89, 247004 (2002); H. 4958 (1996). Schmidt et al., Phys. Rev. Lett. 88, 127002 (2002). [48] S. Mukherji and T. Nattermann, Phys. Rev. Lett. 79, [64] M. Iavarone et al., Phys. Rev. Lett. 89, 187002 (2002); M. 139 (1997). R. Eskildsen et al., Phys. Rev. Lett. 89, 187003 (2002); [49] V. Moshchalkov, M. Menghini, T. Nishio, Q. H. Chen, F. Guibileo et al., Phys. Rev. Lett. 87, 177008 (2001). A. V. Silhanek, V. H. Dao, L. F. Chibotaru, N. D. Zhi- [65] G. Blumberg et al., Physica C 456, 75 (2005); X. K. gadlo, and J. Karpinski, Phys. Rev. Lett. 102, 117001 Chen et al., Phys. Rev. Lett. 87, 157002 (2001). (2009). [66] J. D. Fletcher et al., Phys. Rev. Lett. 95, 097005 (2005). [50] T. Nishio, V. H. Dao, Q. Chen, L. F. Chibotaru, K. Kad- [67] H. Uchiyama et al., Phys. Rev. Lett. 88, 157002 (2002); owaki, and V. Moshchalkov, Phys. Rev. B 81, 020506(R) S. Tsuda et al., Phys. Rev. Lett. 91, 127001 (2003); S. (2010); T. Nishio et al., Physica C, in print. Souma et al., Nature 423, 65 (2003). [51] E. H. Brandt and S.-P. Zhou, Physics 2, 22 (2009). [68] C. D. Dewhurst et al., Physica C 404, 135 (2004). [52] H. Suhl, B. T. Matthias, and L. R. Walker, Phys. Rev. [69] Ch. Niedermayer, Phys. Rev. B 65, 094512 (2002). Lett. 3, 552 (1959); V. A. Moskalenko, Phys. Met. Met- [70] M. Eisterer, Supercond. Sci. Technol. 20, R47 (2007); allogr. 8, 503 (1959); J. Kondo, Prog. Theor. Phys. 29, [71] C. Buzea and T. Yamashita, Supercond. Sci. Technol. 14, 1 (1963). R115 (2001). [53] M. P. Das, H.-X. He, and T. C. Choy, Inter. J. Modern [72] E. Babaev and M. Speight, Phys. Rev. B 72, 180502 Physics B 2, 1513 (1988); T. C. Choy, M. P. Das and H. (2005). He, Phase Transition A20, 1 (1990); For various gener- [73] B. Kang, H.-J. Kim, M. Park and S.-I. Lee, Chinese J. alisations of GL theory see M. P. Das in Studies in High Phys. 43, 583 (2005). Temperature Superconductors, Edited by A. V. Narlikar, [74] A. A. Golubov et al., Phys. Rev. B 66, 054524 (2002). Nova Science Publ. New York (USA) Vol.3, 267 (1989). [75] U. Essmann, private communication. [54] A. J. Leggett, Prog. Theor. Phys. 36, 901 (1966); S. G. [76] P. E. Goa, H. Hauglin. M. Baziljevich, E. Il’yashenko, P. Sharapov et al., Euro. Phys. J. B 30, 45 (2002); G. Blum- L. Gammel and T. H. Johansen, Supercond. Sci. Technol. berg et al., Phys. Rev. Lett. 99, 227002 (2007). 14, 729 (2001). [55] See for example, T. D. Lee, Phys. Rev. D 8, 1226 (1973). [77] A. Gurevich, Phys. Rev. B 67, 184515 (2003). His Lagrangian Eq. (7) contains terms required to satisfy [78] X. X. Xi, Rep. Prog. Phys. 71, 116501 (2008). various invariances. See also M. E. Zhitomirsky and V. [79] I. N. Askerzade, Physica C 390, 281 (2003). H. Dao, Phys. Rev. B 69, 054508 (2004); A. Gurevich, [80] M. Ichioka, K. Machida, N. Nakai, and P. Miranovi´c, Physica C 456, 160 (2007); V. G. Kogan, Phys. Rev. B Phys. Rev. B 70, 144508 (2004). 80, 014507 (2009); I. N. Askerzade, Physics-Uspekhi 49, [81] A. Gumann, S. Graser, T. Dahm, and N. Schopohl, 1003 (2006); H. Doh et al., Phys. Rev. Lett. 83, 5350 Phys. Rev. B 73, 104506 (2006). (1999).