Giant microwave absorption in fine powders of superconductors.

Enhanced microwave absorption, larger than that in the normal state, is observed in fine grains of type-II superconductors (MgB2 and K3C60) for magnetic fields as small as a few % of the upper critical field. The effect is predicted by the theory of vortex motion in type-II superconductors, however its direct observation has been elusive due to skin-depth limitations; conventional microwave absorption studies employ larger samples where the microwave magnetic field exclusion significantly lowers the absorption. We show that the enhancement is observable in grains smaller than the penetration depth. A quantitative analysis on K3C60 in the framework of the Coffey-Clem (CC) theory explains well the temperature dependence of the microwave absorption and also allows to determine the vortex pinning force constant.

Introduction

Electrodynamics of superconductors remains an intensively studied field[1] due to the wealth of attainable fundamental information, including the nature of pairing mechanism and the coupling strength, and also due to the technological importance of these materials. In addition, recent studies focus on novel topological phases including superconductors[2]. As an example, observation of the conductivity coherence peak in conventional superconductors[3] (Nb and Pb) and its absence in high-Tc materials[4] pointed to a BCS mechanism in the former and it was an early indication of non-BCS superconductivity in the latter compounds. Concerning applications, the DC electrodynamic properties in the mixed state of type-II superconductors determine the utility (e.g. loss, permanent field homogeneity, and stability) in superconducting solenoid coils that are widely used in superconducting particle acceleration, solid state spectroscopy, or medical imaging. The AC electrodynamic properties are relevant for applications including e.g. power handling, sound and electromagnetic field detection[5-7], superconducting microwave resonators[8], and in microwave absorbers[9].

The frequency dependent conductivity of superconductors, , is well known for both BCS (i.e. weak-coupled s-wave pairing) and non-BCS superconductors (including strongly coupled s-wave and non s-wave superconductors) in the absence of magnetic field, B = 0. At zero temperature, T = 0, the real part, σ1(ω), is a delta function at ω = 0 followed by σ1(ω) = 0 until the gap edge at ωg = 2Δ/ℏ (ref.[10]) (usually at 0.1–10 THz). According to the Ferrell–Glover–Tinkham (FGT) sum rule[11,12], the spectral weight of the delta function comes from states which are gapped below ωg (the sum rule is discussed in depth in the Supplementary Material). The technologically important radio frequency range spans 9 orders of magnitude in superconductors (from 10 kHz up to 1 THz) with similar characteristic properties, thus measurements in the microwave range (1–100 GHz) are representative.

Conductivity in finite magnetic fields for the mixed state in type-II superconductors was first described by the Bardeen–Stephen (BS) model[13] for the viscous motion of vortices. This was later improved by the Coffey–Clem theory (CC) in a series of seminal papers[14-19], which also includes the effect of pinning force on the vortex motion. The most important prediction of the BS model is a finite σ1 conductivity at ω = 0. However, it is less well-known that the FGT sum rule implies a non-zero σ1 that is larger than in the normal state for ω < ωg. Observation of this effect has been elusive as most contributions study the surface impedance on polycrystalline[3], compacted powder pellet, or thin film samples[20,21]. Surface impedance studies have the advantage that sample geometry is well defined, however the effects of σ1 and σ2 are inevitably intermixed in this type of measurements. Given that σ2 is orders of magnitude larger than σ1 in the superconducting state (due to the small value of the penetration depth, λ, with respect to sample thickness), the surface impedance measurement is less sensitive to changes in σ1 (refs[22-26]).

The effects of σ1 and σ2 are decoupled for fine grains; for a sample placed in a microwave cavity, the loss is due to σ1, whereas the resonance shift is due to σ2 (refs[27-29]). Therefore such samples provide a unique opportunity to test the predictions of the CC theory on σ1,2(T, B). The enigmatic and yet unexplained increase of the electron spin resonance signal in superconductors right below Tc (refs[30,31]) also highlights the need to study further the electromagnetic absorption in superconductors.

