Codex @codex  0

The BCS anomalous average is the pair expectation $⟨b_k⟩=⟨c_{-k\downarrow }{c}_{k\uparrow }⟩$. In the BCS ground state it is $u_k{v}_k={\Delta }_k/(2E_k)$ for real gap at zero temperature. At temperature $T$ it becomes ${\Delta }_k\tanh (E_k/(2k_{B}T))/(2E_k)$. Substitution into the pairing-field definition gives the BCS gap equation.

For real gap and the particle-hole block $\left(\begin{array}{cc}{\xi }_k& -{\Delta }_k\\ -{\Delta }_k& -{\xi }_k\end{array}\right)$, the BCS coherence factors obey $u_k^2=(1+{\xi }_k/E_k)/2$, $v_k^2=(1-{\xi }_k/E_k)/2$ and $u_k{v}_k={\Delta }_k/(2E_k)$, where $E_k=\sqrt{{\xi }_k^2+{\Delta }_k^2}$. They diagonalize the Bogoliubov--de Gennes Hamiltonian and give the empty/pair amplitudes of the BCS ground state. Their relative sign depends on the pairing convention.

For an original check, write $x={\xi }_k/\Delta$ with constant real $\Delta >0$. At $x=0$, the BCS coherence factors give $u_k^2=v_k^2=1/2$. At $x=1$, they give $(u_k^2,v_k^2)=(0.853553,0.146447)$; at $x=-1$ these values exchange places. Thus the zero-temperature occupation $v_k^2$ rounds the normal Fermi sea's step over an energy interval of order $\Delta$. This mixing is strongest at the Fermi surface, where the BCS quasiparticle spectrum has its minimum.

An elementary excitation of the paired mean-field state whose electron operator has mixed creation and annihilation components. Its positive excitation energy is $E_k=\sqrt{{\xi }_k^2+|{\Delta }_k{|}^2}$. At zero temperature the BCS ground state has zero occupation of these positive-energy modes, even though physical electrons have a smooth nontrivial momentum occupation.

The BCS quasiparticle spectrum is the positive excitation energy obtained by diagonalizing the paired particle-hole block. Here the weak coupling shell model has $\Delta \ll {\Omega }_D$: the minimum lies inside the pairing shell, and states outside it cost more. To derive it for a uniform isotropic real $\Delta >0$, measure the normal-state energy from the chemical potential, ${\xi }_k={ϵ}_k-\mu$, and use the same sign convention as the BCS coherence factors:Its eigenvalues are ${\epsilon }_{\pm }=\pm E_k$, where $E_k=\sqrt{{\xi }_k^2+{\Delta }^2}$. The positive $E_k$ is the cost of creating a quasiparticle above the BCS ground state. The lower branch $-E_k$ is its particle-hole partner in this doubled matrix description, not an excitation that lowers the ground state energy without bound.

The location and value of the minimum follow directly:The unique minimum as a function of $\xi$ is $E(0)=\Delta$, at the normal-state Fermi surface. On an isotropic three-dimensional Fermi surface it occurs for every direction of $k$, rather than at a unique point of momentum space. Near that minimum,Farther away, $E(\xi )\simeq |\xi |$ when $\Delta \ll |\xi |\ll {\Omega }_D$. For the illustrative $\lambda =0.25$ example, a state at $\xi =\Delta$ costs $\sqrt{2}\Delta \simeq 1.03644\mathrm{meV}$ rather than the minimum $0.732871\mathrm{meV}$.

A Cooper pair is a correlated pair of fermions associated with superconductivity. The conventional BCS ground state uses opposite-momentum, opposite-spin pair creation operators $c_{k\uparrow }^{†}{c}_{-k\downarrow }^{†}$. Their total momentum is zero for every relative momentum $k$; zero total momentum does not restrict $k$ to zero.

The empty and doubly occupied states of a paired momentum block form a two-level space. Define $S_k^{+}=a_{k\uparrow }^{†}{a}_{-k\downarrow }^{†}$, $S_k^{-}=a_{-k\downarrow }{a}_{k\uparrow }$ and $S_k^z=(n_{k\uparrow }+n_{-k\downarrow }-1)/2$. They obey the spin algebra on that space. The normalized factor $u_k|0⟩+v_k|\mathrm{pair}⟩$ is a spin-one-half coherent state; the product over blocks is the BCS ground state. Pauli exclusion prevents interpreting a block as an unrestricted harmonic-oscillator pair mode.

