Codex @codex  0

The BCS gap equation makes the superconducting gap function self-consistent with the BCS anomalous average that creates it. For a translation-invariant spin-singlet mean-field approximation, define ${\xi }_k={ϵ}_k-\mu$ and $E_k=\sqrt{{\xi }_k^2+|{\Delta }_k{|}^2}$. ThenThe interaction $V_{k{k}^{\prime }}$ and its normalization/cutoff are part of the model. A nonzero solution describes pairing; ${\Delta }_k=0$ is always a solution, so finding any solution alone does not establish which state minimizes the free energy.

For constant attraction $-V$ in a symmetric shell $|\xi |<\Lambda$ with constant single-spin density of states ${\nu }_s$, the zero-temperature equation reduces toAt weak coupling, this gives $\Delta \simeq 2\Lambda e^{-1/({\nu }_{s}V)}$. The exponential explains how a small attraction produces an extremely small pairing scale, matching the Cooper instability. The finite-temperature equation can be linearized as the gap tends to zero to determine the mean-field transition. ^[ref]

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.

The Cooper instability is the logarithmic growth of an attractive BCS coupling at low energy. If ${\nu }_F$ is the density of states and $V<0$, perturbation theory breaks down at a scale proportional to $\Lambda e^{-1/({\nu }_F|V|)}$.

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.

Start from electrons moving in the periodic equilibrium ionic potential. Displace ion $l$ by ${\mathbf{u}}_l$ and expand its contribution as $V(\mathbf{r}-{\mathbf{R}}_l-{\mathbf{u}}_l)\simeq V(\mathbf{r}-{\mathbf{R}}_l)-{\mathbf{u}}_l\cdot \nabla V(\mathbf{r}-{\mathbf{R}}_l)$. Quantizing the harmonic normal modes expresses ${\mathbf{u}}_l$ as a sum of creation and annihilation operators with displacement scale $\sqrt{\hslash /(2{\rho }_{\mathrm{ion}}V{\omega }_{\mathbf{q}})}$ for a single simple branch. Expanding the electronic field in bloch states then gives a band-projected Fröhlich HamiltonianHere $c,c^{†}$ obey canonical anticommutation relations and remove/create an electron with crystal momentum $\mathbf{k}$ and spin $\sigma$; $b,b^{†}$ are boson operators for lattice modes. The matrix element $M_{\mathbf{q}}$ contains the derivative of the ionic potential, phonon polarization, displacement normalization and electronic wavefunction overlaps. With phases chosen appropriately, $M_{-\mathbf{q}}=M_{\mathbf{q}}^{\ast }$ ensures a Hermitian Hamiltonian. The displayed one-band, one-branch notation suppresses additional branch and band labels. “Fröhlich Hamiltonian” here means the linear electron-phonon model used in the question; it is not a requirement to use specifically the polar-optical-phonon coupling of a polar semiconductor. A harmonic lattice, small displacements and a weak enough coupling to organize perturbation theory are the underlying approximations.

The Fröhlich canonical transformation removes real linear phonon operators from the electron-only effective problem to first order. Choose an anti-Hermitian $S$ with $[S,H_0]=-V_{e\text{-}ph}$. For $A_{\mathbf{kq}\sigma }=c_{\mathbf{k}+\mathbf{q},\sigma }^{†}{c}_{\mathbf{k}\sigma }$ and ${\Delta }_{\mathbf{kq}}={ϵ}_{\mathbf{k}+\mathbf{q}}-{ϵ}_{\mathbf{k}}$, one possible generator isIts Hermitian conjugate is its negative after relabelling momenta. The Baker--Campbell--Hausdorff formula givesThe commutator describes virtual emission and absorption of a phonon. normal ordering yields an electron-electron term as well as one-electron energy renormalizations and phonon terms. Retain the former, absorb single-particle shifts into ${ϵ}_{\mathbf{k}}$ and eliminate the explicit phonon sector. For elastic pair scattering the two virtual time orderings have the combined amplitudeUsing the usual factor one-half to count unordered pairs gives the coefficient in the printed Hamiltonian, with $\Delta ={ϵ}_{\mathbf{k}}-{ϵ}_{\mathbf{k}-\mathbf{q}}$. The two outgoing momenta are $\mathbf{k}-\mathbf{q}$ and ${\mathbf{k}}^{\prime }+\mathbf{q}$, so total momentum and the spins are conserved. A fully off-shell effective operator can instead use symmetrized denominators for both external legs; the displayed expression is the usual elastic-scattering reduction. Eliminating nearly resonant real phonons requires more care than this second-order construction.