This motivated us to revisit the microwave conductivity (at about 10 GHz) in the MgB2 and K3C60 superconductors as a function of T and B. We observe an excess microwave loss (or microwave absorption) in small magnetic fields (as low as a few % of the upper critical field, Bc2) for a sample consisting of well-separated fine grains (typical size is a few micrometers). The excess microwave absorption is not observable in a single crystal sample. A quantitative analysis is provided for K3C60, which is a one-gap, cubic superconductor with well known magneto-transport properties[32], whereas MgB2 is a multi-band superconductor with strongly anisotropic Bc2 (refs[33,34]), thus application of the CC model is less straightforward.

Methods and Experimental

We studied fine powder MgB2 samples identical to batches in refs[35] and[30]. Single crystal and powder K3C60 samples were prepared by the conventional K intercalation method; the crystal sample was from the same batch as in ref.[36]. The starting fullerene powder material contains grains with a size of a few micrometers that is retained in the K doped material according to literature studies[37-39]. The powder samples were further ground together with non-conducting SnO2 powder to prevent conducting links between the grains. Samples were sealed in quartz ampules under low pressure helium. Microwave properties were measured with the cavity perturbation method[27,28] as a function of temperature, T, and in various static magnetic fields, B, inside a superconducting solenoid, with zero-field cooling. Zero field measurements (besides the Earth’s magnetic field) were made in another cryostat without a magnet field solenoid to avoid trapped flux (which may amount to 10–20 mT). The unloaded copper cavity has a quality factor, Q0 ~ 10,000 and a resonance frequency, f0 ~ 11.2 GHz, whose temperature dependence is taken into account. The samples were placed in the node of the microwave (or rf) electric field and maximum of the microwave magnetic field inside the TE011 cavity, which is the appropriate geometry to study minute changes in the conductivity[29]. The rf magnetic field is parallel to the DC field of the solenoid, which yields the largest vortex motion induced absorption according to the CC theory[16]. Measurement[40] of the quality factor, Q, and the cavity resonance frequency, f yields the loss: and cavity shift: Δf/f0 = (f − f0)/f0.

Results

Figure 1 shows the microwave cavity loss and cavity shift for a fine powder of MgB2 and for two kinds of K3C60 samples: a single crystal and a fine powder as a function of temperature and for a few magnetic field values. The microwave loss decreases rapidly below Tc in zero magnetic field as expected for superconductors. The most important observation is that the microwave loss becomes significant for a magnetic field as small as 0.1 T for the fine powder samples, whereas even 1 T has little effect on the microwave absorption for the single crystal K3C60 sample. In fact, we observe a giant, about a factor 3 times larger, microwave absorption below Tc than in the normal state. This striking difference between the crystal and fine grain samples is clearly demonstrated for K3C60 where measurements on both kinds of samples are shown. For MgB2, microwave measurements on compacted samples (or surface impedance measurements) supports this observation as therein no enhanced microwave absorption was observed[41-46]. While the absorption becomes significant for the fine powder samples at B = 0.1 T, the shift changes less, which means that the overall superconducting characteristics of the sample are maintained.

Figure 1

Temperature dependent cavity loss, , and cavity shift, Δf/f0 for a fine powder of MgB2 and for the single crystal and powder K3C60 samples. Two magnetic field data are shown for the crystal (0 and 1 T) and three for the powder samples (0, 0.1, and 1 T). Note that the cavity loss changes significantly for the powder sample in contrast to the single crystal sample. Note the different scales for the Δf/f0 data.

We believe that the enhanced microwave absorption is an ubiquitous property of fine powders of type-II superconductors. However, we cannot quantitatively discuss this effect for MgB2 due to the multi-band superconductivity[47,48] and the significant Bc2 anisotropy; Bc2 at 0 K is ~2 T and ~16 T for and , respectively (refs[34,35,49]). We therefore focus on K3C60 in the following. The enhanced microwave loss appears progressively with increasing magnetic field (additional data are shown in the Supplementary Material).

