<div dir="ltr">Thank you, dear Dr. Matt W, How about a system including a water droplet in a big box? we must consider COMVEL_TOL? because in this condition(NVT, BOMD) we have deformation and therefore moving the center of mass.<div><br></div><div>BEST</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Oct 15, 2019 at 2:59 AM Matt W <<a href="mailto:mattwa...@gmail.com">mattwa...@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">No. But if there are no external fields then the centre of mass motion should be zero. <div><br></div><div>If it is not (and it won't be) there are some numerical errors. You need to decide what level of inaccuracy is OK for your simulation.<br><br>On Monday, October 14, 2019 at 4:40:13 AM UTC+1, mh H wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);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></blockquote></div></div>

<p></p>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/530c0a66-9ea8-4baa-a978-faffb71ae7c0%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/530c0a66-9ea8-4baa-a978-faffb71ae7c0%40googlegroups.com</a>.<br>
</blockquote></div>