where \(M_{p}=938.3{\rm \;MeV}/c^{2}\) , \(M_{n}=939.6{\rm \;MeV}/c^{2}\)  are the mass of proton and neutron respectively. \({ V}(\vec{ r}_{ k})\)  is the average field, and \(u({ k}_{1}, { k}_{2})\)  is the residual two-body interaction. The first two terms are kinetic energy of protons and neutrons, respectively. The fourth term is the Coulomb potential between protons. To transform the Hamiltonian from first quantization into second quantization, we choose the single particle orbits for protons and neutrons \(\lbrace \psi _{\pi , i}\rbrace _{i=1}^{N_{\pi }}\)  and \(\lbrace \psi _{\nu , j}\rbrace _{j=1}^{N_{\nu }}\)  [1]} respectively, where \(N_{\pi }\)  and \(N_{\nu }\)  are the number of single particle orbits taken in the calculations.
