LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Energy drift
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Energy drift

From: Axel Kohlmeyer <akohlmey@...24...>
Date: Wed, 16 May 2018 13:42:47 -0400

On Wed, May 16, 2018 at 11:48 AM, Stefan Paquay <stefanpaquay@...24...> wrote:
This sounds like it might be specific to interactive MD simulations then, but I am not sure. Axel is more of an expert on that I think.

​but i cannot say anything, because the description being given is ambiguous. it is not clear what modifications were made to the example input, and whether the restarting was done correctly.​ in the vast majority of these cases, there is a mistake or misunderstanding in there somewhere, but without a crystal ball, it is impossible to detect and to give advice.


It just seems to me that it is weird that your energy should jump such a specific value upon a time step change.
Energy drift, however, is characterized by a very slow and gradual change of the total energy of the system, which is not at all what you are seeing.

On Tue, May 15, 2018 at 3:39 PM, Chuanxing Zhan <chuanxingzhan@...16...> wrote:
I ran a example in.deca-ala_imd under directory lammps-stable_11Aug2017/examples/USER/misc/imd.
I didn't change anything of example script except inserting restart command. I changed the time-step of system at 11000 step, then the energy jumped. I also checked the energy of the system, and is was continued in the whole simulation. Once I increase time-step, the energy  increase.
I just ran the simulation in one file without restart command. The energy of the system also increased when time-step increase. 

On Tuesday, 15 May 2018, 2:04:45 PM GMT-5, Stefan Paquay <stefanpaquay@...24...> wrote:

This most likely is not drift, because it looks like a sudden jump at a nicely round number of steps (10000). Which example did you run? Are all your pair styles and settings restart-compatible? If you do not restart but run twice as long, is the jump present or not?

On Tue, May 15, 2018 at 2:58 PM, Chuanxing Zhan via lammps-users <> wrote:
Dear lammps users,
Sorry to send a wrong reply email before. I have questions about the energy drift of equilibrating system.
I ran a simulation and the equilibrated at different time-step. However, the energy results of same system show the different energy at different time-step. Then I ran a example script in the lammps example, the energy results of the example system  has the similar problem.

I know the "shadow" Hamiltonian may cause various energy results of same system at different time-step. Did the lammps compute the actual Hamiltonian mechanics of the system? Can anybody give any comments?


------------------------------ ------------------------------ ------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites,!
______________________________ _________________
lammps-users mailing list
lammps-users@...396... lists/listinfo/lammps-users

Check out the vibrant tech community on one of the world's most
engaging tech sites,!
lammps-users mailing list

Dr. Axel Kohlmeyer  akohlmey@...12...24...
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.