We also show the B = 0.1 T data for the powder sample (B ≈ 0.002 × Bc2) in Fig. 1; they show a peak in the microwave loss right below Tc followed by a gradual decrease. The zero magnetic field data also shows a small peak (invisible at the scale of Fig. 1) for the powder sample (shown in the Supplementary Material). This small peak is not due to magnetic field and is most probably a tiny conductivity coherence peak (the analogue of the Hebel–Slichter peak[50]) which is known to be strongly suppressed by strong-coupling effects in alkali fullerides[51,52]. While the presence of a coherence peak itself is an interesting physical phenomenon[3,53], it is not relevant for the present discussion.

The fact that the enhanced microwave absorption occurs with the application of the magnetic field hints at a flux motion related phenomenon that is discussed in the framework of the CC theory. The microwave absorption peak occurs above the irreversibility line, i.e. it is related to the physical behavior of the vortex-fluid state; for K3C60 Tirr(B = 0.1 T) ≈ 15 K and Tirr(B = 1 T) < 5 K (ref.[54]).

The strong dependence on the sample morphology is also discussed below. Superconducting fullerides are type-II () and have a short mean free path[32] i.e. they can be described in the local electrodynamics limit as opposed to the non-local (or Pippard) limit, which simplifies the discussion.

Discussion

Conductivity in the superconducting state

The phenomenological CC theory[14-19] is based on a two-fluid model and considers the motion of vortices due to the exciting electromagnetic field in the presence of a viscous background (described by the viscous drag coefficient, η) and a restoring force (described by an effective pinning force constant, κp).

The viscous drag was introduced in the Bardeen–Stephen theory[13] and is determined by the superconducting parameters[10,55]: . The value of κp is unknown and only an upper limit can be estimated from thermodynamic considerations[55-57] for a “perfect pinning center”, i.e. a hollow cylinder with a diameter of about the coherence length, d ≈ ξ. The condensation energy gain per unit length from placing a vortex in this cylinder is about , where the square of the thermodynamic critical field is . This leads to:

For a weaker pinning center, κp can be significantly lower and in the bulk of a perfect superconductor, it would be zero.

The CC theory introduces the concept of the complex penetration depth, :where is the normal fluid skin depth, λ is the usual (real) penetration depth and is the complex effective skin depth[19]. The latter quantity is zero for B = 0 and becomes finite in the mixed state when vortex motion is present. is related to the complex conductivity by . Note that for B = 0 (i.e. when ), we obtain and as expected. The CC theory yields the temperature and magnetic field dependent using explicit expressions for λ(T), δnf(T), Bc2(T). We have implemented the calculation (details are given in the Supplementary Material) and validated our calculations by comparing the results to that published in ref.[19].

Here, we discuss qualitatively the predictions of the CC theory and some typical cases are shown in Fig. 2. In superconductors, at B = 0 the carrier spectral weight below ωg collapses into the function according to the Ferrell–Glover–Tinkham (FGT) sum rule[11,12]. The Kramers–Kronig relation dictates that σ2 = 1/μ0ωλ2. Without vortex pinning, the Meissner state is destroyed for B > Bc1 and the Bardeen–Stephen theory gives . It is worth noting that this result is formally analogous to the AC Drude model as the underlying equation of motion (of electrons or vortices) is the same. Here, we introduced a cut-off frequency . Clearly, σ1 can be larger than σn for ω < ωc.

Figure 2

Illustration of in superconductors for i) B = 0, ii) for finite fields (B > Bc1) with κp = 0 (the Bardeen–Stephen case), and iii) for B ≠ 0 and a finite κp (the case of the CC theory). Conductivity above the superconducting gap, ωg, is not shown. The spectral weight in the delta function is preserved for B ≠ 0. Note that for κp ≠ 0, the conductivity appears as if it were a sum of ’s for the B = 0 and the BS flux-flow regimes (shown with dashed curves). Of the two components, the σ1 ∝ δ(ω) and σ2 ∝ 1/ω is due to vortex pinning.