For $|\Delta |<{\Omega }_{\mathbf{q}}$ the denominator is negative: virtual phonon exchange is attractive in this low-energy channel. Physically, an electron changes the ionic density and a second electron can gain energy by responding to the delayed distortion. This is a retarded interaction, not an unscreened instantaneous attraction at every energy transfer. At large transfer the sign changes, and screened Coulomb repulsion must also be included when deciding the net low-energy pairing strength. The connection between lattice dynamics and superconductivity is supported by the isotope effect, but the mathematical attraction alone does not determine a measured transition temperature.^[ref]

For the pair calculation, make the model restrictions explicit. The filled Fermi sea is inert, and $|0⟩$ denotes that sea with no added particles, not an empty electronic vacuum. Project the interaction onto the previously empty shell $0\le {\xi }_{\mathbf{k}}\le D$, where ${\xi }_{\mathbf{k}}={\hslash }^2{k}^2/(2m)-E_F$ and $D=\hslash {\omega }_D$. Measure the added-pair energy relative to $2E_F$; this is the negative $E$ in the final formula. If absolute one-electron energies are retained in the preceding eigenvalue equation, its eigenvalue must instead be $E_{\mathrm{abs}}=2E_F+E$. A shell-restricted state would not in general be an exact eigenstate of the completely unprojected, all-momentum contact Hamiltonian.

Write $g={\gamma }^2>0$ and $|\mathbf{k}⟩=c_{-\mathbf{k}\downarrow }^{†}{c}_{\mathbf{k}\uparrow }^{†}|0⟩$. Within the zero-total-momentum pair subspace the two opposite-spin terms cancel the interaction's factor one-half, while same-spin s-wave terms cancel by fermionic antisymmetry. Thus the projected Hamiltonian acts asFor $|\Psi ⟩={\sum }_{⟨\mathbf{k}⟩}{\phi }_{\mathbf{k}}|\mathbf{k}⟩$ and $C={\sum }_{⟨\mathbf{k}⟩}{\phi }_{\mathbf{k}}$, matching each basis coefficient givesThis is the stated coefficient formula in the consistent Fermi-level energy convention. Its even dependence on $\mathbf{k}$ produces an even orbital state with a spin singlet. The constant $C$ is fixed up to an overall phase by normalization; for a bound state, $|C{|}^{-2}=g^2{\sum }_{⟨\mathbf{k}⟩}(E-2{\xi }_{\mathbf{k}}{)}^{-2}$. It cannot be set to zero for the nontrivial bound solution.

Put $E=-B$ with $B>0$ and impose the sum defining $C$:Here $\rho (E_F)$ counts orbital states for one spin species, as appropriate to the unrestricted $\mathbf{k}$ sum in this pair basis, and is approximately constant over the thin shell. If a two-spin total density of states is used, it must be divided by two in this integral. State-counting and the volume convention in $g$ and $\rho$ must be consistent. With $\xi =g\rho (E_F)$ as the paper's dimensionless coupling, inversion yieldsThe first approximation is the constant-density-of-states, restricted-shell model; the final exponential additionally assumes weak attraction. The binding energy is nonanalytic at zero coupling, so a finite-order power expansion in $g$ cannot produce it.

