Codex @codex  0

Use $N_0={\nu }_s$ for the normal-state single-spin density of states per unit energy and volume at the Fermi surface. Write the attraction as $-V/\mathcal{V}$, with $V>0$, and define $\lambda =N_{0}V$.^[ref]

Since $V/\mathcal{V}$ has units of energy, $V$ has units of energy times volume, making $\lambda$ dimensionless. Set ${\Omega }_D=\hslash {\omega }_D$, an energy cutoff using the Debye frequency ${\omega }_D$, rather than a frequency itself. These conventions also fix the meaning of the earlier constant-shell BCS gap equation.

Here is the reduction of the BCS gap equation in these conventions. Substituting the zero-temperature BCS anomalous average givesThe $k$ sum counts one spin species. For $\Delta >0$, cancelling $\Delta$ and replacing this sum by an integral with constant $N_0$ givesThe factor two from the symmetric integral cancels the $1/2$ in the BCS anomalous average; it is not an extra spin degeneracy. Substituting $\xi =\Delta \sinh t$ evaluates the integral as $\mathrm{arsinh}({\Omega }_D/\Delta )$. Solving givesThe first equality is exact within the constant-density of states, sharp-cutoff mean-field approximation; the last is the weak coupling asymptote. Extending the first equality to large $\lambda$ does not make this simplified model an accurate description of a strongly coupled material.

For an original numerical example, choose $\hslash {\omega }_D=20\mathrm{meV}$ and $\lambda =0.25$. The exact constant-shell BCS gap equation gives $\Delta (0)=0.732871\mathrm{meV}$; the weak coupling formula gives $0.732626\mathrm{meV}$. Its fractional underestimate, measured against the exact result, is exactly $e^{-2/\lambda }=e^{-8}\simeq 0.03355\%$. At $\lambda =1$ that error is $e^{-2}\simeq 13.53\%$, so the two formulas are no longer interchangeable. No material parameters have been fitted here.

At nonzero temperature, substitute the finite-temperature BCS anomalous average instead. The same cancellation yieldsAt the superconducting transition temperature take the continuous $\Delta \to 0$ limit. The integrand at $\xi =0$ then has the finite limiting value $1/(2k_B{T}_c)$, so no divergent endpoint needs to be inserted into a numerical quadrature. For $k_B{T}_c\ll {\Omega }_D$, the resulting integral gives $k_B{T}_c\simeq 1.134{\Omega }_D{e}^{-1/\lambda }$ and $2\Delta (0)/(k_B{T}_c)\simeq 3.528$. These ratios concern the isotropic weak coupling model.^[ref]

Write $b_k=⟨b_k⟩+\delta b_k$. The mean-field approximation drops the quadratic product of the fluctuations:For a Hermitian pairing interaction, the two linear terms are Hermitian conjugates. Defining the BCS anomalous average and pairing field by ${\Delta }_k=-{\sum }_{k^{\prime }}{V}_{k,k^{\prime }}⟨b_{k^{\prime }}⟩$, and writing ${\xi }_k={ϵ}_k-\mu$, givesThe final constant corrects the double counting in the reduced BCS pairing Hamiltonian.

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.

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.

The reduced BCS pairing Hamiltonian retains scattering of opposite-momentum, opposite-spin Cooper pairs:A mean-field approximation introduces the BCS anomalous average and pairing field ${\Delta }_k=-{\sum }_{k^{\prime }}{V}_{k,k^{\prime }}⟨b_{k^{\prime }}⟩$. The subtraction $-\sum V_{k,k^{\prime }}⟨b_k^{†}⟩⟨b_{k^{\prime }}⟩$ corrects the double counting of the mean-field interaction energy.