<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><div>Thanks. Just wondering what is a 'distal atom' as you mentioned. </div><div><br></div> I am creating a new subject on g_clustsize. <div><br><br>--- On <b>Wed, 7/7/10, chris.neale@utoronto.ca <i><chris.neale@utoronto.ca></i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: chris.neale@utoronto.ca <chris.neale@utoronto.ca><br>Subject: [gmx-users] problem with trjconv -pbc cluster<br>To: gmx-users@gromacs.org<br>Date: Wednesday, 7 July, 2010, 10:26 PM<br><br><div class="plainMail">I wrote my own program to do this and so I never tried that one. We did find that clustering based on a dingle atom in the chain was a very good idea. So take your distal carbon in the surfactant and make in index group of it and cluster it (with a larger cutoff) with
g_clustsize.<br><br>Obviously if there is a problem with g_clustsize then you should also file a bugzilla.<br><br>If you're going to continue with this thread about g_clustsize, then please pick an appropriate subject line so that others can search it usefully later.<br><br>Chris.<br><br>-- original message --<br><br>Hi Chris, Thanks a lot for your responses. I will surely try them. But, before I go to trjconv -pbc cluster, I had one more problem to sort out. This has to do with quantifying the cluster-size distribution from the trajectory. As I mentioned, if I visualize the trajectory in VMD, I see that starting from a random dispersed state of 50 molecules of surfactant, it aggregates gradually and at a long time, two distinct miceller aggregates are formed.But, If I try to quantify the cluster-size distribution using g_clustsize utilty, it gives me bizarre results: the maxclust.xvg for the final time frame shows that I have 1 single cluster
containing all the 50 molecules. But, Vmd shows 2 distinct aggregate.<br>Here is the command line I used: g_clustsize -f traj.xtc -s topol -mol -cut 0.50 I used -mol option so that I get the cluster size distribution in terms of the molecules.Here, executing the commands gives me maxclust.xvg which has maximum size of clusters( in terms of molecule number) as a function of time frame. It shows at the end, maxclustersize is 50 and it has all the particles in it( as obtained from maxclust.ndx)Besides, at any other time, the maximum size of particles predicted by g_clustsize does not match with what VMD shows. I also tried post-processing the trajectory using trjconv with -pbc whole or with -pbc nojump . But, g_clustsize always returns same result:maxcluster contains all the particles.Then I tried playing with -cut option : I started with default value of 0.35, here the result is worse: it shows maxclust = 1 i.e all the
particles are monomeric . But any value beyond 0.35 gives me maxclust = 50. I tried to<br> cross-check whether visualization is misrepresented or not: for this I tried to calculate all the distances among the com of the surfectants using g_dist tool and tried to manually cluster the surfectants and I found that here I can easily come up with two separate cluster based on the distances. So, I think that the problem lies in g_clustsize not in visualization, as -pbc whole or -pbc nojump does not change the result.So, do you have any suggestions ? Any help in pointing out where I am doing wrong will be helpful.Jagannath<br><br>--- On Wed, 7/7/10, chris.neale at utoronto.ca <chris.neale at utoronto.ca> wrote:<br><br>From: chris.neale at utoronto.ca <chris.neale at utoronto.ca><br>Subject: [gmx-users] problem with trjconv -pbc cluster<br>To: gmx-users at gromacs.org<br>Date: Wednesday, 7 July, 2010, 9:27
PM<br><br><br>--<br>gmx-users mailing list <a ymailto="mailto:gmx-users@gromacs.org" href="/mc/compose?to=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 ymailto="mailto:gmx-users-request@gromacs.org" href="/mc/compose?to=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></div></td></tr></table><br>