The way to eliminate the erroneous forces is to use nuclear-position-independent basis functions, to explicitly calculate and then subtract them from the conventionally obtained forces, or to self-consistently optimize the center of localization of the orbitals.
To eliminate erroneous forces, employ nuclear-position-independent basis functions, calculate them explicitly, subtract them from conventionally obtained forces, or self-consistently optimize the center of localization of orbitals.