Hartree-Fock in second quantization and stability of HF solution

The norm of \( |c'\rangle \) is given by (using the intermediate normalization condition \( \langle c' |c\rangle=1 \)) $$ \langle c' | c'\rangle = 1+\sum_{a>F} \sum_{i\le F}|\delta C_{ai}|^2+O(\delta C_{ai}^3). $$ The expectation value for the energy is now given by (using the Hartree-Fock condition) $$ \langle c' |\hat{H} | c'\rangle=\langle c |\hat{H} | c\rangle + \sum_{ab>F} \sum_{ij\le F}\delta C_{ai}^*\delta C_{bj}\langle c |a_{i}^{\dagger}a_{a}\hat{H}a_{b}^{\dagger}a_{j}|c\rangle+ $$ $$ \frac{1}{2!}\sum_{ab>F} \sum_{ij\le F}\delta C_{ai}\delta C_{bj}\langle c |\hat{H}a_{a}^{\dagger}a_{i}a_{b}^{\dagger}a_{j}|c\rangle+\frac{1}{2!}\sum_{ab>F} \sum_{ij\le F}\delta C_{ai}^*\delta C_{bj}^*\langle c|a_{j}^{\dagger}a_{b}a_{i}^{\dagger}a_{a}\hat{H}|c\rangle +\dots $$