Force inconsistencies in VASP 5/6, dependent on NSW
Hi all, attached here are four VASP runs. Two in VASP 6.3.2, two in VASP 5.4.4. The behavior we are observing is briefly summarized as:
- when running MD, we observe very large forces that do not sum to zero (i.e. high force drift)
- taking the exact same INCAR but setting NSW=0 mitigates these large forces, and they sum to zero, but OUTCAR still reports significant force drift (though not the same as for NSW > 0)
- despite the inconsistency in forces, the total energies are nearly identical
Any advice would be appreciated, thanks!