<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=gb2312">
<META content="MSHTML 6.00.2900.2802" name=GENERATOR></HEAD>
<BODY>
<DIV>David:</DIV>
<DIV> </DIV>
<DIV>> >/ Yinghong wrote:<BR>> />/<BR>> />>/
David:<BR>> />>/<BR>> />>/ >/ Dear Dr. David:<BR>>
/>>/ />/<BR>> />>/ />/ According to the shell water model,
I used this method to another<BR>> />>/ kind of molecule, which is
composed of 6 atoms (e.g. benzene).<BR>> />>/ Initially, I put a dummy
and shell particle (a small mass is given to<BR>> />>/ shell, and doing
a normal dynamcis) in the center of this hexagon, in<BR>> />>/ which
shell particle is connected to dummy through your defined<BR>> />>/
isotropic polarization method.<BR>> />>/ />/<BR>> />>/
/>/ Theoretically, polarization can be looked as a spring-like<BR>>
/>>/ connection with constant Kr = sqr(qS)/4*PHI*Epsilon*Alpha, and
the<BR>> />>/ distance between dummy and shell particle can be decided
by rsd =<BR>> />>/ 4*PHI*Epsilon*Alpha * E0 / qS. Is it right?<BR>>
/>>/ />/<BR>> />>/ />/ Now, in my simulation, I applied an
external electric field along<BR>> />>/ Z direction, and the
interactions (vdws + coulomb) between shell<BR>> />>/ particle and all
the other atoms are exclued. (Of course, here, <BR>> What I /<BR>>
>>/ did is only to make a test instead of a real case). Obviously,
for<BR>> />>/ dummy and shell particles, E0 is currently only referred
to the<BR>> />>/ external field, because local field is
excluded.<BR>> />>/ />/<BR>> />>/ />/
Quantitively, I set alpha = 0.3 nm^3, qS = 3.0e and E0 = 1.5<BR>>
/>>/ V/nm, through "mdrun -debug", alpha and qS can be correctly
output,<BR>> />>/ and the calculated value for Kr = 4168 KJ/mol/nm^2 is
also in the<BR>> />>/ right way. After simulation, I used "g_dist" to
check the distance<BR>> />>/ between dummy and shell particle (rsd)
under such electric field. But<BR>> />>/ the calculted value for rsd is
only 10 percent of the theoretical<BR>> />>/ value although I have
tried for many times.<BR>> />>/ />/<BR>> />>/ />/ So,
Could you tell me some possible errors in my defined model,<BR>> />>/
and why rsd can not approach to the theoretical value? What is the<BR>>
/>>/ principle for GMX to calculate this rsd?<BR>> />>/
/>/<BR>> />>/ /<BR>> />>/ > Isn't the problem nm
vs. Ångström?<BR>> />>/ I am very sure it is not that problem. Upon the
parameters mentioned<BR>> />>/ above, rsd should be ~0.1nm
theoretically, but my calculation gave a<BR>> />>/ value of only
0.01nm. So, any other suggestion?<BR>> />>/ <BR>> />>/
<BR>> />/<BR>> />/<BR>> />/ We have<BR>> />/<BR>>
/>/ F = q E = k r or<BR>> />/ r = q E / k<BR>> />/ r = 0.00108
(eV/kj/mol) nm<BR>> />/ = 0.1 nm<BR>> /<BR>> >
Just realized that I repeated your calculation and got the same
result.<BR>> > How about exclusions? Have you checked the tpr file
for that?<BR>> <BR>> Firstly, thanks for your calculating in person.
Which result did you <BR>> get, 0.1nm or 0.01nm?<BR></DIV>
<DIV>>See above.</DIV>
<DIV> </DIV>
<DIV>> <BR>> In my simulation, I did not define the exclusions in
top file. Instead, <BR>> I defined two energy groups in mdp file: SHELL &
Others. Here, I <BR>> wanna check whether the movement of shell particle is
only related to <BR>> external field in the absence of any other non-bonded
interations <BR>> between shell and other atoms. So, I defined
"energygrp_excl = SHELL <BR>> Others" in mdp files. Is that
right?</DIV>
<DIV><BR>> Maybe, but I'm not sure.</DIV>
<DIV> </DIV>
<DIV><BR>> <BR>> Besides, I checked the tpr file, which seems
ok.<BR></DIV>
<DIV>>Does this mean that all exclusions were there?<BR>>Is there
interaction energy between shell and others in the output edr file?</DIV>
<DIV> </DIV>
<DIV>I can not completely understand the meaning of exclusions in tpr file. But,
I found some information about "energygrp_excl" in tpr file, instead.</DIV>
<DIV>For example, I defined four energy groups in my system: SOL
Benzene SHELL DUMMY. </DIV>
<DIV>Then, the following stuff was what I located in tpr file.</DIV>
<DIV> energygrp_excl[ 0]: 0 0 1
0<BR> energygrp_excl[ 1]: 0 0 1
0<BR> energygrp_excl[ 2]: 1 1 1
0<BR> energygrp_excl[ 3]: 0 0 0 0<BR>Which
should mean I have excluded the interactions between SHELL and the other three
groups. So, I think that is no problem to this point.</DIV>
<DIV>Besides, I read the edr file, which also told me the non-bonded
interactions between SHELL and other groups were all equal to zero.</DIV>
<DIV> </DIV>
<DIV>By the way, on the temperature coupling, I set the temperature as the
follows:</DIV>
<DIV>Tcoupl
=
berendsen<BR>tc-grps
= SOL Benzene SHELL DUMMY
<BR>tau_t
= 0.1
0.1
0.001 0.1
<BR>ref_t
= 300
300
10 300<BR></DIV>
<DIV>As I think, 10K of temperature of SHELL will not influence the magnitude of
shell movement from the viewpoint of hydrodynamics. Because a high temperature
(e.g. 300K) will greatly enhance the kinetics of shell particles, which possibly
made shell particle very active and beyond the control of such a weak
spring with Kr = 4100 KJ/mol/nm^2, as well as finally led to the interrupt
of simulation. So, what is your consideration about this point?</DIV>
<DIV> </DIV>
<DIV>Also, I give a mass of 10.0 a.u. to shell particles. As I observed, so
large a mass will not also influence the magnitude of shell movement, instead,
it is only for decreasing the fluctuation of shell particles. Any doubt?</DIV>
<DIV> </DIV>
<DIV>Thanks.</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV>Xie Yinghong</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV> </DIV></BODY></HTML>