In the presence of vortex pinning, the CC theory predicts that is characterized by a mixture of the unperturbed superconducting behavior and that of the Bardeen–Stephen theory with a shared spectral weight which depends on the pinning force constant. Pinning reduces the effect of the vortex flow on σ1. The enhanced σ1(ω) AC conductivity (as compared to the normal state) is a direct consequence of the FGT sum rule for a finite magnetic field. It allows to estimate the maximum possible value of the enhancement as σ1, max(ω) ≈ σn × ωg/ω, that would be realized at T = 0 in the absence of pinning. E.g. for K3C60 and ω/2π = 10 GHz we obtain σ1, max(ω) ≈ 140σn.

The CC theory allows to quantitatively analyze the conductivity in K3C60. The requirement of is satisfied for our magnetic fields of 0.1 … 1 T as Bc1 ~ 10 mT. The CC theory was developed for a superconductor which occupies the total half space. We show in the Supplementary Material that it can be applied for a spherical sample which approximates well finite sized grains containing at least a few hundred/thousand vortices. In addition, the static and rf magnetic fields are parallel in our experiment, which is the standard case for the applicability of the CC theory. Albeit we cannot quantitatively consider the effect of the small particle size on the magnetic properties, we believe that neither the surface barriers (also known as Bean-Livingstone barriers[58]) nor the so-called geometrical barriers[59] affect considerably the applicability of the CC theory. The argument is that both types of barriers would affect the overall number of the vortices under the applied DC magnetic field (or the B value where vortices appear) but not the overall vortex dynamics under the application of the small AC magnetic field, which is the primary reason for the observed microwave absorption.

In Fig. 3, the calculated conductivity is plotted versus the reduced temperature for different force constants, κp, with the parameters ξ0 = 3 nm (which corresponds to Bc2(0 K) = 37.5 T), λ0 = 440 nm, and ρn(Tc) = 2.95 Ωm, the mean values of the corresponding literature parameters for K3C60 (refs[32,37-39,60-62]), which are detailed in Table 1. The pinning force constant is given in the figure with respect to according to Eq. (1). The figure indicates that κp strongly affects the magnetic field dependence of .

Figure 3

Calculated real part and imaginary part of complex rf conductivity vs the reduced temperature for different values of the pinning force constant, κp. The conductivity values are normalized by the normal state conductivity at the critical temperature. The large value of σ2(T = 0)/σn(Tc) is due to a large (δn/λ)2. Note also the different scales for the σ1 values.

Table 1

Transport and magnetic parameters of the K3C60 superconductor: the superconducting transition temperature, Tc; the normal state resistivity at Tc, ρn; the normal state skin depth, δn; the coherence length at T = 0, ξ0; and the magnetic field penetration depth at T = 0, λ0.

Property	Value	Refs
T c	19.5 K	[32]
ρn(Tc)	1.8⋅10^−6,4.1⋅10^−6 Ωm	[60, 61]
δn(11.1 GHz)	9.7, 6.4 μm
ξ 0	2.6, 3.4 nm	[37, 38]
λ 0	240, 480, 600 nm	[38, 39, 62]

The tabulated ξ0 values correspond to an upper critical field, Bc2 at T = 0 of 49 and 28 T, respectively.

Analysis of the experimental data

The sample morphology greatly affects the relation between the material conductivity, , and the microwave parameters, the loss and shift. Two limiting cases are known. 1) The sample is large and the field penetrates into a limited distance from the surface. This approximates the measurement in the large K3C60 single crystal. 2) The sample is a small sphere with radius comparable to the penetration depth. This approximates the K3C60 sample of well divided small grains. We discuss that the experimental observations for the crystal and fine powder K3C60 are explained well by these two regimes.

In the first case, when the rf field penetrates in the skin depth only (known as the skin limit), the following equation holds between the microwave measurement parameters and the material quantities[28]:where the complex penetration depth, , is related to the conductivity as . The dimensionless is the so-called resonator constant[3] and it depends on the sample surface relative to that of the cavity.

The left panels in Fig. 4. show the calculated and measured cavity loss and shift in 0 and 1 T magnetic fields for the single crystal sample. The calculation uses Eq. (3) with κp = κp,max/20. Although this low κp induces a large σ1, there is no visible peak in the cavity loss in this case when excitation is limited to the surface. We discuss in detail in the Supplementary Material that the calculated cavity loss and shift are insensitive to the value of κp in this limit. Clearly, the experimental data for the K3C60 crystal match well the calculations.

