The BCS theory of Bardeen, Cooper and Schrieffer (1957) is the microscopic theory of conventional superconductivity in metals; its authors received the 1972 Nobel Prize in Physics. Its most quoted quantitative predictions are universal numbers: in the weak-coupling limit the zero-temperature energy gap and the transition temperature are related by Δ(0)=(π/eγ)kBTc≈1.764kBTc, the gap closes as Δ(T)/Δ(0)≈1.741−T/Tc, and the specific heat jumps by (Ces−Cen)/Cen≈1.43 at Tc. These numbers are what experiments compare against, and they follow from a single scalar integral equation, the BCS gap equation. This mission formalizes that mathematical core, following the section Erfolge der BCS-Theorie of the German Wikipedia article BCS-Theorie.
Timeline (as summarized in the source):
1950 — isotope effect observed (Lock et al.; Maxwell; Serin et al.): Tc depends on the ionic mass M.
1956 — Cooper shows that an arbitrarily weak attraction binds two electrons above a filled Fermi sea.
1957 — Bardeen, Cooper and Schrieffer publish Theory of Superconductivity, including the gap equation and its consequences.
Setting
Fix the Boltzmann constantkB>0, the Debye energyℏωD>0 (called hbarOmegaD in Lean) and the dimensionless couplingλ=V0Z(EF)>0 (called coupling), where V0 is the attractive pairing potential and Z(EF) is half the electronic density of states at the Fermi energy. With the quasiparticle energy E(η,Δ)=η2+Δ2 define
The zero-temperature gapΔ(0) is the positive solution of 1/λ=I0(ℏωD,Δ); the critical temperatureTc is the positive solution of 1/λ=I(kB,ℏωD,0,T); and for 0<T<Tc the gapΔ(T) is the positive solution of 1/λ=I(kB,ℏωD,Δ,T). The Fermi–Dirac function is f(E,T)=1/(eE/kBT+1), and 1−2f(E,T)=tanh(E/2kBT).
Formalization targets
Goal: the universal gap ratio
λ→0+limkBTcΔ(0)=eγπ=1.764…
with γ the Euler–Mascheroni constant. The source writes Δ(T=0)=2ℏωDe−1/V0Z(EF)=eγπkBTc; both equalities are weak-coupling asymptotics, so the goal is stated as a limit λ→0+ at fixed kB, ℏωD.
Milestones
I0(ℏωD,Δ)=arsinh(ℏωD/Δ).
Δ(0)=ℏωD/sinh(1/λ), the unique positive solution of the T=0 gap equation.
Δ(0)∼2ℏωDe−1/λ as λ→0+.
The Fermi–Dirac and tanh forms of the finite-temperature gap equation agree.
The Tc equation has exactly one positive solution.
kBTc∼π2eγℏωDe−1/λ (the source's 1,14).
Isotope effect: with ℏωD=C/M, Tc∝M−1/2.
Δ(T) exists uniquely for 0<T<Tc, and the gap equation has no positive solution for T≥Tc.
Near Tc: Δ(T)2/(1−T/Tc) has a limit L(λ) as T→Tc−, and L(λ)/Δ(0)2→8e2γ/(7ζ(3))=1.73672 as λ→0+.
Specific-heat jump: the left derivative D(λ) of Δ(T)2 at Tc exists, and −D(λ)/(32π2kB2Tc)→12/(7ζ(3))=1.426….
Significance
The universal ratios are the standard quantitative test that separates weak-coupling (BCS-like) superconductors from strong-coupling ones; deviations from 1.764 and 1.43 are used to diagnose strong electron–phonon coupling. Mathematically, each statement is an assertion about a parameter-dependent integral equation and its asymptotics as λ→0+, with the Euler–Mascheroni constant and ζ(3) arising from the asymptotics of ∫0Xtanh(x)/xdx and from a Sommerfeld-type expansion near Tc. These results are classical in the physics literature; the mission aims at a machine-checked derivation of the textbook constants from the gap equation, together with reusable lemmas on monotone parameter-dependent integrals.
Difficulty
The T=0 statements reduce to an explicit antiderivative. The finite-temperature statements do not: Tc and Δ(T) are defined only implicitly, so the asymptotics require a uniform expansion ∫0Xtanh(x)/xdx=logX+log(4eγ/π)+o(1). The near-Tc and specific-heat statements additionally concern the behaviour of Δ(T)2 at the boundary point Tc, where Δ itself is not differentiable, and the evaluation of the resulting coefficient in the limit ℏωD/kBTc→∞.
Formalization scope
All quantities are real numbers. Δ(0), Tc and Δ(T) are defined as real suprema (sSup) of the corresponding positive solution sets, so they equal 0 when no positive solution exists; existence and uniqueness are separate milestones. In particular Δ(T)=0 for T≥Tc.
The zero-temperature gap is a separate definition: at T=0 the finite-temperature integrand is not used.
The coupling enters only through λ=V0Z(EF); asymptotic claims are limits λ→0+ at fixed kB,ℏωD>0, with the exact constants 2eγ/π, 8e2γ/(7ζ(3)) and 12/(7ζ(3)) in place of the rounded 1,14, 1,742 and 1,43. ζ(3) is written as the real part of riemannZeta 3.
The specific heat is not modeled thermodynamically; the jump milestone formalizes the source's closed-form expression in terms of −dΔ2/dT at Tc.
Needed infrastructure: monotonicity and continuity of parameter-dependent interval integrals, the integral identity for arsinh, asymptotics of ∫0Xtanhx/xdx, and Mathlib's Euler–Mascheroni constant and ζ(3).
In the weak-coupling BCS model with Debye cutoff ℏωD>0 and Boltzmann constant kB>0, the ratio of the zero-temperature gap Δ(0) (positive solution of λ1=∫0ℏωDη2+Δ2dη) to kBTc (where Tc>0 solves λ1=∫0ℏωDηdηtanh2kBTcη) satisfies
λ→0+limkBTcΔ(0)=eγπ=1.764…,
where λ=V0Z(EF) and γ is the Euler–Mascheroni constant.
Frontier · Open leaf nodes
Recent activity
Frontier · Open leaf nodes
No open leaves. Every sub-goal is proved or awaiting decomposition.