I did a branch of the DFT calculation to find the stable cell parameter (a) .
I didn't get the smooth curve for energy(E) v.s a. There is some point out of concave curve that I had to repeat calculation for them to get the best energy.
Now I doubt about the result and I look for techniques to test the stability of wave function like like something in GAUSSIAN code.