<div dir="ltr">Thank you very much! <div>My energy is not shifting, and the forces are fluctuating around 20 kcal/(mol*A) -- 0.016 Ht/Bohr. Is it ok? </div><div><br><br>On Thursday, April 9, 2020 at 3:06:17 AM UTC+3, Lucas Lodeiro wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>Hi, </div><div>In the case of BO framework, the forces calculated with Hellmann-Feynman theorem have some instabilities, due to the theorem requires that <d\phi|\phi> = 0 = <\phi|d\phi>. In most cases this terms are very small, but requires a very accurate convergence of energy and the wfc to ensure a completely negligible value.</div><div>Then the real forces are equal <dH> + E<d\phi|\phi>.</div><div>if <d\phi|\phi> is not negligible, you are avoiding a part of the real forces and the summation  of them must not be 0.</div><div><br></div><div>This in the theory. The implementation of this, uses no infinite numbers and some noise are in the numerical evaluation, another thing to <span lang="en"><span title="">no cancel completely. <br></span></span></div><div><span lang="en"><span title=""><br></span></span></div><div><span lang="en"><span title="">The important thing with this issue is if you are using and NVE system, the total energy (kinetic + potential) must be the conserved quantity, not in a perfect numeral way, but the oscillations must be small and without a biased direction, if the average (along the fluctuations) total energy is incrementing or decreasing you have problems in the forces application, either in the force calculation (more accurate scf cycles) or a too big time step integrator.<br></span></span></div></div><br><div class="gmail_quote"><div dir="ltr">El mié., 8 abr. 2020 a las 13:59, Dmitrii Nikolaev (<<a href="javascript:" target="_blank" gdf-obfuscated-mailto="0j08f4VOAAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">dmi...@gmail.com</a><wbr>>) escribió:<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">Dear collegues,<div><br></div><div>I am running the Born-Oppenheimer molecular dynamics, and I notice that the sum of the forces on all atoms of my molecule is non-zero. However, as I understand for the NVE ensemble the total force should be zero? What is the reason for this force? </div><div><br></div><div>Thank you in advance,</div><div>Dmitrii</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="javascript:" target="_blank" gdf-obfuscated-mailto="0j08f4VOAAAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.<wbr>com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/63f2472a-20a7-4766-a3d3-2cfdacea0f10%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msgid/cp2k/63f2472a-20a7-4766-a3d3-2cfdacea0f10%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;" onclick="this.href='https://groups.google.com/d/msgid/cp2k/63f2472a-20a7-4766-a3d3-2cfdacea0f10%40googlegroups.com?utm_medium\x3demail\x26utm_source\x3dfooter';return true;">https://groups.google.com/d/<wbr>msgid/cp2k/63f2472a-20a7-4766-<wbr>a3d3-2cfdacea0f10%<wbr>40googlegroups.com</a>.<br>
</blockquote></div>
</blockquote></div></div>