[gmx-users] grompp fatal error in GROMACS version 4.5.7
Agnivo Gosai
agnivogromacs14 at gmail.com
Wed Oct 22 00:38:00 CEST 2014
Dear GROMACS users
I received the following error while running grompp to prepare the .tpr
file before the initial energy minimization of my system.
*Fatal error:Group DNA_chain_D referenced in the .mdp file was not found in
the index file.Group names must match either [moleculetype] names or custom
index groupnames, in which case you must supply an index file to the '-n'
optionof grompp.*
I used the following few commands in succession before this :-
(1) pdb2gmx -f 1HAO.pdb -o complex.gro -p complex.top
(2) editconf -f complex.gro -o complex_box.gro -bt dodecahedron -d 1.0
(3) genbox -cp complex_box.gro -cs spc216.gro -p complex.top -o
complex_solv.gro
(4) grompp -f em.mdp -c complex_solv.gro -p complex.top -o ions.tpr
(5) genion -s ions.tpr -o complex_ions.gro -p complex.top -pname NA -nname
CL -neutral -conc 0.1
(6) grompp -f em_real.mdp -c complex_ions.gro -p complex.top -o em.tpr
After command (6) error is coming.
In em.mdp file , at step (4) I have , "energygrps = system" and in
em_real.mdp file at step (6) I have ," nenergygrps = DNA_chain_D
Protein_chain_L Protein_chain_H . ( Note : These files I gathered from Dr.
Lemkuls tutorial on protein-ligand system )
Also my complex.top file has the following :-
; Include chain topologies
#include "complex_DNA_chain_D.itp"
#include "complex_Protein_chain_L.itp"
#include "complex_Protein_chain_H.itp"
...........................................................
...........................................................
[ molecules ]
; Compound #mols
DNA_chain_D 1
Protein_chain_L 1
Protein_chain_H 1
SOL 13578
NA 39
CL 28
and ,
"complex_DNA_chain_D.itp" has the following :-
[ moleculetype ]
; Name nrexcl
DNA_chain_D 3
Similar statements are available in the other two .itp files mentioned in
the complex.top file. Now I am unable to understand what the error
statement means by specifying "index file" ??
Please shed some light upon this. I have also consulted the error section
of the GROMACS documentation but could not understand.
Thanks & Regards,
Agnivo Gosai
Grad. Student , Iowa State University
More information about the gromacs.org_gmx-users
mailing list