[CP2K-user] [CP2K:13088] Forces in BOMD calculations

[Thomas Kühne]

Dear Dmitrii, 

in principle you are obviously right. However, due to the fact that in a numerical calculation the 
number of particles is finite and typically small, the sum of all forces is fluctuating. If everything 
else is correct, the effect is expected to vanish with increasing system size. You can fix the center 
of mass by setting COMVEL_TOL to 0, but in general this is should not be necessary and if so 
likely points to a source of systematic bias in your calculations. 

Greetings, 
Thomas

> Am 08.04.2020 um 19:59 schrieb Dmitrii Nikolaev <dmitrii.... at gmail.com>:
> 
> Dear collegues,
> 
> 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? 
> 
> Thank you in advance,
> Dmitrii
> 
> -- 
> You received this message because you are subscribed to the Google Groups "cp2k" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to cp... at googlegroups.com <mailto:cp... at googlegroups.com>.
> To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/63f2472a-20a7-4766-a3d3-2cfdacea0f10%40googlegroups.com <https://groups.google.com/d/msgid/cp2k/63f2472a-20a7-4766-a3d3-2cfdacea0f10%40googlegroups.com?utm_medium=email&utm_source=footer>.



==============================
Thomas D. Kühne
Dynamics of Condensed Matter
Chair of Theoretical Chemistry
University of Paderborn
Warburger Str. 100
D-33098 Paderborn
Germany
tdku... at mail.upb.de
+49/(0)5251/60-5726

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200408/440bd71f/attachment.htm>