When the energy curve for bulk system is calculated as a function of the lattice parameter, a sudden change of the number of real space grids is a serious problem which produces an erratic discontinuity on the energy curve. In fact, we see the discontinuity in cases of 200 and 290 (Ryd) in Fig. 6 when the cutoff energy is fixed. The discontinuity occurs at the lattice parameter where the number of grids changes. To avoid the discontinuity on the energy curve, a keyword 'scf.Ngrid' is available.
scf.Ngrid 32 32 32 # n1, n2, and n3 for a-, b-, and c-axes