First, the net force on a neutral system must be zero?
yes. But from the
output, I don't see this. I even tried to run one atom in a big periodic cell. To my surprise, the force on this single atom is not zero?
How far is it from zero? You need also a good grid, i.e. high EnergyCutoff to get it closer to zero.