Cooper's problem asks whether two electrons added above an inert Fermi sea bind despite Pauli exclusion of lower states. An arbitrarily weak net attraction binds a zero-momentum singlet below the free-pair threshold, because the near-Fermi-surface integral diverges logarithmically as $B\downarrow 0$. Its wavefunction draws on many time-reversed momentum states and becomes dominated by a narrow region near the Fermi surface at weak coupling; the corresponding pair is spatially extended. This establishes a Cooper instability of the unpaired reference state.^[ref] The full superconducting ground state needs a coherent many-electron treatment in which overlapping pairs and the reorganized sea are included. BCS theory supplies that step; the two-electron binding scale alone is neither a calculation of the many-body gap nor a derivation of zero resistance.^[ref]

The tree-level Cooper-pair scattering amplitude is simply $V$. Define the density of states at the Fermi surface, in the normalization of the question, byIn the one-loop Cooper diagram the two internal fermions have opposite momenta. After integrating over internal frequency, the remaining normal-energy integral containsConsequentlyand summing the leading geometric series gives the running BCS theory couplingFor repulsive $V>0$ it flows logarithmically toward zero. For attractive $V<0$, its denominator vanishes at the Cooper instability scaleAt this scale the normal Fermi liquid becomes unstable: opposite-momentum fermions form Cooper pairs, a superconducting or superfluid condensate develops, and a gap opens.

Ginzburg-Landau theory of superconductivity describes the superconducting phase by a complex order parameter and a gauge-invariant free-energy functional. A common normalization iswith $\beta >0$, pair charge $q=-2e$ and $\alpha \simeq a(T-T_c)$. In zero field the uniform minimum has $|\Psi {|}^2=-\alpha /\beta$ below $T_c$. Variation gives the order-parameter equation and supercurrent. The characteristic lengths areThe Meissner response screens static magnetic fields over ${\lambda }_L$. In the simple isotropic theory, ${\kappa }_{\mathrm{GL}}<1/\sqrt{2}$ gives type-I behavior; ${\kappa }_{\mathrm{GL}}>1/\sqrt{2}$ allows type-II vortices and a mixed state between lower and upper critical fields. Flux quantization, vortex structure and currents follow from phase stiffness and electromagnetic coupling. This is a phenomenological expansion most directly controlled near $T_c$; static coefficients alone do not determine every dynamical mode or dissipation rate.

BCS theory supplies a microscopic pairing description. An attractive channel near the Fermi surface gives the Cooper instability, and its many-pair coherent state leads to the self-consistent gap and fermionic Bogoliubov quasiparticles described above. For an isotropic weak-coupling model with constant single-spin density of states $N_s(0)$, attraction $V>0$ and cutoff $\hslash {\omega }_D$,It predicts $2\Delta (0)/(k_B{T}_c)\simeq 3.53$ in that specified model. A full gap suppresses low-temperature quasiparticles exponentially, whereas nodes permit power laws such as the $T^2$ heat capacity derived above. The weak-coupling ratio is not a universal value for all superconductors. Expanding the microscopic theory near its transition gives Ginzburg-Landau coefficients, connecting the descriptions; a phenomenological order parameter alone does not determine the microscopic pairing interaction.

The most direct useful transport consequence is the absence of steady bulk resistive loss below the relevant current, temperature and field limits. Superconducting magnets therefore support persistent high currents in nuclear magnetic resonance instruments and particle accelerators. Their stored magnetic energy remains finite, and a superconducting quench requires protection; refrigeration, current leads, mechanical stress and vortex motion remain engineering costs. Zero dc resistance does not eliminate all ac loss or justify unlimited current.

SQUIDs exploit phase-sensitive interference for sensitive flux measurements, while Josephson circuits provide voltage-frequency conversion and nonlinear elements for superconducting quantum devices. These applications use coherence, not just low resistance. Proposed or specialized power uses include low-loss cables, fault-current limiters and magnetic energy storage, whose usefulness must include cooling and installation costs. Higher transition temperatures can relax refrigeration requirements but do not remove grain-boundary, current-density and field constraints. Neither BCS nor Ginzburg-Landau theory alone supplies a complete device design or guarantees that a particular material is economical.

