<div dir="ltr"><div><div><div>Dear Wei,<br><br></div>  One should indeed use the unfolded coordinates to calculate the centre of mass, otherwise you get jumps in its position - that is probably indeed why you see any change in the VASP results, as if the (mass-weighted) total force on the atoms is removed (good to know about that behaviour in the VASP code), the centre of mass cannot in principle drift (except for thermostats etc). If there is no option in VASP, one can always do it in post-processing step: If the change in a coordinate of atom I, R_aI(t+dt) - R_aI(t) > L/2 (a = x, y, z), R_aI(t+dt) += L, and then use these to calculate the centre of mass. But, like said, it should remain constant...<br><br></div>    Greetings,<br><br></div>       apsi<br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-01-11 18:05 GMT+01:00 Wei Lai <span dir="ltr"><<a href="mailto:wei.w...@gmail.com" target="_blank">wei.w...@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Ari,<div><br></div><div>As far as saving unfolded coordinates in VASP, I don't know the parameter to control this.  In my opinion, center of mass should be calculated on folded coordinates due to PBC.</div><div><br></div><div>As far as center of mass force correction, thanks for making a good point.  I looked at the source code of VASP 5.4.1, it seems that drift forces (summation of all forces but not exactly forces on the center of mas) are removed at each step.</div><div><br></div><div>Wei<span class=""><br><br>On Wednesday, January 11, 2017 at 11:18:26 AM UTC-5, Ari Paavo Seitsonen wrote:</span><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Dear Wei,<br><br></div><span class="">  Sorry, I did not understand: Did you know correct the coordinates from VASP so that they are not folded back, and still the centre of mass oscillates? Very strange if indeed so... According to my experience there is _always_ a drift in the DFT-based MD simulation (unless one removes the force on the centre of mass at each step), by increasing the accuracy it can made smaller (also in VASP), but basically it never vanishes.<br><br></span></div>    Greetings,<br><br></div>       apsi<br></div><div><div class="h5"><div><br><div class="gmail_quote">2017-01-11 16:27 GMT+01:00 Wei Lai <span dir="ltr"><<a rel="nofollow">wei....@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Ari,<div><br></div><div>Yes, by default VASP folds the coordinates back to the cell and they are fractional between 0 and 1.  In the plot I showed earlier, I simply multiplied the fractional coordinates by lattice parameters.  With VASP, center of mass was oscillating but not drifting.  In cp2k, it was drifting because forces were not calculated accurately enough due to the low cutoff energy.</div><div><br></div><div>Wei <span><br><br>On Wednesday, January 11, 2017 at 10:05:16 AM UTC-5, Ari Paavo Seitsonen wrote:</span><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><span><div dir="ltr"><div><div><div>Dear Wei,<br><br></div>  Just to confirm, a stupid question: When processing the results from VASP did you take into account that by default the code folds the coordinates back into the cell from origin to the lattice vectors, or, the relative coordinates between 0 and 1? That the centre of mass would bounce back and forth is not really a common kind of "drift".<br><br></div>    Greeting from Paris,<br><br></div>       apsi<br></div></span><div><br><div class="gmail_quote"><span>2017-01-10 3:14 GMT+01:00 Wei Lai <span dir="ltr"><<a rel="nofollow">wei....@gmail.com</a>></span>:<br></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><div dir="ltr">Matt,<div><br></div><div>Thanks for the suggestions.  Examining the force noise suggests that higher cutoff values are needed.  Increasing the cutoff values does stabilize the center of mass, although the statistics of fluctuation are different from those from VASP.  But I guess this is probably due to different types of basis sets.</div><span><font color="#888888"><div><br></div></font></span><div><span><font color="#888888">Wei</font></span><div><div><br><br>On Thursday, January 5, 2017 at 10:21:24 AM UTC-5, 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">Dear Wei,<div><br></div><div>this drift indicates that you are not calculating forces accurately enough. You can see forces in the output file by adding <a href="https://manual.cp2k.org/cp2k-4_1-branch/CP2K_INPUT/FORCE_EVAL/PRINT/FORCES.html" rel="nofollow" target="_blank">https://manual.cp2k.org<wbr>/cp2k-4_1-branch/CP2K_INPUT/<wbr>FORCE_EVAL/PRINT/FORCES.html</a></div><div><br></div><div>Possible solutions are:</div><div><br></div><div>(i) use the GAPW method (add METHOD GAPW in the &QS section) - might be slower</div><div>(ii) Increase the CUTOFF - try 600 or 800 ry, at some point GAPW might be better</div><div>(iii) Use smoothing on the XC function</div><div>(iv) just subtract the COM motion to remove the drift (<a href="https://manual.cp2k.org/cp2k-4_1-branch/CP2K_INPUT/MOTION/MD.html#list_COMVEL_TOL" rel="nofollow" target="_blank">https://manual.cp2k.org/cp2k-<wbr>4_1-branch/CP2K_INPUT/MOTION/M<wbr>D.html#list_COMVEL_TOL</a>)<br><br>Matt<br><br>On Wednesday, January 4, 2017 at 1:29:39 AM UTC, Wei Lai wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><p class="MsoNormal" style="line-height:normal"><span style="font-size:10pt;font-family:Arial,sans-serif;background-image:initial;background-position:initial;background-repeat:initial">Dear
cp2k users and developers,</span><span style="font-size:12.0pt;font-family:"Times New Roman",serif"></span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif"> </span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif">I am testing a NVE
simulation with velocity scaling for crystalline CaO (300 K) and having
questions regarding the drift of center of mass.</span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif"> </span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif">In a regular NVE
simulation, I am expecting that (1) the simulation box is fixed and positions
of atoms can be in fractional/scaled or Cartesian units, (2) if atoms moved out
of the simulation box, they can be folded back to the box or left unchanged,
(3) the mass-averaged position centers oscillate around x, y, z direction.  I am showing the center of mass positions
from VASP, which are consistent with the expectation.</span></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif"><br></span></p><p style="text-align:center;clear:both"></p><p style="text-align:center;clear:both"><a href="https://lh3.googleusercontent.com/-tRDdhXVYalE/WGxP2ScXNBI/AAAAAAAAhnk/8bGLtoxRwxINxHoXywkIJFBPOkR4IipMwCLcB/s1600/cp2k-vs-VASP.TIF" style="margin-left:1em;margin-right:1em" rel="nofollow" target="_blank"><img width="400" height="170" border="0"></a></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif"><br></span></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif"><br></span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif"> </span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif">If I do the same for
cp2k (input file enclosed), I am clearly seeing drift of center of mass.  If I analyze the dynamics of the trajectory
(now shown), it appears that both Ca and O atoms are diffusing instead of
vibrating, as obtained from VASP trajectory.</span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif"> </span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif">I am wondering if you
have seen similar behaviors or if I have the wrong settings in the input file.</span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif"> </span></p>

