<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Dear Lukasz:</DIV>
<DIV> Thank you for your reply, I have checked my *.ndx file very carefully and there are the follow words in the file:</DIV>
<DIV>[AU]</DIV>
<DIV>1</DIV>
<DIV>[CL]</DIV>
<DIV>2 3 4 5</DIV>
<DIV>[OW]</DIV>
<DIV>.........</DIV>
<DIV>so I think atom indexes for AU & CL really present there, I have produced the yqq.ndx file by the "make_ndx" command from a *.tpr file, and choiced the option "t" (atom type) when defining the defferent groups, By the way, I built the *.ndx file after my simulation carried out, and use the "grompp" command to rebuilt the md.tpr file to add the index information in it, I think "make_ndx" won't write anything into the *.trr file, is there any problem in my operation? Can anybody help me? Thank you very much.</DIV>
<DIV> Best Wishes!</DIV>
<DIV> Sincerely:YeQing <BR>--- <B>09年11月27日，周五, Lukasz Cwiklik <I><lukasz.cwiklik@uochb.cas.cz></I></B> 写道operation：<BR></DIV>
<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid"><BR>发件人: Lukasz Cwiklik <lukasz.cwiklik@uochb.cas.cz><BR>主题: Re: [gmx-users] An unexpected behavior of rdf calculated by the g_rdf command<BR>收件人: "Discussion list for GROMACS users" <gmx-users@gromacs.org><BR>日期: 2009年11月27日,周五,下午10:06<BR><BR>
<DIV class=plainMail>On Fri, Nov 27, 2009 at 6:40 AM, 青 叶 <<A href="http://cn.mc159.mail.yahoo.com/mc/compose?to=yeqing_121@yahoo.com.cn" ymailto="mailto:yeqing_121@yahoo.com.cn">yeqing_121@yahoo.com.cn</A>> wrote:<BR>><BR>> Hello all gmx-users:<BR>> 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:<BR>> g_rdf -f yqqmd.trr -s yqqmd.tpr -n yqq.ndx -o rdf.xvg -cn cn.xvg -rdf atom<BR>> 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):<BR>> # 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>[...]<BR>> 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.<BR>><BR>> 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?<BR>><BR>> Any suggestions will be very grateful.<BR>> My English is not good, please forgive my awkward expression.<BR>> Thank you!<BR>>
Sincerely: YeQing<BR>> ________________________________<BR>> 好玩贺卡等你发，邮箱贺卡全新上线！<BR><BR><BR>Dear YeQing,<BR>Regarding zeros in rdf results - can you please check how your yqq.ndx<BR>file looks like, I mean, are atom indexes for AU & CL really present<BR>there?<BR>(Btw, it would be more convenient if you do not paste so many zeros<BR>from your rdf results.)<BR><BR>Best,<BR>Lukasz<BR><BR>--<BR>Lukasz Cwiklik<BR><A href="http://cwiklik.wordpress.com/" target=_blank>http://cwiklik.wordpress.com</A><BR>--<BR>gmx-users mailing list <A href="http://cn.mc159.mail.yahoo.com/mc/compose?to=gmx-users@gromacs.org" ymailto="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</A><BR><A href="http://lists.gromacs.org/mailman/listinfo/gmx-users"
target=_blank>http://lists.gromacs.org/mailman/listinfo/gmx-users</A><BR>Please search the archive at <A href="http://www.gromacs.org/search" target=_blank>http://www.gromacs.org/search</A> before posting!<BR>Please don't post (un)subscribe requests to the list. Use the<BR>www interface or send it to <A href="http://cn.mc159.mail.yahoo.com/mc/compose?to=gmx-users-request@gromacs.org" ymailto="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</A>.<BR>Can't post? Read <A href="http://www.gromacs.org/mailing_lists/users.php" target=_blank>http://www.gromacs.org/mailing_lists/users.php</A><BR></DIV></BLOCKQUOTE></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>