[gmx-users] Definition of "polarization"

[Yinghong]

Dear Gmx'ers:

As mentioned in the publication about shell water model, a shell particle is connected to a dummy atom by a spring-like connection with the following relationship:

1. K=sqr(qS)/(4*PHI*Epsilon*alpha), 
2. rsd=(4*PHI*Epsilon*alpha)*E / qS;  (rsd is the distance between shell and dummy particles at any moment, and E is refered to electrical field strength).

Anyone can tell me where is the definition of two equations above, especially in the case of "isotropic polarization", in GMX source code.
Because in my simulation, what I want to polarized is not a water molecule. So I need a definite instruction, any suggestion?

Secondly, how can I read the calculated value of spring constant K by GMX after I defined "alpha" and "qS"? Can "debug" realize it? If yes, how? Although I have used GMX for a long time, but that is the first time for me to encounter such a problem.

Waiting for your response, any instruction will be greatly appreciated.


Xie Yinghong
Hong Kong Univ.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20051214/46672f93/attachment.html>