The existing calculations were performed by the method of projector augmented-wave pseudopotentials with density functional theory coded while in the Vienna ab inito simulation package31, 32. For the electrons’ exchange correlation Strength, the Perdew–Burke–Ernzerhof function was used33. The kinetic Electrical power cutoff is chosen as 550?e