Skip to content

Superconductivity

Superconductivity is not perfect conductivity. The distinction is the whole subject, and it is the first thing to get straight: a perfect conductor and a superconductor respond differently to a magnetic field, and only one of them expels it. Everything in this page is downstream of that single experimental fact.

Picture current in a copper wire as electrons scattering off vibrations and impurities, losing a little energy each time — that is resistance. Now picture the electrons paired up and moving as a single coherent object, described by one collective wavefunction that extends through the whole piece of metal. To scatter, an electron would have to break its pair, and at low temperature there is not enough thermal energy to do it. Resistance does not become small; it becomes exactly zero, and current set flowing in a superconducting ring has been observed to persist for years without decay.

The Meissner effect is stranger. Bring a magnet near a superconductor and the superconductor pushes back, levitating the magnet. That is not simply “zero resistance” doing its work: a perfect conductor would freeze in whatever field was present when it cooled, because with ρ=0\rho = 0 the flux through a loop cannot change. A superconductor does the opposite — it expels the field regardless of history. Expulsion is a thermodynamic statement about the ground state, and that is why the London equations, which describe it, are not a consequence of ρ=0\rho = 0.

Discovered by Onnes in 1911, in mercury at Tc=4.2 KT_c = 4.2\,\mathrm{K}.

  1. Zero resistance. ρ=0\rho = 0 for T<TcT < T_c.
  2. Meissner effect. Magnetic flux is expelled from the interior: B=0\mathbf{B} = 0 inside for T<TcT < T_c and B<BcB < B_c.
  3. Critical field. Superconductivity is destroyed above Bc(T)=Bc(0)[1−(T/Tc)2]B_c(T) = B_c(0)\big[1 - (T/T_c)^2\big].
  4. Critical current. Destroyed above a critical current density JcJ_c.

The critical field is thermodynamic: it is where the free energy gained by condensing into the superconducting state no longer pays for expelling the field. The condensation energy density is

12μ0Hc2=Fn−FsV,\frac{1}{2}\mu_0 H_c^2 = \frac{F_n - F_s}{V},

the difference between normal and superconducting free energies per unit volume, which is why Bc(T)B_c(T) has a definite shape and is not a fitting curve.

The London brothers wrote the phenomenology down as two equations for the supercurrent density Js\mathbf{J}_s, with nsn_s the density of superconducting electrons:

∂Js∂t=nse2me E,∇×Js=−nse2me B.\frac{\partial \mathbf{J}_s}{\partial t} = \frac{n_s e^2}{m_e}\,\mathbf{E}, \qquad \nabla \times \mathbf{J}_s = -\frac{n_s e^2}{m_e}\,\mathbf{B}.

The first is Newton’s second law for carriers that never scatter: me ∂v/∂t=−eEm_e\,\partial\mathbf{v}/\partial t = -e\mathbf{E} with Js=−nsev\mathbf{J}_s = -n_s e\mathbf{v}. The second is the one that carries new physics, and it is assumed, not derived from ρ=0\rho = 0. To see that the two are genuinely different, take the first London equation and Maxwell’s equation for a conductor with σ→∞\sigma \to \infty; what follows is that ∂B/∂t=0\partial\mathbf{B}/\partial t = 0 — flux is frozen, whatever it was when the material went superconducting. The second London equation says B=0\mathbf{B} = 0 outright, which is a different claim and the one experiment supports.

Deriving the penetration depth. Take ∇×\nabla \times of the second London equation and use Ampere’s law ∇×B=μ0Js\nabla \times \mathbf{B} = \mu_0 \mathbf{J}_s (neglecting the normal current):

∇×(∇×Js)=−nse2me ∇×B=−nse2me μ0Js.\nabla \times (\nabla \times \mathbf{J}_s) = -\frac{n_s e^2}{m_e}\,\nabla \times \mathbf{B} = -\frac{n_s e^2}{m_e}\,\mu_0 \mathbf{J}_s.