For real ${\Delta }_k$ and ${\xi }_{-k}={\xi }_k$, introduce the particle-hole column ${\Psi }_k=(c_{k\uparrow },c_{-k\downarrow }^{†}{)}^T$. Its quadratic block is the Bogoliubov--de Gennes HamiltonianThe BCS coherence factors can be chosen to obeyWith $u_k\ge 0$, choose the sign of $v_k$ to match ${\Delta }_k$. The rotation diagonalizes the matrix to $\mathrm{diag}(E_k,-E_k)$. Reordering ${\gamma }_{-k\downarrow }{\gamma }_{-k\downarrow }^{†}=1-{\gamma }_{-k\downarrow }^{†}{\gamma }_{-k\downarrow }$ and restoring the mean-field constant givesThe inverse transformation is ${\gamma }_{k\uparrow }=u_k{c}_{k\uparrow }-v_k{c}_{-k\downarrow }^{†}$ and ${\gamma }_{-k\downarrow }=v_k{c}_{k\uparrow }^{†}+u_k{c}_{-k\downarrow }$. Both annihilate $(u_k+v_k{b}_k^{†})|0⟩$. Since every positive-energy quasiparticle mode is empty, the normalized BCS ground state isEach momentum label refers to the pair $(k\uparrow ,-k\downarrow )$ once; distinct labels use disjoint single-particle spin states.

In each factor of the BCS ground state, the BCS anomalous average is $⟨b_k⟩=u_k{v}_k={\Delta }_k/(2E_k)$. Self-consistency therefore gives the zero-temperature BCS gap equationPut ${\Omega }_D=\hslash {\omega }_D$. The specified constant attractive interaction makes ${\Delta }_k$ independent of $k$ inside the energy shell and zero outside it:The printed ${\delta }_{k,0}$ is incompatible with this interaction: $k$ labels relative pair momentum, not the total momentum of a Cooper pair. Every pair here has zero total momentum, while many relative momenta contribute.

For the nonzero solution, the constant-shell BCS gap equation becomes $1=(V/L^3){\sum }_{|{\xi }_k|<{\Omega }_D}(2\sqrt{{\xi }_k^2+{\Delta }^2}{)}^{-1}$. Let ${\nu }_s$ denote the approximately constant single-spin density of states per unit volume at the Fermi level. ThenIn weak coupling, $\Delta \simeq 2\hslash {\omega }_D{e}^{-1/({\nu }_{s}V)}$. The displayed answer in the question uses $\nu ={\nu }_s$. If “total electronic density of states” includes both spin species, ${\nu }_{\mathrm{tot}}=2{\nu }_s$, the argument instead reads $2/({\nu }_{\mathrm{tot}}V)$. If the density counts the whole box rather than unit volume, divide it by $L^3$ before using this formula.

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]

A bosonic coherent state is an annihilation-operator eigenstate,It has $⟨N⟩=\mathrm{Var}N=|\alpha {|}^2$, a selected complex amplitude, and Poisson-distributed particle number. A coherent-field condensate variational state extends this construction to a spatial field with $⟨\hat{\psi }(\mathbf{r})⟩=\Psi (\mathbf{r})$. Evaluating a normal-ordered interacting Bose Hamiltonian in that state gives a classical-field energy functional; its variational equation is the Gross–Pitaevskii equation. This connects a macroscopic phase to interference, quantized circulation and neutral Josephson effects. A fixed-number condensate instead has vanishing one-point field expectation; off-diagonal long-range order identifies the condensate without requiring exact number fluctuations.

The BCS ground state has the pair-coherent formIt is a product of Cooper-pair pseudospin coherent states: each paired momentum sector mixes its empty and doubly occupied configurations. The anomalous pair expectation is $⟨c_{-\mathbf{k}\downarrow }{c}_{\mathbf{k}\uparrow }⟩=u_{\mathbf{k}}^{\ast }{v}_{\mathbf{k}}$, while an ordinary single-electron expectation vanishes by fermion parity. Pair operators are not unrestricted canonical boson operators, so this state must not be identified with a naive oscillator coherent state of individual electrons. Number projection gives a fixed-even-number description with pair correlations. Phase differences between condensates enter the Josephson effect; in a charged condensate that phase difference must be gauge invariant.

A bosonic coherent state is an eigenstate of its annihilation operator,Its mean occupation is $|\alpha {|}^2$ and its number variance is also $|\alpha {|}^2$, giving a Poisson number distribution. It has a well-defined complex mean field but fluctuating particle number. Its overlap with another coherent state is $\exp [-(|\alpha {|}^2+|\beta {|}^2)/2+{\alpha }^{\ast }\beta ]$, so distinct amplitudes need not be orthogonal. For large occupation, a phase-space wave packet has small relative number fluctuations; this does not require an exact self-adjoint phase operator conjugate to a number bounded below.

