AW: [gmx-users] How do forces from PME get calculated?
dommert at icp.uni-stuttgart.de
Tue Jan 15 23:15:13 CET 2013
There are different possibilities to derive the forces from the PME method.
Usually with SPME, the so called analytical differentiation scheme is
applied, where the gradient of the reciprocal sum is directly calculated.
This is computationally very efficient, because it requires only one iFFT,
and conserves energy. The other method is to perform the differentiation in
reciprocal space, referred to as ik differentiation. This requires 3 iFFTs,
but conserves momentum. The Ewald sum is using the latter method.
There are several articles describing the methods and testing their
accuracy. If you do an error estimate with the GROMACS tool g_pme_error, you
will get a hint to a paper, which discusses the different methods and error
estimates and refers to further literature.
Institut für Computerphysik
> -----Ursprüngliche Nachricht-----
> Von: gmx-users-bounces at gromacs.org [mailto:gmx-users-
> bounces at gromacs.org] Im Auftrag von Schlesier, Thomas
> Gesendet: Dienstag, 15. Januar 2013 21:56
> An: gmx-users at gromacs.org
> Betreff: [gmx-users] How do forces from PME get calculated?
> Maybe a somewhat dumb question:
> How do forces from PME get calculated?
> In the manual (and also some textbooks) i have only found expressions for
> potential for PME (or ewald-summation), but no information how the forces
> One idea i have would be
> dF_i = - [ V(r_i(t)) - V(r_i(t+dt)) ] / [ r_i(t) - r_i(t+dt) ] meaning,
> difference in the potential energy for particle i at two integration steps
> divide this by the distance the particle move times (-1).
> Would this be the right approach?
> gmx-users mailing list gmx-users at gromacs.org
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> * Please don't post (un)subscribe requests to the list. Use the www
> send it to gmx-users-request at gromacs.org.
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
More information about the gromacs.org_gmx-users