[gmx-users] editconf not centering, and it is cutting off my molecule
Xu Dong Huang
xudongh at eden.rutgers.edu
Mon Dec 31 22:53:39 CET 2012
Dear advanced gromacs users,
I have created an arbitrary molecule .gro file containing coordinates of my martini-beads. basically, there are 50 beads sitting on each axis of (x,y, -x and -y) forming like a cross. (Each bead with distance 0.43). However, when I try to use editconf on my .gro file to generate a box (editconf -f star.gro -o star_box.gro -c -d 1.0 -bt cubic) , it gives me back a new .gro file with a bad molecule. In VMD viewing, I see that the negative axis part of my molecules are cut off. Closer inspection shows the following for a particular segment of the molecule coordinates:
**original .gro coordinates:
1star O1 170 0.000 -8.170 0.000
1star O1 171 0.000 -8.600 0.000
2star O2 172 0.000 -9.030 0.000
2star O2 173 0.000 -9.460 0.000
2star O2 174 0.000 -9.890 0.000
2star O2 175 0.000 -10.320 0.000
2star O2 176 0.000 -10.750 0.000
2star O2 177 0.000 -11.180 0.000
2star O2 178 0.000 -11.610 0.000
2star O2 179 0.000 -12.040 0.000
2star O2 180 0.000 -12.470 0.000
2star O2 181 0.000 -12.900 0.000
2star O2 182 0.000 -13.330 0.000
2star O2 183 0.000 -13.760 0.000
*new .gro coordinates:
1star O1 169 12.421 4.681 16.695
1star O1 170 12.421 4.251 16.695
1star O1 171 12.421 3.821 16.695
2star O2 172 12.421 3.391 16.695
2star O2 173 12.421 2.961 16.695
2star O2 174 12.421 2.531 16.695
2star O2 175 12.421 22.741 16.695
2star O2 176 12.421 23.171 16.695
2star O2 177 12.421 23.601 16.695
2star O2 178 12.421 24.031 16.695
2star O2 179 12.421 24.461 16.695
2star O2 180 12.421 24.891 16.695
2star O2 181 12.421 25.321 16.695
2star O2 182 12.421 25.751 16.695
2star O2 183 12.421 26.181 16.695
Notice the jump on atom number 174 and 175, why is there a sudden coordinate jump?
All i can see in VMD is my negative axis part of the molecule beads are cut off.
Any input is appreciated. Thank you!
Xu Dong Huang
Chemical & Biochemical Engineering-Undergraduate Researcher
Rutgers School of Engineering
xudongh at eden.rutgers.edu
More information about the gromacs.org_gmx-users
mailing list