<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Hello all gmx-users:</DIV>
<DIV> I have run a 100ns long MD simulation, the system is HAuCl4 acid solution, which contains one [AuCl4]- cluster and 503 water moleculars, and calculated the rdf by using the g_rdf command as follow:</DIV>
<DIV>g_rdf -f yqqmd.trr -s yqqmd.tpr -n yqq.ndx -o rdf.xvg -cn cn.xvg -rdf atom</DIV>
<DIV>I have defined the groups of AU CL OW H in the yqq.ndx file, and calculated the rdf of Au-Cl, but strangely, all the function value is zero (my box vector is 2.5nm):</DIV>
<DIV># This file was created Thu Nov 26 10:58:07 2009<BR># by the following command:<BR># g_rdf -f yqqmd.trr -s yqqmd.tpr -n yqq.ndx -rdf atom -o rdf.xvg -cn cn.xvg <BR>#<BR># g_rdf is part of G R O M A C S:<BR>#<BR># Gallium Rubidium Oxygen Manganese Argon Carbon Silicon<BR>#<BR>@ title "Cumulative Number RDF"<BR>@ xaxis label "r"<BR>@ yaxis label "number"<BR>@TYPE xy<BR>@ subtitle "AU-CL"<BR> 0 0<BR> 0.002 0<BR> 0.004 0<BR> 0.006 0<BR> 0.008 0<BR>
0.01 0<BR> 0.012 0<BR> 0.014 0<BR> 0.016 0<BR> 0.018 0<BR> 0.02 0<BR> 0.022 0<BR> 0.024 0<BR> 0.026 0<BR> 0.028 0<BR>
0.03 0<BR> 0.032 0<BR> 0.034 0<BR> 0.036 0<BR> 0.038 0<BR> 0.04 0<BR> 0.042 0<BR> 0.044 0<BR> 0.046 0<BR> 0.048 0<BR>
0.05 0<BR> 0.052 0<BR> 0.054 0<BR> 0.056 0<BR> 0.058 0<BR> 0.06 0<BR> 0.062 0<BR> 0.064 0<BR> 0.066 0<BR> 0.068 0<BR>
0.07 0<BR> 0.072 0<BR> 0.074 0<BR> 0.076 0<BR> 0.078 0<BR> 0.08 0<BR> 0.082 0<BR> 0.084 0<BR> 0.086 0<BR> 0.088 0<BR>
0.09 0<BR> 0.092 0<BR> 0.094 0<BR> 0.096 0<BR> 0.098 0<BR> 0.1 0<BR> 0.102 0<BR> 0.104 0<BR> 0.106 0<BR> 0.108 0<BR>
0.11 0<BR> 0.112 0<BR> 0.114 0<BR> 0.116 0<BR> 0.118 0<BR> 0.12 0<BR> 0.122 0<BR> 0.124 0<BR> 0.126 0<BR> 0.128 0<BR>
0.13 0<BR> 0.132 0<BR> 0.134 0<BR> 0.136 0<BR> 0.138 0<BR> 0.14 0<BR> 0.142 0<BR> 0.144 0<BR> 0.146 0<BR> 0.148 0<BR>
0.15 0<BR> 0.152 0<BR> 0.154 0<BR> 0.156 0<BR> 0.158 0<BR> 0.16 0<BR> 0.162 0<BR> 0.164 0<BR> 0.166 0<BR> 0.168 0<BR>
0.17 0<BR> 0.172 0<BR> 0.174 0<BR> 0.176 0<BR> 0.178 0<BR> 0.18 0<BR> 0.182 0<BR> 0.184 0<BR> 0.186 0<BR> 0.188 0<BR>
0.19 0<BR> 0.192 0<BR> 0.194 0<BR> 0.196 0<BR> 0.198 0<BR> 0.2 0<BR> 0.202 0<BR> 0.204 0<BR> 0.206 0<BR> 0.208 0<BR>
0.21 0<BR> 0.212 0<BR> 0.214 0<BR> 0.216 0<BR> 0.218 0<BR> 0.22 0<BR> 0.222 0<BR> 0.224 0<BR> 0.226 0<BR> 0.228 0<BR>
0.23 0<BR> 0.232 0<BR> 0.234 0<BR> 0.236 0<BR> 0.238 0<BR> 0.24 0<BR> 0.242 0<BR> 0.244 0<BR> 0.246 0<BR> 0.248 0<BR>
0.25 0<BR> 0.252 0<BR> 0.254 0<BR> 0.256 0<BR> 0.258 0<BR> 0.26 0<BR> 0.262 0<BR> 0.264 0<BR> 0.266 0<BR> 0.268 0<BR>
0.27 0<BR> 0.272 0<BR> 0.274 0<BR> 0.276 0<BR> 0.278 0<BR> 0.28 0<BR> 0.282 0<BR> 0.284 0<BR> 0.286 0<BR> 0.288 0<BR>
0.29 0<BR> 0.292 0<BR> 0.294 0<BR> 0.296 0<BR> 0.298 0<BR> 0.3 0<BR> 0.302 0<BR> 0.304 0<BR> 0.306 0<BR> 0.308 0<BR>
0.31 0<BR> 0.312 0<BR> 0.314 0<BR> 0.316 0<BR> 0.318 0<BR> 0.32 0<BR> 0.322 0<BR> 0.324 0<BR> 0.326 0<BR> 0.328 0<BR>
0.33 0<BR> 0.332 0<BR> 0.334 0<BR> 0.336 0<BR> 0.338 0<BR> 0.34 0<BR> 0.342 0<BR> 0.344 0<BR> 0.346 0<BR> 0.348 0<BR>
0.35 0<BR> 0.352 0<BR> 0.354 0<BR> 0.356 0<BR> 0.358 0<BR> 0.36 0<BR> 0.362 0<BR> 0.364 0<BR> 0.366 0<BR> 0.368 0<BR>
0.37 0<BR> 0.372 0<BR> 0.374 0<BR> 0.376 0<BR> 0.378 0<BR> 0.38 0<BR> 0.382 0<BR> 0.384 0<BR> 0.386 0<BR> 0.388 0<BR>
0.39 0<BR> 0.392 0<BR> 0.394 0<BR> 0.396 0<BR> 0.398 0<BR> 0.4 0<BR> 0.402 0<BR> 0.404 0<BR> 0.406 0<BR> 0.408 0<BR>
0.41 0<BR> 0.412 0<BR> 0.414 0<BR> 0.416 0<BR> 0.418 0<BR> 0.42 0<BR> 0.422 0<BR> 0.424 0<BR> 0.426 0<BR> 0.428 0<BR>
0.43 0<BR> 0.432 0<BR> 0.434 0<BR> 0.436 0<BR> 0.438 0<BR> 0.44 0<BR> 0.442 0<BR> 0.444 0<BR> 0.446 0<BR> 0.448 0<BR>
0.45 0<BR> 0.452 0<BR> 0.454 0<BR> 0.456 0<BR> 0.458 0<BR> 0.46 0<BR> 0.462 0<BR> 0.464 0<BR> 0.466 0<BR> 0.468 0<BR>
0.47 0<BR> 0.472 0<BR> 0.474 0<BR> 0.476 0<BR> 0.478 0<BR> ......</DIV>
<DIV>then I used the ngmx command to observe the trajectory file of the whole simulation, and indeed see there are four Cl atom around Au atom, (generally speaking, four Cl atoms are at the same plane and form a square, Au is in the center of this square, with a Au-Cl bond length of about 0.228 nm), but I have no idea why the rdf can not reflect the existence of Cl atom around Au?? Can anyone help me on this problem? I will very grateful for your help.</DIV>
<DIV> </DIV>
<DIV>Besides this problem, my Au-OW rdf seems strange also, the second peak, which corresponding to the second hydrated shell, is bigger than the first one. To my konwledge, the second hydrated shell is more "disorded" than the first one so the second rdf peak, generally, is lower and broader than the first one. I guess I had made some mistake in the calculation, but I do not know where is the mistake, does anyone meet this situation ever before?</DIV>
<DIV> </DIV>
<DIV>Any suggestions will be very grateful.</DIV>
<DIV>My English is not good, please forgive my awkward expression.</DIV>
<DIV>Thank you!</DIV>
<DIV> Sincerely: YeQing</DIV></td></tr></table><br>
<hr size=1><a href="http://cn.rd.yahoo.com/mail_cn/tagline/card/*http://card.mail.cn.yahoo.com/"> 好玩贺å¡ç‰ä½ å‘,邮箱贺å¡å…¨æ–°ä¸Šçº¿ï¼</a>