 Solution

Group the nonzero modes into independent pairs $p>0$. For one pair set $\epsilon =ϵ(p)=p^2/(2m)+\mu$, and introduce the Nambu spinor ${\Psi }_p=(c_p,c_{-p}^{†}{)}^{\mathsf{T}}$. Normal ordering using $c_{-p}{c}_{-p}^{†}=1-c_{-p}^{†}{c}_{-p}$ givesThis exactly contains both kinetic terms and the pairing terms once; there is no additional factor of one half because only $p>0$ is summed.

Insert ${\Psi }_p=U_p{\mathcal{A}}_p$, where ${\mathcal{A}}_p=({\alpha }_p,{\alpha }_{-p}^{†}{)}^{\mathsf{T}}$ and $U_p$ is the Bogoliubov transformation matrix. Direct matrix multiplication yieldsChoose the diagonalizing branchThen the off-diagonal terms vanish and the upper diagonal entry is $E_p>0$. Extend the angle to negative momenta as in part b; the sine relation just written applies to the positive representative of each pair. The BCS coherence factors can be chosen asHere $u_p$ denotes $\cos {\theta }_p$, and $v_p$ denotes $\sin {\theta }_p$; these expressions are nonsingular because $\mu >0$ and $m>0$. The minus sign of $v_p$ for positive $\Delta$ is fixed by the pairing phase in this Hamiltonian.

A second use of the canonical anticommutation relations givesThe zero-momentum term was never paired and is simply $\mu {\alpha }_0^{†}{\alpha }_0$. HencewithThe constant is the ground-state energy. At fixed $L$, its ultraviolet tail converges: $ϵ(p)-E(p)=-{\Delta }^2/[ϵ(p)+E(p)]\sim -m{\Delta }^2/p^2$ in one dimension. Replacing $E(0)$ by $\sqrt{{\mu }^2+{\Delta }^2}$ would incorrectly pair a spinless mode with itself.