<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=utf-8">
<META content="MSHTML 6.00.2900.2802" name=GENERATOR></HEAD>
<BODY>
<DIV>David:</DIV>
<DIV><BR>><I> Dear Dr. David:<BR></I>><I> <BR></I>><I> According to the
shell water model, I used this method to another kind of molecule, which is
composed of 6 atoms (e.g. benzene). Initially, I put a dummy and shell particle
(a small mass is given to shell, and doing a normal dynamcis) in the center of
this hexagon, in which shell particle is connected to dummy through your defined
isotropic polarization method.<BR></I>><I> <BR></I>><I> Theoretically,
polarization can be looked as a spring-like connection with constant Kr =
sqr(qS)/4*PHI*Epsilon*Alpha, and the distance between dummy and shell particle
can be decided by rsd = 4*PHI*Epsilon*Alpha * E0 / qS. Is it
right?<BR></I>><I> <BR></I>><I> Now, in my simulation, I applied an
external electric field along Z direction, and the interactions (vdws + coulomb)
between shell particle and all the other atoms are exclued. (Of course, here,
What I did is only to make a test instead of a real case). Obviously, for dummy
and shell particles, E0 is currently only referred to the external field,
because local field is excluded. <BR></I>><I> <BR></I>><I>
Quantitively, I set alpha = 0.3 nm^3, qS = 3.0e and E0 = 1.5 V/nm, through
"mdrun -debug", alpha and qS can be correctly output, and the calculated value
for Kr = 4168 KJ/mol/nm^2 is also in the right way. After simulation, I used
"g_dist" to check the distance between dummy and shell particle (rsd) under such
electric field. But the calculted value for rsd is only 10 percent of the
theoretical value although I have tried for many times. <BR></I>><I>
<BR></I>><I> So, Could you tell me some possible errors in my defined model,
and why rsd can not approach to the theoretical value? What is the principle for
GMX to calculate this rsd?<BR></I>><I> <BR></I><BR>> Isn't the problem nm
vs. Ångström?<BR></DIV>
<DIV>I am very sure it is not that problem. Upon the parameters mentioned above,
rsd should be ~0.1nm theoretically, but my calculation gave a value of only
0.01nm. So, any other suggestion?</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV>Thanks</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV>Xie Yinghong</DIV></BODY></HTML>