[gmx-users] Electrostatic forces
dommert at icp.uni-stuttgart.de
Wed Jun 17 07:43:35 CEST 2009
* Mark Abraham <Mark.Abraham at anu.edu.au> [2009-06-17 15:31:43 +1000]:
> Florian Dommert wrote:
>> * Mark Abraham <Mark.Abraham at anu.edu.au> [2009-06-17 14:14:22 +1000]:
>>> Florian Dommert wrote:
>>>> However I am very confident and in case of success, that there will be
>>>> soon an error estimate for the Ewald Sum available, which will be
>>>> the first
>>>> step to the an implementation a tuning routine for the SPME
>>>> paramters to achieve optimal
>>>> balance between performance and accuracy ;)
>>> I've already implemented a version of mdrun that actually computes
>>> the RMS error in the force components under PME, and am planning to
>>> release it soon.
>> That is very nice to hear, how do you compute the error ? By comparing
>> to an
>> Ewald Sum ?
> Holding beta fixed, I compare force components with those from a
> converged real-space summation and high Fourier grid density &
> interpolation order.
So you have to perform a very costly simulation for every system, when
you gather the reference force ? And which beta do you choose, because
if you take the right choice you can decrease the computational cost
extremely. So theoretically at first you have to find the right beta by
sampling through the corresponding parameter space with a fixed
Interpolation order and grid size. In the optimal range a change of beta
within 0.1 will yield a difference in the error of about 10-1 this trend
continues around +/- 0.5 of the optimal value for beta.
> gmx-users mailing list gmx-users at gromacs.org
> Please search the archive at http://www.gromacs.org/search before posting!
> Please don't post (un)subscribe requests to the list. Use the www
> interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Institute for Computational Physics
Tel: +49 - 711 / 6856-3613
Fax: +49 - 711 / 6856-3658
EMail: dommert at icp.uni-stuttgart.de
!! PGP-ENCODED emails preferred !!
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Size: 194 bytes
Desc: not available
More information about the gromacs.org_gmx-users