<div dir="ltr">without considering  COMVEL_TOL in MD section, calculation is wrong? <br><br>On Tuesday, February 9, 2016 at 6:06:40 PM UTC+8, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi Tao,<div><br></div><div>The keyword COMVEL_TOL in MD section will allow you to remove center of mass motion.</div><div><br></div><div>However, if it is drifting very far it probably indicates that your setup is not quite good enough in conserving energy/momentum etc.</div><div><br></div><div>If particles are just diffusing out of the central periodic image, this is not a problem, as you note above. You, or VMD or whatever visualizer you use, can just wrap them back.</div><div><br></div><div>Remember, in CP2K a nice thing is you can run a quick classical calculation to check how it handles COM motion etc with a nicely conserving setup, rather than waiting for a DFT run to finish.</div><div><br></div><div>Matt<br><br>On Tuesday, February 9, 2016 at 9:38:42 AM UTC, tao liu wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi All,<br><br>I have tested the keyword &CENTER_COORDINATES, but this keyword just works for the first step, CP2k will center all the atoms to the default position cell/2.. But in the following steps, the center of atoms will change along with simulation time. Do I understand correctly?  If yes, is there anyway to carry out this &CENTER_COORDINATES keyword in every timestep? Thanks.<br><br>Tao<br></div></blockquote></div></div></blockquote></div>