6.4: Electrodynamics of Superconductivity, and the Gauge Invariance
( \newcommand{\kernel}{\mathrm{null}\,}\)
The effect of superconductivity20 takes place (in certain materials only, mostly metals) when temperature T is reduced below a certain critical temperature Tc, specific for each material. For most metallic superconductors, Tc is of the order of typically a few kelvins, though several compounds (the so-called high-temperature superconductors) with Tc above 100 K have been found since 1987. The most notable property of superconductors is the absence, at T<Tc, of measurable resistance to (not very high) dc currents. However, the electromagnetic properties of superconductors cannot be described by just taking σ=∞ in our previous results. Indeed, for this case, Eq. (33) would give δs=0, i.e., no ac magnetic field penetration at all, while for the dc field we would have the uncertainty σω→? Experiment shows something substantially different: weak magnetic fields do penetrate into superconductors by a material-specific distance δL∼10−7−10−6 m, the so-called London’s penetration depth,21 which is virtually frequency-independent until the skin depth δs, measured in the same material in its “normal” state, i.e. the absence of superconductivity, becomes less than δL. (This crossover happens typically at frequencies ω∼1013−1014 s−1.) The smallness of δL on the human scale means that the magnetic field is pushed out of macroscopic samples at their transition into the superconducting state.
This Meissner-Ochsenfeld effect, discovered experimentally in 1933,22 may be partly understood using the following classical reasoning. The discussion of the Ohm law in Sec. 4.2 implied that the current’s (and hence the electric field’s) frequency ω is either zero or sufficiently low. In the classical Drude reasoning, this is acceptable while ωτ<<1, where τ is the effective carrier scattering time participating in Eqs. (4.12)-(4.13). If this condition is not satisfied, we should take into account the charge carrier inertia; moreover, in the opposite limit ωτ>>1, we may neglect the scattering at all. Classically, we can describe the charge carriers in such a “perfect conductor” as particles with a non-zero mass m, which are accelerated by the electric field, following the 2nd Newton law (4.11),
m˙v=F=qE,
so that the current density j=qnv they create, changes in time as
˙j=qn˙v=q2nmE.
In terms of the Fourier amplitudes of the functions j(t) and E(t), this means
−iωjω=q2nmEω.
Comparing this formula with the relation jω=σEω implied in the last section, we see that we can use all its results with the following replacement:
σ→iq2nmω.
This change replaces the characteristic equation (29) with
−iω=κ2mωiq2nμ, i.e. κ2=μq2nm,
i.e. replaces the skin effect with the field penetration by the following frequency-independent depth:
δ≡1κ=(mμq2n)1/2.
Superficially, this means that the field decay into the superconductor does not depend on frequency:
H(x,t)=H(0,t)e−x/δ,
thus explaining the Meissner-Ochsenfeld effect.
However, there are two problems with this result. First, for the parameters typical for good metals (q=−e,n∼1029 m−3,m∼me,μ≈μ0), Eq. (44) gives δ∼10−8 m, one or two orders of magnitude lower than the experimental values of δL. Experiment also shows that the penetration depth diverges at T→Tc, which is not predicted by Eq. (44).
The second, much more fundamental problem with Eq. (44) is that it has been derived for ωτ>>1. Even if we assume that somehow there is no scattering at all, i.e. τ=∞, at ω→0 both parts of the characteristic equation (43) vanish, and we cannot make any conclusion about κ. This is not just a mathematical artifact we could ignore. For example, let us place a non-magnetic metal into a static external magnetic field at T>Tc. The field would completely penetrate the sample. Now let us cool it. As soon as temperature is decreased below Tc, the above calculations would become valid, forbidding the penetration into the superconductor of any change of the field, so that the initial field would be “frozen” inside the sample. The experiment shows something completely different: as T is lowered below Tc, the initial field is being pushed out of the sample.
The resolution of these contradictions has been provided by quantum mechanics. As was explained in 1957 in a seminal work by J. Bardeen, L. Cooper, and J. Schrieffer (commonly referred to as the BSC theory), the superconductivity is due to the correlated motion of electron pairs, with opposite spins and nearly opposite momenta. Such Cooper pairs, each with the electric charge q=−2e and zero spin, may form only in a narrow energy layer near the Fermi surface, of certain thickness Δ(T). This
parameter Δ(T), which may be also considered as the binding energy of the pair, tends to zero at T→Tc, while at T<<Tc it has a virtually constant value Δ(0)≈3.5kBTc, of the order of a few meV for most superconductors. This fact readily explains the relatively low spatial density of the Cooper pairs: np∼nΔ(T)/εF∼1026 m−3. With the correction n→np, Eq. (44) for the penetration depth becomes
δ→δL=(mμq2np)1/2.London’s penetration depth
This result diverges at T→Tc, and generally fits the experimental data reasonably well, at least for the so-called “clean” superconductors (with the mean free path l=νFτ, where νF∼(2mεF)1/2 is the r.m.s. velocity of electrons on the Fermi surface, much longer than the Cooper pair size ξ – see below).
The smallness of the coupling energy Δ(T) is also a key factor in the explanation of the Meissner-Ochsenfeld effect, as well as several macroscopic quantum phenomena in superconductors. Because of Heisenberg’s quantum uncertainty relation δrδp∼ℏ, the spatial extension of the Cooper-pair’s wavefunction (the so-called coherence length of the superconductor) is relatively large: ξ∼δr∼ℏ/δp∼ℏνF/Δ(T)∼10−6 m. As a result, npξ3>>1, meaning that the wavefunctions of the pairs are strongly overlapped in space. Due to their integer spin, Cooper pairs behave like bosons, which means in particular that at low temperatures they exhibit the so-called Bose-Einstein condensation onto the same ground energy level εg.23 This means that the frequency ω=εg/ℏ of the time evolution of each pair’s wavefunction Ψ=ψexp{−iωt} is exactly the same, and that the phases φ of the wavefunctions, defined by the relation
ψ=|ψ|eiφ,
are equal, so that the electric current is carried not by individual Cooper pairs but rather their Bose-Einstein condensate described by a single wavefunction (47). Due to this coherence, the quantum effects (which are, in the usual Fermi-gases of single electrons, masked by the statistical spread of their energies, and hence of their phases), become very explicit – “macroscopic”.
To illustrate this, let us write the well-known quantum-mechanical formula for the probability current density of a free, non-relativistic particle,24
jw=iℏ2m(ψ∇ψ∗− c.c. )≡12m[ψ∗(−iℏ∇)ψ− c.c. ],
where c.c. means the complex conjugate of the previous expression. Now let me borrow one result that will be proved later in this course (in Sec. 9.7) when we discuss the analytical mechanics of a charged particle moving in an electromagnetic field. Namely, to account for the magnetic field effects, the particle’s kinetic momentum p≡mv (where v≡dr/dt is particle’s velocity) has to be distinguished from its canonical momentum,25
P≡p+qA.
where A is the field’s vector potential defined by Eq. (5.27). In contrast with the Cartesian components pj=mνj of the kinetic momentum p, the canonical momentum’s components are the generalized momenta corresponding to the Cartesian components rj of the radius-vector r, considered as generalized coordinates of the particle: Pj=∂L/∂νj, where L is the particle’s Lagrangian function. According to the general rules of transfer from classical to quantum mechanics,26 it is the vector P whose operator (in the coordinate representation) equals −iℏ∇, so that the operator of the kinetic momentum p=P−qA is equal to −iℏ∇qA. Hence, to account for the magnetic field27 effects, we should make the following replacement,
−iℏ∇→−iℏ∇−qA,
in all field-free quantum-mechanical relations. In particular, Eq. (48) has to be generalized as
jw=12m[ψ∗(−iℏ∇−qA)ψ− c.c. ].
This expression becomes more transparent if we take the wavefunction in form (47):
jw=ℏm|ψ|2(∇φ−qℏA).
This relation means, in particular, that in order to keep jw gauge-invariant, the transformation (8)-(9) has to be accompanied by a simultaneous transformation of the wavefunction’s phase:
φ→φ+qℏχ.
It is fascinating that the quantum-mechanical wavefunction (or more exactly, its phase) is not gauge-invariant, meaning that you may change it in your mind – at your free will! Again, this does not change any observable (such as jw or the probability density ψψ∗, i.e. any experimental results.
Now for the electric current density of the whole superconducting condensate, Eq. (52) yields the following constitutive relation:
j≡jwqnp=ℏqnpm|ψ|2(∇φ−qℏA)Supercurrent density
The formula shows that this supercurrent may be induced by the dc magnetic field alone and does not require any electric field. Indeed, for the simplest, 1D geometry shown in Fig.2a, j(r)=j(x)nz, A(r)=A(x)nz, and ∂/∂z=0, so that the Coulomb gauge condition (5.48) is satisfied for any choice of the gauge function χ(x), and for the sake of simplicity we can choose it to provide φ(r)≡const,28 so that
j=−q2npmA≡−1μδ2LA.
where δL is given by Eq. (46), and the field is assumed to be small, and hence not affecting the probability |ψ|2 (normalized to 1 in the absence of the field). This is the so-called London equation, proposed (in a different form) by F. and H. London in 1935 for the Meissner-Ochsenfeld effect’s explanation. Combining it with Eq. (5.44), generalized for a linear magnetic medium by the replacement μ0→μ, we get
∇2A=1δ2LA,
This simple differential equation, similar to Eq. (23), for our 1D geometry has an exponential solution similar to Eq. (32):
A(x)=A(0)exp{−xδL},B(x)=B(0)exp{−xδL},j(x)=j(0)exp{−xδL},
which shows that the magnetic field and the supercurrent penetrate into a superconductor only by London’s penetration depth δL, regardless of frequency.29 By the way, integrating the last result through the penetration layer, and using the vector potential’s definition, B=∇×A (for our geometry, giving B(x)=dA(x)/dx=−δLA(x)) we may readily verify that the linear density J of the surface supercurrent still satisfies the universal coarse-grain relation (38).
This universality should bring to our attention the following common feature of the skin effect (in “normal” conductors) and the Meissner-Ochsenfeld effect (in superconductors): if the linear size of a bulk sample is much larger than, respectively, δs or δL, than B=0 in the dominating part of its interior. According to Eq. (5.110), a formal description of such conductors (valid only on a coarse-grain scale much larger than either δs or δL), may be achieved by formally treating the sample as an ideal diamagnet, with μ=0. In particular, we can use this description and Eq. (5.124) to immediately obtain the magnetic field’s distribution outside of a bulk sphere:
B=μ0H=−μ0∇ϕm, with ϕm=H0(−r−R32r2)cosθ, for r≥R.
Figure 3 shows the corresponding surfaces of equal potential ϕm. It is evident that the magnetic field lines (which are normal to the equipotential surfaces) bend to become parallel to the surface near it.

This pattern also helps to answer the question that might arise at making the assumption (24): what happens to bulk conductors placed into in a normal ac magnetic field – and to superconductors in a normal dc magnetic field as well? The answer is: the field is deformed outside of the conductor to sustain the following coarse-grain boundary condition:30
Coarse-grain boundary conditionBn|surface =0,
which follows from Eq. (5.118) and the coarse-grain requirement B|inside =0.
This answer should be taken with reservations. For normal conductors it is only valid at sufficiently high frequencies where the skin depth (33) is sufficiently small: δs<<a, where a is the scale of the conductor’s linear size – for a sphere, a∼R. In superconductors, this simple picture is valid not
only if δs<<a, but also only in sufficiently low magnetic fields, because strong fields do penetrate into superconductors, destroying superconductivity (completely or partly), and as a result violating the Meissner-Ochsenfeld effect – see the next section.
Reference
20 Discovered experimentally in 1911 by Heike Kamerlingh Onnes.
21 Named so to acknowledge the pioneering theoretical work of brothers Fritz and Heinz London – see below.
22 It is hardly fair to shorten the name to just the “Meissner effect”, as it is frequently done, because of the reportedly crucial contribution by Robert Ochsenfeld, then a Walther Meissner’s student, into the discovery.
23 A quantitative discussion of the Bose-Einstein condensation of bosons may be found in SM Sec. 3.4, though the full theory of superconductivity is more complicated because it has to describe the condensation taking place simultaneously with the formation of effective bosons (Cooper pairs) from fermions (single electrons). For a detailed, but very readable coverage of the physics of superconductors, I can refer the reader to the monograph by M. Tinkham, Introduction to Superconductivity, 2nd ed., McGraw-Hill, 1996.
24 See, e.g., QM Sec. 1.4, in particular Eq. (1.47).
25 I am sorry to use traditional notations p and P for the momenta – the same symbols which were used for the electric dipole moment and polarization in Chapter 3. I hope there will be no confusion, because the latter notions are not used in this section.
26 See, e.g., CM Sec. 10.1, in particular Eq. (10.26).
27 The account of the electric field is easier, because the related energy qϕ of the particle may be directly included in the potential energy operator.
28 This is the so-called London gauge; for our simple geometry, it is also the Coulomb gauge (5.48).
29 Since at T>0 not all electrons in a superconductor form Cooper pairs, at any frequency ω≠0 the unpaired electrons provide energy-dissipating Ohmic currents, which are not described by Eq. (54). These losses become very substantial when the frequency ω becomes so high that the skin-effect length δs of the material (as measured with superconductivity suppressed, say by high magnetic field) becomes less than δL. For typical metallic superconductors, this crossover takes place at frequencies of a few hundred GHz, so that even for microwaves, Eq. (57) still gives a fairly accurate description of the field penetration.
30 Sometimes this boundary condition, as well as the (compatible) Eq. (38), are called “macroscopic”. However, this term may lead to confusion with the genuine macroscopic boundary conditions (5.117)-(5.118), which also ignore the atomic-scale microstructure of the “effective currents” jef=∇×M, but (as was shown earlier in this section) still allow explicit, detailed accounts of the skin-current (34) and supercurrent (55) distributions.