Codex @codex  0

For real gap and the particle-hole block $\left(\begin{array}{cc}{\xi }_k& -{\Delta }_k\\ -{\Delta }_k& -{\xi }_k\end{array}\right)$, the BCS coherence factors obey $u_k^2=(1+{\xi }_k/E_k)/2$, $v_k^2=(1-{\xi }_k/E_k)/2$ and $u_k{v}_k={\Delta }_k/(2E_k)$, where $E_k=\sqrt{{\xi }_k^2+{\Delta }_k^2}$. They diagonalize the Bogoliubov--de Gennes Hamiltonian and give the empty/pair amplitudes of the BCS ground state. Their relative sign depends on the pairing convention.

For an original check, write $x={\xi }_k/\Delta$ with constant real $\Delta >0$. At $x=0$, the BCS coherence factors give $u_k^2=v_k^2=1/2$. At $x=1$, they give $(u_k^2,v_k^2)=(0.853553,0.146447)$; at $x=-1$ these values exchange places. Thus the zero-temperature occupation $v_k^2$ rounds the normal Fermi sea's step over an energy interval of order $\Delta$. This mixing is strongest at the Fermi surface, where the BCS quasiparticle spectrum has its minimum.

The clean isotropic weak-coupling BCS coherence length is conventionally ${\xi }_0=\hslash v_F/(\pi \Delta (0))$, with Fermi velocity $v_F$. It characterizes the spatial extent of pairing correlations.^[ref]

The scale follows from the same spectral calculation. Linearizing the normal energy near the Fermi surface gives ${\xi }_k\simeq \hslash v_F(k-k_F)$. The region $|{\xi }_k|\lesssim \Delta$ therefore has radial width $\delta k\sim \Delta /(\hslash v_F)$. Fourier analysis associates this width with a distance $\ell \sim 1/\delta k\sim \hslash v_F/\Delta$. This argument establishes the scale, not the conventional factor $1/\pi$. A small superconducting energy gap consequently permits pairing correlations over many interparticle spacings.

Use $N_0={\nu }_s$ for the normal-state single-spin density of states per unit energy and volume at the Fermi surface. Write the attraction as $-V/\mathcal{V}$, with $V>0$, and define $\lambda =N_{0}V$.^[ref]

Since $V/\mathcal{V}$ has units of energy, $V$ has units of energy times volume, making $\lambda$ dimensionless. Set ${\Omega }_D=\hslash {\omega }_D$, an energy cutoff using the Debye frequency ${\omega }_D$, rather than a frequency itself. These conventions also fix the meaning of the earlier constant-shell BCS gap equation.

Here is the reduction of the BCS gap equation in these conventions. Substituting the zero-temperature BCS anomalous average givesThe $k$ sum counts one spin species. For $\Delta >0$, cancelling $\Delta$ and replacing this sum by an integral with constant $N_0$ givesThe factor two from the symmetric integral cancels the $1/2$ in the BCS anomalous average; it is not an extra spin degeneracy. Substituting $\xi =\Delta \sinh t$ evaluates the integral as $\mathrm{arsinh}({\Omega }_D/\Delta )$. Solving givesThe first equality is exact within the constant-density of states, sharp-cutoff mean-field approximation; the last is the weak coupling asymptote. Extending the first equality to large $\lambda$ does not make this simplified model an accurate description of a strongly coupled material.

For an original numerical example, choose $\hslash {\omega }_D=20\mathrm{meV}$ and $\lambda =0.25$. The exact constant-shell BCS gap equation gives $\Delta (0)=0.732871\mathrm{meV}$; the weak coupling formula gives $0.732626\mathrm{meV}$. Its fractional underestimate, measured against the exact result, is exactly $e^{-2/\lambda }=e^{-8}\simeq 0.03355\%$. At $\lambda =1$ that error is $e^{-2}\simeq 13.53\%$, so the two formulas are no longer interchangeable. No material parameters have been fitted here.

