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
The question for today
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.
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.
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.
Work it through
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.
The exponent is 2/g=8. The exact cutoff-model binding is 2Ec/(exp8−1)=20/2980 meV.
This gives ε≈0.00671 meV. The added pair lies ε below the threshold 2EF, so its energy relative to that threshold is −0.00671 meV.
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 pair binding is approximately 0.00671 meV, giving a threshold-relative energy −0.00671 meV and ε/Ec≈0.000671.
The living notebook
Explore the model
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.
- Sum at x = 0
- 0.866
- Resultant amplitude
- 1.732
- Phase difference
- 60°
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.
Pause & reason
Check your understanding
Your turn
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.
Keep the idea
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.