Using ∇×(∇×Js)=∇(∇⋅Js)−∇2Js\nabla \times (\nabla \times \mathbf{J}_s) = \nabla(\nabla \cdot \mathbf{J}_s) - \nabla^2\mathbf{J}_s and ∇⋅Js=0\nabla \cdot \mathbf{J}_s = 0 for steady state:

∇2Js=1λL2 Js,λL=meμ0nse2,\nabla^2 \mathbf{J}_s = \frac{1}{\lambda_L^2}\,\mathbf{J}_s, \qquad \lambda_L = \sqrt{\frac{m_e}{\mu_0 n_s e^2}},

and identically for B\mathbf{B}:

∇2B=BλL2.\nabla^2 \mathbf{B} = \frac{\mathbf{B}}{\lambda_L^2}.

For a half-space x>0x > 0 with field B0B_0 outside, the decaying solution is

B(x)=B0 e−x/λL,\mathbf{B}(x) = B_0\,e^{-x/\lambda_L},

so the field falls by 1/e1/e over the London penetration depth λL\lambda_L, typically 1010—100 nm100\,\mathrm{nm}. The Meissner effect is not a surface with field exactly zero beyond it; it is field decaying over a finite depth, and B=0\mathbf{B} = 0 in the bulk is the large-body limit of that decay.

The London equations describe the electromagnetic response but say nothing about why there is a superconducting order parameter, or what happens near TcT_c. Ginzburg-Landau theory introduces a complex order parameter ψ(r)\psi(\mathbf{r}) with ∣ψ∣2=ns|\psi|^2 = n_s and writes the free energy as

F=Fn+∫d3r[α∣ψ∣2+β2∣ψ∣4+12m∗∣(−iℏ∇−2eA)ψ∣2+∣B∣22μ0].F = F_n + \int d^3r \left[ \alpha|\psi|^2 + \frac{\beta}{2}|\psi|^4 + \frac{1}{2m^*}\left|(-i\hbar\nabla - 2e\mathbf{A})\psi\right|^2 + \frac{|\mathbf{B}|^2}{2\mu_0} \right].

Minimising gives two equations. Varying in ψ∗\psi^*:

αψ+β∣ψ∣2ψ+12m∗(−iℏ∇−2eA)2ψ=0,\alpha\psi + \beta|\psi|^2\psi + \frac{1}{2m^*}(-i\hbar\nabla - 2e\mathbf{A})^2\psi = 0,

which is a Schrödinger equation for ψ\psi with charge 2e2e — the first hint of pairing, written down in 1950, seven years before BCS explained it. Varying in A\mathbf{A} gives

Js=e∗ℏ2m∗(ψ∗∇ψ−ψ∇ψ∗)−e∗2m∗∣ψ∣2A,e∗=2e,\mathbf{J}_s = \frac{e^*\hbar}{2m^*}\left(\psi^*\nabla\psi - \psi\nabla\psi^*\right) - \frac{e^{*2}}{m^*}|\psi|^2\mathbf{A}, \qquad e^* = 2e,

which reproduces both London equations in the limit ∣ψ∣|\psi| constant, and additionally shows the flux through a ring is quantised in units of h/2eh/2e — the other hint.

Two lengths fall out of the theory, and their ratio controls everything:

ξ=ℏ22m∗∣α∣,κ=λLξ.\xi = \sqrt{\frac{\hbar^2}{2m^*|\alpha|}}, \qquad \kappa = \frac{\lambda_L}{\xi}.

ξ\xi is the coherence length, the smallest distance over which ∣ψ∣|\psi| can vary — the size of a Cooper pair. κ\kappa is the GL parameter, and the single criterion

κ<12  type I,κ>12  type II\kappa < \frac{1}{\sqrt{2}} \;\text{type I}, \qquad \kappa > \frac{1}{\sqrt{2}} \;\text{type II}

is the whole of the type I/II distinction. That a material’s response to a field is decided by one dimensionless ratio is the theory’s central result.

BCS theory (Bardeen, Cooper, Schrieffer, 1957) explains the order parameter. The mechanism is Cooper pairing: two electrons just above the Fermi surface with opposite momenta and spins form a bound state through the lattice-mediated effective attraction. The pair has charge 2e2e and spin 00, so it is a boson and can condense.