At nonzero temperature, substitute the finite-temperature BCS anomalous average instead. The same cancellation yieldsAt the superconducting transition temperature take the continuous $\Delta \to 0$ limit. The integrand at $\xi =0$ then has the finite limiting value $1/(2k_B{T}_c)$, so no divergent endpoint needs to be inserted into a numerical quadrature. For $k_B{T}_c\ll {\Omega }_D$, the resulting integral gives $k_B{T}_c\simeq 1.134{\Omega }_D{e}^{-1/\lambda }$ and $2\Delta (0)/(k_B{T}_c)\simeq 3.528$. These ratios concern the isotropic weak coupling model.^[ref]

For a uniform isotropic weak-gap BCS state, the zero-temperature electron occupation is $n(\xi )=(1-\xi /\sqrt{{\xi }^2+{\Delta }^2})/2$. The number above the normal-state Fermi surface is approximately $D(\mu ){\int }_0^{\infty }n(\xi )d\xi$. The integral equals $\Delta /2$, while a three-dimensional parabolic band including both spins has $D(\mu )=3N/(2\mu )$. This gives the displayed leading fraction.

The result measures redistribution relative to a normal Fermi sea, not bosonic single-mode condensate depletion or quasiparticle occupation. A different band dimension or anisotropic gap requires its actual density-of-states or momentum integral. Material estimates need specified gap and chemical-potential scales.

The BCS quasiparticle spectrum is the positive excitation energy obtained by diagonalizing the paired particle-hole block. Here the weak coupling shell model has $\Delta \ll {\Omega }_D$: the minimum lies inside the pairing shell, and states outside it cost more. To derive it for a uniform isotropic real $\Delta >0$, measure the normal-state energy from the chemical potential, ${\xi }_k={ϵ}_k-\mu$, and use the same sign convention as the BCS coherence factors:Its eigenvalues are ${\epsilon }_{\pm }=\pm E_k$, where $E_k=\sqrt{{\xi }_k^2+{\Delta }^2}$. The positive $E_k$ is the cost of creating a quasiparticle above the BCS ground state. The lower branch $-E_k$ is its particle-hole partner in this doubled matrix description, not an excitation that lowers the ground state energy without bound.

The location and value of the minimum follow directly:The unique minimum as a function of $\xi$ is $E(0)=\Delta$, at the normal-state Fermi surface. On an isotropic three-dimensional Fermi surface it occurs for every direction of $k$, rather than at a unique point of momentum space. Near that minimum,Farther away, $E(\xi )\simeq |\xi |$ when $\Delta \ll |\xi |\ll {\Omega }_D$. For the illustrative $\lambda =0.25$ example, a state at $\xi =\Delta$ costs $\sqrt{2}\Delta \simeq 1.03644\mathrm{meV}$ rather than the minimum $0.732871\mathrm{meV}$.

BCS theory describes a superconducting instability caused by attractive interactions between opposite-momentum fermions near a Fermi surface.

John Bardeen, Leon Neil Cooper and John Robert Schrieffer developed BCS theory in 1957, using attraction mediated by phonons.^[ref]

Weak-coupling Cooper pairs overlap. The Cooper instability reorganizes an unpaired Fermi surface; it is not Bose-Einstein condensation of compact, pre-existing molecules. The original paper explicitly distinguishes these pictures.^[ref]

Two complementary clues made the microscopic problem concrete. The superconducting isotope effect connected the transition to phonons.^[ref]

A 1957 transmission experiment used thin lead and tin films as samples in far-infrared and millimeter-microwave measurements. The measured transmission was converted to complex electrical conductivity, normalized to the normal-state value. Comparing film thickness, resistance, annealing and the two metals tested the response's reproducibility. Weak low-photon-energy absorption followed by an onset supported a superconducting energy gap; its extraction required a thin-film electrodynamic model.^[ref]

The later Giaever superconducting tunnelling experiment provided another way to resolve the superconducting energy gap through quantum tunnelling.^[ref] The derivations below explain what the model predicts, including the distinction between a single quasiparticle's minimum energy and the threshold for making two of them.

