Dear all,<br><br>recently I have noticed the same problem with FIST simulations. I took the H2O-32_ewald.inp input, equilibrated it for 1 ns and run a 1ns NVE simulation from the thermalized configuration. I have tried to refine many simulation parameters (e.g. time step, RCUT, Ewald pars) and also (thinking that there was a bad implementation of constraints) I run the simulations for a flexible model. I always get a non-negligible drift (see attachment); the graph shows block averaged energies of 1 ns simulations performed with cp2k and obeli.x (a program of mine). The values are shifted for a better comparison. I have been looking into this, thinking that the problem was in the input file...but since other people are experiencing it with QS, I wonder if there is some bug in the velocity verlet routines. The results are platform independent (I tried both on mac and linux).<br><br>best<br><br>marco<br><br>On Wednesday, September 5, 2012 3:32:26 AM UTC-4, albe...@virgilio.it wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div>Dear all,</div><div><br></div><div>I'm carrying out some tests to reproduce the VDOS of water and I used as</div><div>starting points the input files 'H2O-32.inp' and 'H2O-64.inp' in</div><div>tests/QS/benchmark directory.</div><div><br></div><div>In both cases I ran a NVT simulation at 300K for 15 ps followed by an NVE</div><div>simulation on the equilibrated system.</div><div><br></div><div>I'm puzzled about the results in the energies: in both NVT and NVE  the plot of</div><div>"Cons Qty" show a constant and linear drift while temperature, kinetic and</div><div>potential energy are equilibrated and stable.</div><div><br></div><div>Since "Cons Qty" should be conserved and this is not the case I think that</div><div>there are some problems in the setup of the simulation.</div><div>I did not modify the parameters of the benchmark inputs but I just added the</div><div>thermostat section in NVT and I modified only the printing option.</div><div><br></div><div>Thank you very much for your help and suggestions,</div><div><br></div><div>Alberto</div><div><br></div></blockquote>