advanced28 min studyLesson 534 in Physics

A Fermi surface changes the weak-attraction pairing problem

Derive the logarithmic Cooper binding scale and state exactly what the frozen-sea model establishes.

On this page

What you’ll understand

  • Construct the two-added-particle Cooper problem above a filled Fermi sea.
  • Derive its exponentially small binding scale with a declared DOS convention.
  • Separate a two-particle instability from a superconducting many-body state.

Weak attraction does not bind two particles in every three-dimensional vacuum potential. Near a filled Fermi sea, however, Pauli exclusion changes the available states and a finite density of states produces a logarithmic enhancement. Cooper’s problem isolates this effect by adding two particles while freezing the sea beneath them. The resulting bound state motivates superconducting pairing, but it does not yet establish a self-consistent gap, superfluid stiffness, or the transition temperature of the many-body system.

Freeze the sea and retain a shell of pair states

Take two added fermions in opposite-momentum, opposite-spin states above the Fermi energy. The occupied sea blocks states below it. Let ξ denote the single-particle energy measured upward from the Fermi energy, and retain an attractive matrix element −G only for 0<ξ<Ec. Use a box-normalized convention in which G has energy units and N₀ is the number of allowed pair labels per unit ξ, equal to the single-spin density of states in this simple model. Then g=N₀G is dimensionless.

Assume this density of states is nearly constant in the retained shell, and seek a pair energy below the two-particle threshold by a positive amount ε. The Schrödinger equation makes each pair amplitude proportional to one over 2ξ+ε. Summing those amplitudes gives a scalar consistency condition. The sea has no dynamics here: no other pairs rearrange, no screening is derived, and no phonon mechanism is computed. The attractive interaction and cutoff are inputs to the model, not consequences of its solution.

1=g∫0Ecdξ2ξ+ε=g2ln⁡ ⁣(1+2Ecε),g=N0G>01=g\int_0^{E_c}\frac{d\xi}{2\xi+\varepsilon}=\frac g2\ln\!\left(1+\frac{2E_c}{\varepsilon}\right),\qquad g=N_0G>0
Ec, ξ, ε, and G have the same energy units; N₀ has inverse-energy units in the chosen finite-volume normalization. ε>0 means the two added particles have energy 2EF−ε relative to the unchanged sea. The continuum approximation and constant N₀ are explicit assumptions.↔ Scroll wide equations horizontally.

The logarithm creates an exponentially small scale

As the proposed binding tends to zero, the integral diverges logarithmically because the density of pair states remains nonzero at the threshold. For every positive g the continuum equation therefore has a positive solution. Solving the logarithm gives the exact cutoff-model expression below. At weak coupling it is exponentially small compared with the shell energy. This essential dependence cannot be reproduced by a finite polynomial expansion in g around zero; weak attraction can produce a qualitatively new scale while giving a numerically tiny binding energy.

The constant density of states is central to that conclusion. Near the bottom of a three-dimensional vacuum band the density of states vanishes as the square root of energy, and a binding threshold can remain. In a small finite system, discrete level spacing cuts off the continuum logarithm. A vanishing density of states, a different pairing channel, or a retarded frequency-dependent interaction also changes the calculation. The useful statement is an instability of this Fermi-surface model, not an unrestricted theorem that every arbitrarily weak attractive force binds in every setting.

ε=2Ece2/g−1≃2Ece−2/g,aξ∝12ξ+ε(g≪1)\varepsilon=\frac{2E_c}{e^{2/g}-1}\simeq2E_ce^{-2/g},\qquad a_\xi\propto\frac1{2\xi+\varepsilon}\quad(g\ll1)
aξ is an unnormalized coefficient of a time-reversed pair state. The approximation requires 2/g large. The factor two in the exponent follows from the two added particles and the explicitly single-spin DOS convention.↔ Scroll wide equations horizontally.

A pair instability is the beginning of a many-body theory

In BCS theory, many overlapping pair states reorganize together, and occupied states on both sides of the Fermi surface participate in a self-consistent ground state. Its weak-coupling gap has a different exponential from this two-added-particle binding scale when the same g convention is used. Equating those two quantities would miss that collective rearrangement. Likewise a bound pair alone does not guarantee phase rigidity across a sample; phase fluctuations and dimensionality can separate pair formation from the onset of macroscopic superfluid response.

The waves lab provides a limited picture of coherent amplitude addition. Reset and set both controls to phase difference zero and amplitude ratio one; predict a summed sine wave of amplitude two. Reset and set phase difference 180 degrees and ratio one; predict cancellation. This demonstrates that relative signs matter when amplitudes are combined. It does not represent Pauli blocking, an attractive interaction, a pair wavefunction, or the logarithmic energy integral. The actual Cooper result comes from the allowed-state spectrum and the self-consistency equation, rather than from classical interference alone.

Estimate the small Cooper binding

Use the constant-DOS shell model with g=0.25 and Ec=10 meV. Take exp(8)≈2981. Find the binding energy ε, the pair energy relative to the two-particle threshold, and the ratio ε/Ec.

  1. The exponent is 2/g=8. The exact cutoff-model binding is 2Ec/(exp8−1)=20/2980 meV.

  2. This gives ε≈0.00671 meV. The added pair lies ε below the threshold 2EF, so its energy relative to that threshold is −0.00671 meV.

  3. Divide by Ec=10 meV to obtain ε/Ec≈6.71×10⁻⁴. The small ratio supports the weak-binding shell approximation if the actual level spacing is also much smaller than ε.

The result

The pair binding is approximately 0.00671 meV, giving a threshold-relative energy −0.00671 meV and ε/Ec≈0.000671.

Explore the model

Interactive lab

Interactive lab · waves

Inspect coherent addition of two amplitudes

Reset and set BOTH controls to phase difference 0° and amplitude ratio 1, predicting summed amplitude 2. Reset and set phase difference 180° and ratio 1, predicting zero. These equal-frequency classical waves do not solve a Cooper pair.

Inspect coherent addition of two amplitudesTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forRelative signs affect coherent sums, while the pairing instability specifically requires the fermionic spectrum and attraction.

Check your understanding

What is directly established by the constant-DOS continuum Cooper calculation?

Try it yourself

Keep Ec=10 meV but take g=0.20 and exp(10)≈22026. Estimate the continuum binding. If a small sample has a relevant level spacing 0.01 meV, explain whether that continuum prediction is controlled.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Compare the predicted binding with the spacing after solving the logarithmic model; the existence of a continuum is an independent approximation.

Compare with the solution

The continuum formula gives ε=20/(22026−1)≈0.000908 meV. This is about eleven times smaller than the stated spacing. Replacing the discrete pair spectrum by a constant-DOS integral is therefore not controlled on the predicted binding scale. One must solve the actual discrete-state problem before assigning that binding energy; the continuum logarithm alone is insufficient.

Take it with you

  • Pauli blocking and a finite Fermi-surface DOS produce the Cooper logarithm.
  • The weak binding depends exponentially on the inverse coupling.
  • A two-particle bound state is not yet a theory of superconducting phase rigidity.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.