<div>Thanks Mark,</div>
<div>here are the links.</div>
<div><a href="http://www.psc.edu/general/software/packages/gromacs/online/water.html">http://www.psc.edu/general/software/packages/gromacs/online/water.html</a></div>
<div> </div>
<div><a href="http://www.mail-archive.com/gmx-users@gromacs.org/msg07653.html">http://www.mail-archive.com/gmx-users@gromacs.org/msg07653.html</a></div>
<div> </div>
<div>I tried your advice with 1 water molecule and since it worked I tried it with SOL.ndx, but I got</div>
<div>D[SOL] -0.0333(+/- 0.0370) 1e-5 cm^2/s</div>
<div>It seems not reasonable to me (negative and its error is bigger than itself). Any idea?</div>
<div>Thanks in Advance/Jamie<br><br></div>
<div class="gmail_quote">On Sun, Sep 20, 2009 at 7:02 PM, Mark Abraham <span dir="ltr"><<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div class="im">Jamie Seyed wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hi Mark,<br>Thanks for the puzzling answer. I tried the index file counting the number of water molecule but I got an error.<br>
"Fatal error:<br>The index group does not consist of whole molecules"<br></blockquote><br></div>OK so apparently it wants atom indices for whole molecules. That seems much more plausible to me than David's advice you quoted (please give URLs to such quotes in future, perhaps you took it out of context...)<br>
<br>So try giving it atom indices for a whole water molecule as an experiment. Seems like an obvious try to me :-) Then, scale it up if it seems to work.
<div class="im"><br><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">I also did not understand what you pointed out about "One of these sets of numbers is larger than the other" ...<br>
</blockquote><br></div>One of the set of atom indices in a protein and the set of molecule indices in a solvated protein will have more members than the other. Instead of asking which is right, you might try out a case that will fail under one assumption and work under the other :-) Then perhaps ask a question with some evidence behind it.
<div class="im"><br><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">questions:<br>1) for msd calculation for water molecules, do I need really make index file for all waters,<br>
2) What I made (count the water molecules in a [mymolecules] index file), did not work. How I suppose to do make that?<br>3) In getting started-Water page that I found in google it only say "g_msd -n index" and did not explain which index... how about if I use oxygen.ndx file??<br>
</blockquote><br></div>Provide links, please. We're not going to duplicate your googling.<br><br>Mark<br><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div class="im">I really appreciate your help. Many Thanks in Advance/Jamie<br><br><br><br></div>
<div class="im"> On Sun, Sep 20, 2009 at 3:33 PM, Mark Abraham <<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> <mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>>> wrote:<br>
<br> Jamie Seyed wrote:<br><br> Dear all,<br> I want to calculate the msd using g_msd. I have a macromolecule<br> and some (let say 100) water molecules and I want to use msd for<br> water only. from the manual and list I think I need to make an<br>
index for molecules (water) and use -mol option.<br> From David's answer on the list I found<br> -----------------------------------------------------<br> "If you know which molecules they are then you just write in a<br>
text editor:<br> [ mymolecules ]<br> 3 5 19 23 45 109"<br> ----------------------------------------------------<br> My question is: let say the macromolecule is coming first in the<br> pdb file and then there are water molecules. If I want to make<br>
an index for only waters, according to this advice, should I<br> write 2 3 4... 100 (just remove 1 that goes for macromolecule)<br> I am wondering if there is a confliction between these<br> numbers(refer to molecules) and atom numbers of<br>
macromolecule.... Would you please clarify this for me??<br><br><br> Try it and see. One of these sets of numbers is larger than the<br> other, and if you use a number in that interval you'll either get an<br>
error or success, either of which will answer your question.<br><br> Mark<br> _______________________________________________<br> gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
</div> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>>
<div class="im"><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<br>
posting!<br> Please don't post (un)subscribe requests to the list. Use the www<br> interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.
<div class="im"><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><br><br><br></div>------------------------------------------------------------------------
<div class="im"><br><br>_______________________________________________<br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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 www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">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>
<div>
<div></div>
<div class="h5">_______________________________________________<br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">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 www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">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></div></blockquote></div><br>