In principle, all the asymptotic atomic parameters (multipoles, polarizabilities, dispersion coefficients etc.) are conformation dependent, even for small molecules like water. In our previous work, we usually compute these parameters at equilibrium geometry, and keep them fixed throughout the simulation. In this work, we find that much better accuracy can be achieved if we allow charges and