The equilibrium conditions are indeed fulfilled. I'm convinced that there were no convergence problems in this simulation. The oscillations do reduce for larger systems as expected.