 Solution

Use the reduced attractive reduced BCS pairing Hamiltonian with volume V, ${\xi }_k={\epsilon }_k-\mu$ and attraction $\lambda >0$:A specified pairing shell or an ultraviolet renormalization is part of the model. The trace is grand canonical. Insert fermionic coherent states into its time-sliced trace. The resulting Grassmann variables are antiperiodic, and with inverse-energy imaginary time $0\le \tau \le \beta$ the dimensionless action isThis boundary condition is essential for Fermi-Dirac statistics; replacing it by a bosonic periodic trace changes the answer. A complex Hubbard–Stratonovich transformation uses the local identitywith time-slice factors understood. The pairing field renders the fermion action quadratic. For a uniform static field the independent Nambu spinor block $({\psi }_{k\uparrow },{\bar{\psi }}_{-k\downarrow }{)}^T$ hasIts determinant is $-({\omega }_n^2+E_k^2)$. The phase-independent determinant sign cancels under consistent reference normalization. Integrating the fermions gives the BCS pairing-field effective action, ${\mathcal{S}}_{\mathrm{eff}}={\int }_0^{\beta }V|\Delta {|}^2/\lambda d\tau -\mathrm{Tr}\ln {\mathcal{G}}^{-1}$, with the normal-ordering/reference term retained. The fermionic pairing-block determinant normalization uses one independent pairing block for every k; a fully doubled spin-particle-hole trace instead needs one half. The subtracted grand potential of the reduced BCS model is particularly unambiguous:Performing the Matsubara frequency sum yieldsup to a common constant. At zero gap this reduces exactly to the two-spin free-gas result. The Matsubara derivation of the BCS saddle follows by differentiating with respect to ${\Delta }^{\ast }$ and using ${\beta }^{-1}{\sum }_n({\omega }_n^2+E^2{)}^{-1}=\tanh (\beta E/2)/(2E)$ givesThe second equation fixes the chemical potential at specified number. This stationary action is BCS mean-field theory, not an exact integration of pairing-field fluctuations. With constant single-spin density of states ${\nu }_F$ per volume and a symmetric energy cutoff $\Lambda$, weak coupling gives $\Delta (0)\simeq 2\Lambda e^{-1/(\lambda {\nu }_F)}$ and $k_B{T}_c\simeq (2e^{{\gamma }_E}/\pi )\Lambda e^{-1/(\lambda {\nu }_F)}$. The zero-temperature condensation grand-potential density is $-{\nu }_F{\Delta }^2/2$ to leading weak-coupling order; using a spin-summed density would change these displayed factors. The BCS coherence factors are $u_k^2=(1+{\xi }_k/E_k)/2$ and $v_k^2=(1-{\xi }_k/E_k)/2$, and the associated BCS ground state is the normalized product of empty and occupied pair states. ^[ref]

Allowing a slowly varying complex pairing field generates phase stiffness and, near the transition, the Ginzburg-Landau theory of superconductivity. A neutral paired system has a gapless phase sound mode; in a charged bulk superconductor electromagnetic coupling produces the Anderson-Higgs mechanism in a superconductor and the Meissner effect. These collective fluctuations are distinct from the gapped fermionic excitations. A microscopic dynamical calculation is needed to fix collective-mode dispersion and damping; the static saddle alone does not determine them. ^[ref]