<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal;background-image:initial;background-position:initial;background-repeat:initial"><span style="font-size:10pt;font-family:Arial,sans-serif">Thanks, Wei</span></p></div></blockquote></div></div></blockquote></div></div></div></div></span><div><div><span>

<p></p>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br></span>
To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp2k+...@googlegroups.com</a>.<br>
To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>.<span><br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank">https://groups.google.com/grou<wbr>p/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank">https://groups.google.com/d/op<wbr>tout</a>.<br>
</span></div></div></blockquote></div><br><br clear="all"><span><font color="#888888"><br>-- <br></font></span><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><span><font color="#888888">-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=<wbr>*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=<wbr>-=*=-=*=-=*=-=*=-<br>  Ari Paavo Seitsonen / <a rel="nofollow">Ari.P...@iki.fi</a> / <a href="http://www.iki.fi/~apsi/" rel="nofollow" target="_blank">http://www.iki.fi/~apsi/</a></font></span><span><br>    Ecole Normale Supérieure (ENS), Département de Chimie, Paris<br>    Mobile (F) : +33 789 37 24 25    (CH) : +41 79 71 90 935<br></span></div></div></div></div></div></div>
</div>
</blockquote></div></div><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 rel="nofollow">cp2k+...@googlegroups.com</a>.<br>
To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank">https://groups.google.com/grou<wbr>p/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank">https://groups.google.com/d/op<wbr>tout</a>.<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=<wbr>*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=<wbr>-=*=-=*=-=*=-=*=-<br>  Ari Paavo Seitsonen / <a rel="nofollow">Ari.P...@iki.fi</a> / <a href="http://www.iki.fi/~apsi/" rel="nofollow" target="_blank">http://www.iki.fi/~apsi/</a><br>    Ecole Normale Supérieure (ENS), Département de Chimie, Paris<br>    Mobile (F) : +33 789 37 24 25    (CH) : +41 79 71 90 935<br></div></div></div></div></div></div>
</div>
</div></div></blockquote></div></div><div class="HOEnZb"><div class="h5">

<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:cp2k+uns...@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank">https://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/<wbr>optout</a>.<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-<br>  Ari Paavo Seitsonen / <a href="mailto:Ari.P.S...@iki.fi" target="_blank">Ari.P.S...@iki.fi</a> / <a href="http://www.iki.fi/~apsi/" target="_blank">http://www.iki.fi/~apsi/</a><br>    Ecole Normale Supérieure (ENS), Département de Chimie, Paris<br>    Mobile (F) : +33 789 37 24 25    (CH) : +41 79 71 90 935<br></div></div></div></div></div></div>
</div>