Dear gromacs users,
I am conducting an analysis on my system consisting of a DMPC bilayer with
small organic molecules inserting into it. The small molecule consists of a
ring with two substituents opposite each other (like para for a benzene ring).
Below I have tried to sketch the ring of the molecule with the two 'para'
substituents:
2
|
/ \
| |
\ /
|
1
I want to create a vector starting in '1' and having the arrow head in '2' and
then calculate how the angle between this vector and the z-axis varies during
the simulation.
To do this I created an .ndx file consisting of two groups, namely the atom '1'
and the atom '2'. Afterwards I used g_bundle (version 4.0.4 since my
simulations are done in this version):
g_bundle -f xxx.xtc -s xxx.tpr -n xxx.ndx -na 1 -z -ot xxx.xvg
When asked for two groups I first give group 0 corresponding to atom '1' and
the I give group 1 corresponding to atom '2'.
The result is rotated 180 degrees from what I would expect from looking at the
simulation. Am I giving the groups in the wrong order? Will the order used
correspond to a vector starting in '2' and with arrow head in '1'?
By the way, is 1 the correct value for -na in this situation? I find it a bit
difficult to understand that paragraph of the manual.
Thank you,
Sarah
--
gmx-users mailing list [email protected]
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 [email protected].
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