The gap equation. With Ek=ξk2+Δ2E_{\mathbf{k}} = \sqrt{\xi_{\mathbf{k}}^2 + \Delta^2} the quasiparticle energy, ξk\xi_{\mathbf{k}} the normal-state energy measured from EFE_F, and VpairV_{\mathrm{pair}} the effective attraction, self-consistency requires

Δ=Vpair∑kΔ2Ektanh⁡ ⁣(Ek2kBT).\Delta = V_{\mathrm{pair}} \sum_{\mathbf{k}} \frac{\Delta}{2E_{\mathbf{k}}} \tanh\!\left(\frac{E_{\mathbf{k}}}{2k_B T}\right).

Dividing out Δ\Delta and replacing the sum by an integral over the density of states N(EF)N(E_F) (constant near the Fermi level, over the Debye window ℏωD\hbar\omega_D) gives

1=N(EF)Vpair∫0ℏωDdξξ2+Δ2tanh⁡ ⁣(ξ2+Δ22kBT).1 = N(E_F) V_{\mathrm{pair}} \int_0^{\hbar\omega_D} \frac{d\xi}{\sqrt{\xi^2 + \Delta^2}} \tanh\!\left(\frac{\sqrt{\xi^2+\Delta^2}}{2k_B T}\right).

At T=0T = 0 the tanh⁡\tanh is 11 and the integral is elementary:

∫0ℏωDdξξ2+Δ2=arsinh⁡ ⁣(ℏωDΔ)≈ln⁡ ⁣2ℏωDΔ,\int_0^{\hbar\omega_D} \frac{d\xi}{\sqrt{\xi^2 + \Delta^2}} = \operatorname{arsinh}\!\left(\frac{\hbar\omega_D}{\Delta}\right) \approx \ln\!\frac{2\hbar\omega_D}{\Delta},

so

Δ(0)=2ℏωD e−1/N(EF)Vpair.\Delta(0) = 2\hbar\omega_D\, e^{-1/N(E_F)V_{\mathrm{pair}}}.

At Δ→0\Delta \to 0, setting T=TcT = T_c:

kBTc=1.13 ℏωD e−1/N(EF)Vpair.k_B T_c = 1.13\,\hbar\omega_D\, e^{-1/N(E_F)V_{\mathrm{pair}}}.

Dividing the two eliminates the unknown coupling entirely:

2Δ(0)kBTc=4ℏωD1.13 ℏωD⋅12⋅2≈3.53,\frac{2\Delta(0)}{k_B T_c} = \frac{4\hbar\omega_D}{1.13\,\hbar\omega_D} \cdot \frac{1}{2} \cdot 2 \approx 3.53,

a universal number with no free parameter, and the cleanest test of the theory. Measured values in conventional superconductors cluster around 3.53.5.

The exponential in −1/(NV)-1/(N V) is the physical content. A weak attraction produces an exponentially small gap, which is why superconductivity was not understood for 46 years: perturbation theory in VV produces a series in NVNV, and the answer is not analytic at NV=0NV = 0. Every term of the perturbation series is finite and they sum to zero.

Type IType II
GL parameterκ<1/2\kappa < 1/\sqrt{2}κ>1/2\kappa > 1/\sqrt{2}
Critical fieldone, BcB_ctwo, Bc1<Bc2B_{c1} < B_{c2}
Intermediate statenonemixed, Bc1<B<Bc2B_{c1} < B < B_{c2}
ExamplesPb, Hg, AlNb, YBCO, all high-TcT_c

In the mixed state, magnetic flux penetrates as vortices — tubes of normal material carrying one flux quantum Φ0=h/2e\Phi_0 = h/2e each, with supercurrent circulating around them. Vortices repel, so they form a lattice; pinning them to defects is what lets a type II material carry large currents, which is why superconducting magnets are built from type II materials and why JcJ_c is an engineering quantity rather than a fundamental one.

