I performed a simple simulation with 750 water molecules (lammps 2001). Ewald method and LJ (12-6) potential energy were used for non-bond interactions between the water molecules. However, I found different values of “thermo flag “, which is used to control the frequency of printing the thermodynamic information to screen, had some influence on the results, including the thermodynamic info, and the atom positions at the same timesteps. For example, the “thermo flag” values were chose to be 1 and 50 for two test simulations with the same frequency of dumping the atoms (every 100 timesteps), the differences between the two thermodynamic info printed from the two simulations (at the same timesteps, such as the 2000th steps, etc ) were found. Also there are some differences for the atom positions in the two dump files. The differences become great with the simulation progress.
I check the original procedure “thermo.f” and it seems that this command has no influence on performing simulation (compute force and integrate the equation).
Any comments about this will be highly appreciated.