It is the sketch of answer, can be changed further
When one diagonalized the Hamiltonian of BCS theory, the resulting form is
$$H\sim \sum_{k,\sigma}\sqrt{\Delta^2+\xi_k^2}\alpha^{\dagger}_{k,\sigma}\alpha_{k,\sigma}+...,$$
where I have omitted some terms (they are not interesting for the answer), $\sigma$ denotes spin, $\xi_k\equiv\epsilon_k-\mu$ and $\alpha_{k,\sigma}$ is Bogoliuobv quasi-particle.
From this structure it is visible that the minimum energy is simply $\Delta$. This Hamiltonian also tells us that it is difficult to excite the quasi-particle at low temperature: there is the gap $\Delta$ between filled and empty quasi-particle states. It means that ground state (with minimum energy) is rigid.
Next, we should remind the Landau criterion for superfluidity (for details, check vol. 9, chapter 3, par. 23 of Landau Course in Theoretical Physics), which tells us if $\min(E_k/k)>0$, then the superfluid behavior is possible. The Cooper pairs appear when one introduce new particle creation operators for Hamiltonian diagonalization. So, we have: 1) effective bosons appeared, 2) resulting spectrum obeys Landau criterion for superfluidity. Nice.
Next step is to demonstrate somehow that $\Delta$ appears due to the 2nd order phase transition. There are at leas two possibilities: 1) if you are familiar with field-theoretical formalism and Hubbard transformation, you can immediately perform four fermion interaction in Cooper channel and obtain famous Ginzburg-Landau (GL) effective action, 2) you can also verify (from Bogoliubov transformation) that $\Delta$ depends on temperature by power-law. It tells us that there is a continuous phase transition w.r.t to temperature and we can write down expansion of free energy near $T_c$ in terms of $|\Delta|^2$. It is the mentioned GL-expansion.
Next step is to add external electromagnetic field. To do it, you should couple the this field to $\Delta$ as usual. By minimization of free energy, you can derive the GL-equations. From GL-equations, one can immediately obtain so-called Londons' equation, which demonstrates that magnetic field does not penetrate into the superconductor except the thing layer. Nice.
I have seen a lot of books where Londons' equation is written down only for magnetic field, but to illustrate superconductivity one should use electric field. This derivation gives
$$\frac{\partial\bf{j}}{\partial t}=\frac{1}{\delta^2}{\bf E},$$
where relation ${\bf E}=-\partial_t{\bf A}$ is taken into account and $\delta$ is the penetration length. Then, you shold rewrite this expression in $(\omega,{\bf q})$-space,
$${\bf j}(\omega,{\bf q}=0)=\frac{i}{\omega\delta^2}{\bf E}(\omega,{\bf q}=0).$$
Therefore, the conductivity has pole at $\omega=0$. Nice, let us keep it in mind.
Let me remind that similar expression for normal metal the simplest framework to derive conductivity is Drude model. This model results the following expression for AC-conductivity,
$${\bf j}(t)=\frac{ne^2\tau}{m}\frac{1}{1-i\omega\tau}{\bf E}(t)$$
with $E(t)=Ee^{-i\omega t}$ and $j(t)=je^{-i\omega t}$. The conductivity looks like
$$\text{Re}\,\sigma(\omega)=\frac{ne^2\tau}{m}\frac{1+\omega^2\tau^2},$$
$$\text{Im}\,\sigma(\omega)=\frac{ne^2\tau}{m}\frac{\omega\tau}{1+\omega^2\tau^2}.$$
In case of an ideal conductor, we can just take limit $\tau\rightarrow\infty$, which results
$$\text{Re}\,\sigma(\omega)\xrightarrow{\tau\rightarrow \infty}\frac{ne^2}{m}\delta(\omega),$$
$$\text{Im}\,\sigma(\omega)\xrightarrow{\tau\rightarrow\infty}=\frac{ne^2}{m\omega}.$$
So, we coclude that imaginary part of AC-conductivity of an ideal conductor has the pole at $\omega=0$. Moreover, real part of AC-conductiviy diverging at $\omega=0$. The DC-conductivity can be obtained from AC-conductivity by simply setting $\omega\rightarrow 0$.
Therefore, in order to see diverging conductivity, you should:
- Convince yourself somehow that anomalous average $\Delta$ appears in the Fermi system with attractive interaction
- Verify that there is continuous phase transition and $\Delta$ is the order parameter
- Add external electromagnetic field to consideration and obtain corresponding equations of motion for EM field (GL-equations, Londons' equation)
- Perform comparison between Londons' equation with electric field and Drude model of normal metal
- Consider the case of an ideal conductor and observe that conductivity is divergent.
Last point: it is worth mentioning that case of ideal conductor differs from superconductor. One can demonstrate that in case of ideal conductor we have $\partial{\bf B}/\partial t=0$, whereas in case of superconductor we have ${\bf B}=0$, where ${\bf B}$ is the magnetic field.