Jennifer Williams wrote:

Hi Justin

Thanks for the suggestions. However, I don't really want to stop the molecule moving into the next image-I start off with an alkyl chain where a large portion of the structure overlaps the pbc.

i.e C1-C2-C3 \\ C4-C5-C6

where || is the unit cell boundary. To see the whole alkyl chain I use the graphics/representations/periodic tab in VMD. The problem is that VMD will not allow me to draw a bond between C3 and C4 meaning I get odd looking movies. VMD prefers to join the C3 to the C4 in the same unit cell rather than the C4 sitting right next to it on the other side of the pbc. This means long bonds stretching the length of my unit cell and missing bonds between atoms sitting next to each other! I have spent a while trying to solve this and have posted to the VMD forum (as it is actually a problem with VMD and not gromacs).


I still think trjconv -pbc mol would be useful to you. It will make the "broken" molecules whole, such that they are not split over the periodic boundaries, which seems like it would eliminate your problem. I do the same with lipid bilayers all the time.

If I could somehow do the same thing to the traj.xtc as I did to confgro.out using genconf -nbox 2 2 2 that would solve my visualisation problem. Do you know of a way to do this? Even a round-about way? I suppose if I dumped each frame in the trajectory as a pdb, used genbox on them to get a 2x2x2 box and then somehow stuck them together into an .xtc this might work (but any suggestions which are less messy are greatly appreciated!)


The only way to do that is indeed the messy way.

-Justin

--
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to