[gmx-users] grompp error
Nur Syafiqah Abdul Ghani
pqah123 at gmail.com
Tue Oct 23 02:46:09 CEST 2012
Dear All,
Right now i want to do another simulation which is water and co-solvent,HFIP.
I get the small compound file from
http://compbio.biosci.uq.edu.au/atb/download.py?molid=6187#files
I need to used the all atom one and the force field that i used is gromos96.
Right now i'm done packing the HFIP and the water by using packmol and
change all the pdb file into gro file.
Therefore I want to minimize them by using the command like below:
grompp -f em1.mdp -c mixture.gro -p topol.top -o minimize.tpr -v
and the error that coming out is :
Fatal error:
Syntax error - File solvent.itp, line 1
Last line read:
'[ atoms ]'
Invalid order for directive atoms
So, I already checked the problem by gedit the solvent.itp.Here is my itp file :
[ atoms ]
; nr type resnr resid atom cgnr charge mass total_charge
1 CH0 1 hfi C1 1 0.384 12.0110
2 F 1 hfi F1 1 -0.128 18.9984
3 F 1 hfi F2 1 -0.128 18.9984
4 F 1 hfi F3 1 -0.128 18.9984 ; 0.000
.............
9 CH0 1 hfi C3 3 0.384 12.0110
10 F 1 hfi F4 3 -0.128 18.9984
..........
; total charge of the molecule: 0.000
[ bonds ]
; ai aj funct c0 c1
1 2 2 0.1360 7.2300e+06
1 3 2 0.1360 4.7700e+06
...........
9 11 2 0.1360 4.7700e+06
9 12 2 0.1360 7.2300e+06
[ pairs ]
; ai aj funct ; all 1-4 pairs but the ones excluded in GROMOS itp
1 8 1
1 10 1
1 11 1
1 12 1
2 6 1
.......
6 8 1
6 10 1
6 11 1
.......
7 11 1
7 12 1
8 9 1
[ angles ]
; ai aj ak funct angle fc
2 1 3 2 107.60 507.00
2 1 4 2 107.60 507.00
2 1 5 2 111.40 532.00
3 1 4 2 107.60 507.00
......
5 9 10 2 111.40 532.00
5 9 11 2 111.40 532.00
5 9 12 2 109.50 618.00
..........
11 9 12 2 107.60 507.00
[ dihedrals ]
; GROMOS improper dihedrals
; ai aj ak al funct angle fc
[ dihedrals ]
; ai aj ak al funct ph0 cp mult
2 1 5 9 1 0.00 3.77 3
........
[ exclusions ]
; ai aj funct ; GROMOS 1-4 exclusions
and the topology file is :
; Include topology for solvent
#include "solvent.itp"
;
; Include water topology
#include "gromos53a6.ff/spc.itp"
;
#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
; i funct fcx fcy fcz
1 1 1000 1000 1000
#endif
;
[ system ]
; Name
solvent in water box
;
;
[ molecules ]
; Compound #mols
hfi .......
Is thereany typing error that i seem cant see?
I already compare with other problem but seems that I'm lost.
Thanks in advance,
Best Regards,
Nur Syafiqah Abdul Ghani,
Theoretical and Computational Chemistry Laboratory,
Department of Chemistry,
Faculty of Science,
Universiti Putra Malaysia,
43400 Serdang,
Selangor.
013-7188131
alternative email : syafiqahabdulghani at gmail.com
More information about the gromacs.org_gmx-users
mailing list