<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">

<html>

<head>

<meta http-equiv="content-type" content="text/html; charset=ISO-8859-7">

</head>

<body bgcolor="#ffffff" text="#000066">

<font size="-1"><font face="Helvetica, Arial, sans-serif">Hi

gmx_developers.<br>

<br>

In manual (v4.5) p. 100 ewald_rtol parameter is referred as  "The

ewald_rtol parameter is the relative strength of the electrostatic

interaction at the cut-off"<br>

According to this (and some older posts in lists), beta parameter is

obtained by solving equation <br>

ewald_rtol = erfc( beta * r_c) / r_c     [1]<br>

In gmxlib/ewald_utils.c, function calc_ewaldcoeff()  (versions 3.3.4 -

4.5.4) beta is obtained by solving eq. <br>

ewald_rtol = erfc( beta * r_c)             [2]<br>

(division by r_c omitted ).<br>

These 2 expressions give the same beta only if r_c=1.<br>

I tried to investigate convergence of real/recip. sums as function of

beta. For certain (beta, r_c), I calculate ewald_rtol using [1], put

this ewald_rtol to mdp.<br>

In md.log the reported beta is not the same with the one that was used

to calculate ewald_rtol (not only last few digits deviation), except if

r_c=1.<br>

If I calculate ewald_rtol using [2] and give this value in mdp, mdrun

reports the expected beta.<br>

In almost all systems this doesn't affect results since real/recip.

sums converge for values of beta slightly different (few %) than the

supposed beta,<br>

and typical rcoulomb is close to 1.<br>

This might affect results in some cases that exhibit difficult (if not

impossible) real/recip. sums convergence, that for various reasons

require large cut-offs.<br>

<br>

Is this a bug ?<br>

<br>

DD<br>

<br>

<br>

<br>

<br>

<br>

</font></font>

</body>

</html>