[Pw_forum] Energy conservation in Verlet-MD?
Norbert.Nemec.List at gmx.de
Thu Jun 24 14:51:36 CEST 2010
I am trying a very simple MD calculation (optical phonon mode of Silicon, similar to the original Car-Parinello paper). I am surprised that the total energy is not constant but oscillates with about 10% of the amplitude of the potential energy.
This oscillation is *independent* of the MD time step. Furthermore it did not change with the k-grid, the energy cutoff, the functional or the pseudopotentials.
Attached is one of the input files that I used.
Does anyone know what might cause this behavior?
GMX DSL: Internet-, Telefon- und Handy-Flat ab 19,99 EUR/mtl.
Bis zu 150 EUR Startguthaben inklusive! http://portal.gmx.net/de/go/dsl
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 705 bytes
Desc: not available
More information about the users