why are you using such tight threshold? setting such small value and don't care about the convergence (if you can obtain one) makes no sense. Try with a loose conv_thr such as 1D-4 > &electrons > electron_maxstep=100, > conv_thr=7.35D-8, > scf_must_converge=.false. > /