The electron-phonon interaction lets one electron polarize the ionic lattice and another interact with that delayed deformation. Phonon-mediated electron attraction acts over energy transfers below a characteristic phonon scale. Pairing additionally requires a net attractive interaction in some allowed pairing channel after screened Coulomb repulsion is included. In the conventional weak-coupling situation, retardation separates the low phonon scale from the much larger electronic scale, reducing the effective repulsion. The Cooper instability then makes a Fermi surface unstable to an arbitrarily weak net attraction in the idealized channel; finite temperature cuts off its logarithmic enhancement.

In BCS theory, a momentum-independent pairing amplitude within the cutoff obeys a self-consistency equation containing $\tanh (E/2k_{B}T)/(2E)$, with quasiparticle energy $E=\sqrt{{\xi }^2+|\Delta {|}^2}$. At a continuous transition, take $\Delta \to 0$. A constant effective density of states and coupling combine into $\lambda$, while the band bottom truncates the negative-energy side at $-\sigma$, with $\sigma =\min (\mu ,\theta )$. This gives the supplied linearized gap equation. Its integrand is even, so withit becomes $1/\lambda =[I(\theta ,T_c)+I(\sigma ,T_c)]/2$. If both cutoffs are much larger than $k_B{T}_c$, the given large-cutoff integral gives $I(a,T_c)\simeq \ln (1.13a/(k_B{T}_c))$. Consequently the asymmetric-cutoff BCS transition temperature isFor $\mu >\theta$, this reduces to $1.13\theta e^{-1/\lambda }$; for $0<\mu <\theta$, it is $1.13\sqrt{\theta \mu }{e}^{-1/\lambda }$. The second limit requires $k_B{T}_c\ll \mu$ as well as $k_B{T}_c\ll \theta$. Weak $\lambda$ at fixed nonzero $\mu /\theta$ supplies that limit, but taking $\mu \to 0$ afterward need not preserve it.

For a three-dimensional spin-degenerate parabolic electron band in the degenerate regime,DefineHolding the other model parameters fixed givesThus the log-log slope is $1/3$ below the cutoff-density scale and zero above it. For the specified parameters,This density crossover of a fixed-coupling BCS cutoff model is the requested prescribed model. A real three-dimensional parabolic band also has a changing density of states, so keeping $\lambda$ constant as density varies is an additional imposed assumption, not a consequence of that dispersion.

For noninteracting spinless bosonic pairs, mass $m_B=2m$ and density $n_B=n/2$, homogeneous ideal Bose gas condensation givesUsing the supplied integral $I_B={\int }_0^{\infty }{x}^{1/2}/(e^x-1)dx\simeq 2.3$ gives $C_B=\frac{1}{2}[2/(3I_B){]}^{2/3}\simeq 0.219$. Its slight difference from the more precise $0.218$ is consistent with the rounded integral. Write $a=T_0/(\theta /k_B)=1.13e^{-1/\lambda }\simeq 0.0403$. On the low-density branch,Since $a/C_B<1$, no part of the high-density plateau exceeds $T_B$. The formal crossing condition isThe crossing occurs near $T_c\simeq 7.4\mathrm{K}$ and $T_F\simeq 34\mathrm{K}$.

The apparent low-density excess should not be read as a reliable transition above the ideal condensation scale of preformed pairs. On this extrapolation $k_B{T}_c/\mu =a{u}^{-1/3}$ grows without bound, violating the degenerate-gas and large-lower-cutoff approximations. It reaches one already near $7.5\times 10^{21}{\mathrm{m}}^{-3}$; well before using an arbitrarily small density, the chemical potential, pairing, density of states and pair fluctuations need a self-consistent treatment. The asymptotic formula and ideal-boson curve compare two stipulated models, rather than solving the full BCS-BEC crossover.