Flux quantisation in units of h/2eh/2e rather than h/eh/e was the experimental confirmation that the carriers have charge 2e2e, and it is independent evidence for pairing.

Two superconductors separated by a thin insulating barrier exchange a supercurrent. The DC Josephson effect: with no applied voltage, a current

J=Jcsin⁡δJ = J_c \sin\delta

flows, where δ=θ1−θ2\delta = \theta_1 - \theta_2 is the phase difference across the junction and JcJ_c depends on the barrier. Current flows with zero voltage — not a small resistance, no resistance.

The AC Josephson effect: with a DC voltage VV across the junction, the phase evolves as

dδdt=2eVℏ,\frac{d\delta}{dt} = \frac{2eV}{\hbar},

so the current oscillates at frequency ω=2eV/ℏ\omega = 2eV/\hbar. Since ee and ℏ\hbar are the two most precisely known constants in physics, this turns a voltage measurement into a frequency measurement and is the basis of the volt standard.

Both effects follow from one assumption: that ψ\psi is a coherent order parameter with a well-defined phase. Neither has a normal-state analogue, which is why the Josephson junction is the cleanest probe of superconducting phase coherence.

Discovered in 1986 by Bednorz and Müller. Cuprates such as YBCO (YBa2Cu3O7−δ\mathrm{YBa_2Cu_3O_{7-\delta}}) reach Tc≈135 KT_c \approx 135\,\mathrm{K}, above liquid nitrogen’s 77 K77\,\mathrm{K} — the reason they transformed applications.

They are type II with large κ\kappa, layered, and not explained by BCS. Three differences matter:

  • dd-wave pairing. The gap has dx2−y2d_{x^2-y^2} symmetry, Δ(k)=Δ0(cos⁡kx−cos⁡ky)/2\Delta(\mathbf{k}) = \Delta_0(\cos k_x - \cos k_y)/2, which changes sign under a 90°90° rotation and vanishes along the diagonals kx=±kyk_x = \pm k_y. Conventional superconductors are ss-wave, with Δ\Delta constant. The sign change is measurable in phase-sensitive experiments and is established.
  • Short coherence length. ξ∼1\xi \sim 1—2 nm2\,\mathrm{nm} against ∼100 nm\sim 100\,\mathrm{nm} conventionally, so κ\kappa is large and the materials are extreme type II. It also makes them sensitive to defects at the scale of a unit cell.
  • Strong anisotropy. Properties differ sharply between the abab-planes and the cc-axis, a consequence of the layered perovskite structure.

The pairing mechanism is still open. Phonons cannot be the whole story: the isotope effect, which in conventional superconductors confirms lattice involvement through Tc∝M−1/2T_c \propto M^{-1/2}, is much reduced or absent in cuprates. Candidate mechanisms involve antiferromagnetic spin fluctuations, and no consensus has formed.

  • Calling it perfect conductivity. A perfect conductor freezes in flux; a superconductor expels it. Only the second is observed, and only the London equations give it.
  • Treating B=0\mathbf{B} = 0 as exact everywhere. The field decays over λL\lambda_L, which is finite; the bulk statement is a limit for samples much larger than λL\lambda_L.
  • Reading the gap equation as an equation to solve for VV. It is a self-consistency condition on Δ\Delta: Δ\Delta appears on both sides, and Δ=0\Delta = 0 is always a solution. Superconductivity is the non-trivial solution, and it exists only below TcT_c.
  • Expecting perturbation theory to find the gap. The answer is exp⁡(−1/NV)\exp(-1/NV), which is non-analytic at NV=0NV = 0: every term of the series in NVNV vanishes. This is why the problem sat open for 46 years.
  • Quoting 2Δ(0)/kBTc=3.532\Delta(0)/k_BT_c = 3.53 as universal. It is universal within weak-coupling BCS. Strong coupling, and all the cuprates, deviate — and the deviation is diagnostic.
  • Using κ\kappa as a small correction. κ\kappa decides type I versus type II outright at κ=1/2\kappa = 1/\sqrt{2}, which decides whether a material is useful for magnets.