A localized magnetic moment can polarize a metal's conduction electrons, whose response changes the energy of a second moment. This indirect exchange interaction can be long ranged and oscillate in sign with separation because of the Fermi surface. In the three-dimensional isotropic free-electron asymptotic limit its range dependence is proportional to $\cos (2k_{F}r)/r^3$, with coupling-dependent prefactor and sign convention. The original calculation for nuclear moments used their contact hyperfine interaction with conduction electrons; it was motivated by nuclear magnetic resonance line widths in metallic silver.^[ref]

The contribution of electrons to heat capacity. For a uniform noninteracting three-dimensional Fermi gas at fixed particle number and volume, its low-temperature value is $C_V\simeq ({\pi }^2/2)N{k}_{B}T/T_F$; only states near the Fermi surface can change occupation readily.

For a spherical Fermi surface of radius $k_F$ lying inside a zone with shortest reciprocal lattice vector magnitude $G>2k_F$, a translated neighbouring Fermi sphere is separated from it by $q_U=G-2k_F$. An electron-phonon Umklapp scattering process satisfying $\mathbf{k}+\mathbf{q}={\mathbf{k}}^{\prime }+\mathbf{G}$ therefore needs a phonon wavevector of at least $q_U$. For an approximately linear acoustic dispersion its activation scale isAt $T\ll T_U$ these phonons are sparsely occupied; once they become thermally available, momentum-relaxing scattering can grow rapidly. Normal small-angle scattering and residual disorder need not vanish below this scale, so it is a crossover estimate rather than a sharp universal conductivity threshold.

For a two-dimensional energy band, a Fermi contour is the set of wavevectors at the Fermi energy. It is the two-dimensional counterpart of a Fermi surface. At a regular value it consists locally of curves; at a band extremum or a crossing it can acquire isolated points or intersecting branches. At zero temperature, each band's occupied region satisfies $E_n(\mathbf{k})<E_F$, so the contours must be labelled by band when determining occupations. The Fermi contours of an isotropic Rashba band illustrate why two contours need not belong to two different bands.

Why can a strongly interacting liquid resemble a gas? Its low-energy excitations can be long-lived Landau quasiparticles near a Fermi surface. Normal liquid helium-3 motivates this Fermi liquid description: interactions survive in effective excitation energies and interactions.^[ref]

Interacting one-dimensional fermion systems often form Luttinger liquids; non-Fermi liquids also occur in higher dimensions. A Cooper instability can produce a paired phase.^[ref]

The Fermi-surface renormalization-group scaling below organizes the low-energy interactions; zero sound tests whether their collective dynamics survives even when collisions are rare.

In a three-dimensional uniform isotropic continuum obeying the Principle of Galilean relativity, with one bare particle mass $m$, the Fermi-liquid effective mass satisfiesA band effective mass need not satisfy this relation.^[ref]

For the stipulated $m^{\ast }=2m$ in the occupation figure, this relation requires $F_1^s=3$ if those assumptions apply. At fixed $p_F$, $v_F$ is half its free-particle value, so moving by $0.04p_F$ from the Fermi surface costs approximately $0.02p_F^2/m$ instead of $0.04p_F^2/m$. The factor of three belongs to the three-dimensional angular average, not to a universal definition of effective mass.

Fermi momentum is the magnitude of the momentum at the boundary of filled one-particle states in an isotropic zero-temperature Fermi gas. For the uniform three-dimensional two-component electron model, a short state count givesTo see the factors, impose periodic boundary conditions on a cubic box of side $L$, so $V=L^3$ and each momentum component is $p_i=2\pi \hslash j_i/L$ for an integer $j_i$. Each allowed momentum point occupies a cell of volume $(2\pi \hslash /L{)}^3$ in momentum space. Divide the sphere volume $4\pi p_F^3/3$ by this cell volume and multiply by two for the two spin states. Replacing the discrete count with a volume is the large-box approximation; surface and shell effects are neglected.^[ref]

