> Date: Thu, 19 Feb 2009 01:53:47 +1100
> From: mark.abra...@anu.edu.au
> To: gmx-users@gromacs.org
> Subject: Re: [gmx-users] ewald_rtol parameter question
>
> Claus Valka wrote:
> > Hello,
> >
> > I would like to ask you a (simple) clarification question:
> > Reading the manual and equations 4.158 up to 4.160 there is a
> > convergence parameter b (with hellenic letters β).
> > I would like to ask you whether this parameter is the same as the
> > ewald_rtol that I come accross in the manual.
>
> No.
>
> > If it is not, is there any
> > function between those two quantities?
>
> Yes. It's something like erfc(R)/R <= ewald_rtol for R = beta *
> rcoulomb. That is, beta is controlled by the value of the two parameters
> so that the contributions to the direct-space sum from charges outside
> rcoulomb are negligibly small. Check the Darden papers referenced from
> the GROMACS manual for details. GROMACS reports the value of 1/beta in
> the .log file, IIRC.
>
> Mark
ewald_rtol sets the relative magnitude of the direct space eletrostatics part
at the cut-off (as described briefly in the mdp options manual pages), i.e.:
erfc(rcoulomb)/rcoulomb = ewald_rtol
We should add a few lines on this to the manual.
Berk
_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
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-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
