[gmx-users] .gro => atom number

[Justin A. Lemkul]

Chih-Ying Lin wrote:
> Hi
> 
> .gro file
> 
> MD of 2 waters, t= 0.0
> 6
>     1ABC OW1    1   0.126   1.624   1.679  0.1227 -0.0580  0.0434
>     1ABC HW2    2   0.190   1.661   1.747  0.8085  0.3191 -0.7791
>     1ABC HW3    3   0.177   1.568   1.613 -0.9045 -2.6469  1.3180
>     2ABC OW1    4   1.275   0.053   0.622  0.2519  0.3140 -0.1734
>     2ABC HW2    5   1.337   0.002   0.680 -1.0641 -1.1349  0.0257
>     2ABC HW3    6   1.326   0.120   0.568  1.9427 -0.8216 -0.0244
>    1.82060   1.82060   1.82060
> 
>                         ^
>                         ||
>    atom number (5 positions, integer)
> 
> 
> Since .top only define the force field by the atom number set (1,2,3).
> 

The topology refers to the numbering of a single molecule, so yes, each one has 
three atoms.

> 
> Should this atom number set become =>
> 1                           1
> 2                           2
> 3        =>               3
> 4                           1
> 5                           2
> 6                           3
>

No, atom numbering should be continuous and non-overlapping.  How is your 
topology set up?  If you have something like:

#include "spc.itp"

[ molecules ]
SOL 	2

...then everything should be fine.

-Justin

> 
> 
> Thank you
> Lin
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
> 

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

Justin A. Lemkul
Graduate Research Assistant
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================