[gmx-users] Shift Electrostatic Summation

[David van der Spoel]

Shuangxing Dai wrote:
> Hi, all,
>     I was wondering to modify Shift electrostatic summation to a new 
> one. The origin Shift form in Gromacs is like 1/r-1/Rc (Rc is cut-off), 
> I need to modify to erfc(alpha*r)/r-erfc(alpha*Rc)/Rc, and the original 
> long-range correction of Shift in code shift_util.c is q^2/(2*Rc) ( 
> something like that ), I need to 
> q^2*(erfc(alpha*Rc)/(2*Rc)+alpha/sqrt(pi)). So question is:
> How the force to Shift was calculated? Gromacs use the differentiation 
> to energy, got the force formula, then calculated the force? Or just use 
> the numerical differentiation to the energy and got it? If is the latter 
> one, I do not need to modify the force calculation.
> And also, the alpha for Ewald summation, is there an option or 
> parameters to adjust this value?
> Hi, Mark and Berk, if I want do modify shift electrostatic summation as 
> I want, what other files or place do I need to modify, except 
> shift_util.c and force.c?
> Shuangxing Dai

It would be wise to test everything with a table look up. Much easier, 
and the better performance than analytical potentials.

> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmx-users mailing list    gmx-users at 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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php


-- 
David van der Spoel, Ph.D., Professor of Biology
Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:	+46184714205. Fax: +4618511755.
spoel at xray.bmc.uu.se	spoel at gromacs.org   http://folding.bmc.uu.se