- I would suggest first checking convergence of the total energy with respect to energy cutoff, k-point sampling in a systematic manner for a static configuration, i.e., without doing relaxation; you should aim at convergence of the total energy within 1-10 meV per atom.
- You may also consider using a different pseudopotential, e.g., more modern one such as PseudoDojo.
- Another thing, you could try to enable Preconditioner in Iteration control in the Calculator. That might help improving SCF convergence for slab/surface/2D material PW calculations.