Figure 4

Comparison of measured and calculated cavity loss and shift parameters in the skin (left panels) and penetration limit (right panels). Calculation details are given in the text. Note that the calculated curves and the experimental data agree well for both sample types.

A suitable was chosen to match the calculation to the experiment. We find that for both the calculation and experiment, the cavity loss parameter drops rapidly below Tc, although σ1/σn is around unity due to the vortex motion. This effect is due to the development of a significant σ2/σn ~ 100, which limits the penetration of microwaves into the sample and thus reduces the loss. This means that the microwave surface impedance measurement is not capable of providing information about σ1 in the presence of vortex motion. We note that the experimental curves do not show such a rapid change as a function of temperature as the calculation. This may be related to the finite size and surface roughness of the single crystal sample.

Second, we discuss the opposite limit, when the microwave field penetrates into the sample (known as the penetration limit), the cavity measurables depend differently on the sample parameters. It was shown[63] for a sphere with radius, a:where is the complex wavenumber, with being the complex index of refraction. The dimensionless γ is a sample volume dependent constant.

The right panels in Fig. 4, show the measured cavity loss and shift data for the fine powder sample together with a fit according to Eq. (5). To obtain these fits, we fixed the transport and magnetic parameters (ρn, ξ0, λ0) of K3C60 to the respective mean values as given in Table 1. We assumed that the sample consists of spheres with a uniform diameter, a. The zero magnetic field data depends only on γ and a when the other parameters, δn and λ, are fixed. A fit to the B = 0 data yields , and a = 6.2(2) μm. We then proceed to fit the magnetic field dependent data with κp as the only free parameter and we obtain  N/m2, which is about κp,max/20. As shown in Fig. 4, the calculation agrees well with the experimental data. The presence of a finite microwave absorption at B = 1 T down to the lowest temperature appears to be surprising; its origin is therefore discussed qualitatively. Although the penetration depth due to superconductivity becomes smaller than the particle size, it is accompanied by a large (σ1/σn ≈ 30) conductivity due to the vortex motion. It means that a substantial microwave absorption occurs on the surface of the sample. The important observation is that the loss remains proportional to σ1 in this case, although its effect is reduced. This is further elaborated in the last part of the Supplementary Material.

Somewhat better fits could be obtained when letting λ, ρn, and Bc2 differ from the mean literature values. In addition, Eq. (5) is valid for spheres only, it thus fixes the ratio between the real and imaginary parts (cavity loss and shift). A different particle shape or particle size distribution would allow for a different scaling factor for the loss and shift data which could also improve the fits. Although improved fits could be attained, we believe that the simplest model explains well the experimental observation of an enhanced microwave absorption. In addition it allows to determine an effective pinning force constant, which is an important parameter to describe the electrodynamics of type-II superconductors. However, we note that κp determined herein may overestimate the bulk pinning force constant; it is known that the presence of a substantial surface-volume ratio may give rise to additional vortex pinning[58,59], with a strength that is difficult to estimate.

Summary

We demonstrated that moderate magnetic fields, which are small compared to the upper critical field, induces a large microwave absorption in fine powders of type II superconductors, like MgB2 and K3C60. The effect is absent for samples containing larger grains or compacted powder pellets. The Bardeen–Stephen model of flux-flow predicts that the real part of the AC conductivity can be enhanced in the microwave range, but this effect has not been observed. We analyze the conductivity using the Coffey–Clem theory which also accounts for vortex pinning effects. It is applied to calculate the microwave properties for two kinds of samples: when the electromagnetic field penetration is limited to the surface (skin limit) or when it fully penetrates into the fine grain samples (penetration limit). We show that microwave absorption in the skin limit is little affected by the vortex-motion enhanced σ1 but in the penetration limit, the effect is clearly observable. A quantitative analysis for K3C60 yields the vortex pinning force constant that can be hardly determined by other means. Our observation allowed us to explain long-standing microwave anomalies in superconductors[30,31] and it may lead to pertinent applications in microwave communication techniques.