Abstract: A modified variable charge method based on iterative fluctuation charge model is proposed which is suitable for molecular dynamics simulation of a canonical ensemble.It compensates kinetic energy loss in simulation with thermal bath in the canonical ensemble.Molecular dynamics simulations show that this method has no energy leaking.Number of iteration in the modified method is reduced and the computational efficiency is increased.

Key words: molecular-dynamics (MD), variable charge method, canonical ensemble