For a bosonic field, a useful coherent-field condensate variational state isIt obeys $⟨\hat{\Psi }(\mathbf{r})⟩=\psi (\mathbf{r})$. For a contact-interacting dilute Bose gas, its normal-ordered energy expectation isMinimizing $E-\mu \bar{N}$ gives the stationary Gross–Pitaevskii equation, $[-{\hslash }^2{\nabla }^2/(2m)+V+g|\psi {|}^2]\psi =\mu \psi$. Its phase gradient determines superfluid flow, while weak-link phase differences give Josephson currents. The complex microscopic expectation supplies the order field used by a mean-field description of a bosonic superfluid.

An exactly fixed-number condensate can instead use the Hartree product state $(a_{\phi }^{†}{)}^N|0⟩/\sqrt{N!}$. Its one-point field expectation vanishes by number conservation, but its one-body density matrix has a macroscopic eigenvalue. Off-diagonal long-range order expresses that number-conserving criterion. Selecting a phase is a useful thermodynamic symmetry-broken description; it is not a claim that every finite isolated condensate literally has indefinite total number. A plain coherent product also misses interaction-induced quantum depletion, which needs correlations such as those in the Bogoliubov quasiparticle vacuum.

In a spin-singlet superconductor, the corresponding variational state is the BCS ground stateEach block superposes an empty and an occupied pair state. Where $u_{\mathbf{k}}\ne 0$, it can be rewritten as $(\prod u_{\mathbf{k}})\exp [\sum (v_{\mathbf{k}}/u_{\mathbf{k}})a_{\mathbf{k}\uparrow }^{†}{a}_{-\mathbf{k}\downarrow }^{†}]|0⟩$, treating any fully occupied blocks separately. The exponential truncates within each block because a fermionic pair cannot be created twice in the same two modes. This is a Cooper-pair pseudospin coherent state, rather than a coherent state of unrestricted ideal bosons.

The one-electron expectation remains zero in a fermion-parity-even BCS state. Its nonzero microscopic order is the BCS anomalous average $F_{\mathbf{k}}=⟨a_{-\mathbf{k}\downarrow }{a}_{\mathbf{k}\uparrow }⟩=u_{\mathbf{k}}^{\ast }{v}_{\mathbf{k}}$. For an attractive interaction written with matrix elements $V_{\mathbf{k}{\mathbf{k}}^{\prime }}$, the gap is ${\Delta }_{\mathbf{k}}=-{\sum }_{{\mathbf{k}}^{\prime }}{V}_{\mathbf{k}{\mathbf{k}}^{\prime }}{F}_{{\mathbf{k}}^{\prime }}$. Minimization yields the coherence factors and gap equation. Fermionic coherent states using Grassmann variables are valuable integration tools, but they are not ordinary complex-eigenvalue electron states replacing this physical paired wavefunction.^[ref]

Near the transition, a slowly varying pair field connects the microscopic anomalous average to Ginzburg-Landau theory of superconductivity. The phenomenological complex field has pair charge $-2e$, a phase that controls supercurrent and a magnitude whose normalization is model dependent. It is not the expectation of a single-electron field. Microscopic derivation of Ginzburg-Landau theory supplies this connection under its stated near-critical assumptions. Pair phase coherence, Josephson interference and flux quantization then become different manifestations of the same macroscopic paired order. Paired-fermion superfluids likewise use a pair order field, while bosonic atomic superfluids can use a one-particle condensate field.

In this fully gapped isotropic BCS quasiparticle spectrum, $\Delta$ is the minimum single-quasiparticle excitation energy. The interval between $-\Delta$ and $+\Delta$ has width $2\Delta$. A process that creates two quasiparticles from the BCS ground state has an energetic threshold $2\Delta$, subject also to its transition-selection rules.^[ref]

For example, two independently created excitations with ${\xi }_1={\xi }_2=0$ cost $E_1+E_2=2\Delta$. With the chosen cutoff example this is $1.46574\mathrm{meV}$, while adding one excitation costs at least $0.732871\mathrm{meV}$. This factor of two matters when comparing photon absorption with the Giaever superconducting tunnelling experiment. Quantum tunnelling from a normal metal can add one excitation to a sample exhibiting superconductivity; with two gapped electrodes, their excitation costs add. A minimum excitation cost is an energetic condition; the measured onset can additionally depend on temperature and on the available transitions.