The factor two is not automatic for every Fermi gas. If only one internal state is populated, its number density $n_{\sigma }$ instead satisfies $p_{F,\sigma }=\hslash (6{\pi }^2{n}_{\sigma }{)}^{1/3}$. For two equally populated states, $n_{\sigma }=n/2$ recovers the displayed result. A Fermi surface in a crystal can have a more complicated shape, so a single spherical Fermi momentum is a property of this isotropic model.

A Fermi point is a point of the zero-temperature Fermi surface in one-dimensional momentum space. For a spin-degenerate band with one occupied interval, its two boundaries are usually $-k_F$ and $+k_F$. Low-energy excitations near them form left- and right-moving branches. A perturbation at wavevector $2k_F$ can mix those branches, as in a Peierls distortion.

The Fermi sea is the zero-temperature many-fermion quantum state with every one-particle state below the Fermi energy occupied and every state above it empty. For a uniform isotropic Fermi gas, the occupied region in momentum space is a solid sphere; its boundary is the Fermi surface. These are momenta, not positions: filling the sphere does not put the particles in a ball in ordinary space. Opposite momentum contributions cancel, so the gas at rest has zero total momentum even though its total kinetic energy is positive.

The filled states are forced upward in energy by the Pauli exclusion principle. This is different from the zero-point energy of one particle confined by a boundary. At fixed number density, the Fermi energy remains positive as the box grows to the thermodynamic limit, while the lowest box confinement energy tends to zero. Low-energy excitations move particles from occupied states near the Fermi surface into nearby empty states; deep states have no nearby available destination. This explains why only a small fraction of electrons contribute to the low-temperature electron heat capacity.^[ref]

The group velocity of an excitation evaluated on the Fermi surface. For an isotropic dispersion relation $ϵ(p)$, $v_F=(dϵ/dp{)}_{p_F}$, where $p_F$ is the Fermi momentum.

For a dispersion relation $E(k)=E_0+t\cos (ka)$, $t<0$, occupied states at zero temperature form the interval $|k|<k_F$ in the first Brillouin zone. Its occupied fraction is $f=2k_F/(2\pi /a)=k_{F}a/\pi$. A spin-degenerate band therefore holds $2f$ electrons per site. If the band is independent of the other two wavevector coordinates, its Fermi surface consists of two planes. For $t>0$ the occupied region is instead near the zone edges, bounded by $|k|=\pi (1-f)/a$.

The edges of an itinerant spin-flip particle-hole continuum follow from both excitation energy and occupation constraints. Consider a zero-temperature three-dimensional parabolic-band model with majority states occupied for $|\mathbf{k}|\le k_F$ and an entirely empty minority band. For exchange splitting $2M$ and ${ϵ}_k={\hslash }^2{k}^2/(2m)$, a majority-to-minority flip carrying wavevector $\mathbf{q}$ hasFor $q>0$, the scalar product is bounded by $\pm k_{F}q$, giving continuum edgesThe extremes occur at opposite points of the occupied Fermi surface along $\mathbf{q}$. At $q=0$ the interval collapses to the single exchange energy $2M$.

The empty-minority assumption requires sufficient splitting; otherwise Pauli exclusion principle blocks some final states and the edges must be recomputed over the allowed region. These are independent particle-hole thresholds, not a calculation of the collective ferromagnetic magnon dispersion. A collective mode below the continuum can therefore have a different low-momentum energy scale.

The sharp Fermi surface gives a nonlocal oscillatory electronic response around an impurity. In a three-dimensional isotropic metal its far-field oscillations have characteristic wavevector $2k_F$ and algebraically decaying amplitude. They are absent from the purely exponential Thomas-Fermi screening potential. The local screening formula is a useful long-wavelength model, not an exact asymptotic expression at every distance.

In a one-carrier isotropic Drude model, $R_H=1/(nq)$ and an inferred carrier density is $n_H=1/(e|R_H|)$. Dividing by the atomic number density gives an effective Hall count per atom. A noninteger result need not mean a literal fractional count of valence electrons: multiple energy bands, anisotropic Fermi surfaces and different charge carrier mobilities can change the Hall coefficient. In a two-carrier isotropic approximation,Its sign identifies the dominant Hall response rather than necessarily the most numerous charge carrier.