This method provides the approximate solution for the many-body Schrödinger equation in the ground state. This method was initially introduced by D. R. Hartree at the end of 1920s and further developed by V. Fock in 1930, just after the discovery of the Schrödinger equation . HF method serves as the backbone of the molecular orbital theory which suggests that the motion of each electron can be described with the single-particle function and this particle function does not depend explicitly on instantaneous motions of other electron 57 . The electron-electron interaction is approximated through the Hartree approximation where Nelectron wave-function (Ψ) is considered as the product of single-particle orbitals (𝜓 𝑖 (𝒓 𝑖 𝜎 𝑖 )).