[gmx-users] stuck
Adam Fraser
adam.n.fraser at gmail.com
Fri Oct 12 17:51:44 CEST 2007
Hello,
I'm trying to get started working with a system of 2 hexadecane molecules in
water but I'm having some trouble.
I'm getting the following error when I run:
grompp -f simulate.mdp -c 2Hex_solv.gro -p 2Hex.top -o md.tpr
Program grompp, VERSION 3.3.2
Source code file: readir.c, line: 794
Fatal error:
Group protein not found in indexfile.
Maybe you have non-default goups in your mdp file, while not using the '-n'
option of grompp.
In that case use the '-n' option.
So I built an .ndx file with make_ndx and passed it to grompp with -n and it
still doesn't work.
Up until this point I've done the following:
// convert pdb to gro
editconf -f 2Hex.pdb -o 2Hex.gro
// create topology file from .gro
x2top -f 2Hex.gro -o 2Hex.top
// build hexadecane.itp and include it in HexA.top (seems unneccessary)
#include "hexadecane.itp"
// create md input
grompp -f minim.mdp -c 2Hex.gro -p 2Hex.top -o input.tpr
// minimize structure
mdrun -v -s input.tpr -o 2Hex_min.trr -c 2Hex_min.gro
// set up box
editconf -f 2Hex_min.gro -o 2Hex_min_box.gro -d 0.75 -bt cubic
// solvate
genbox -cp 2Hex_min_box.gro -cs spc216.gro -o 2Hex_solv.gro -p HexA.top
#include "spc.itp" (in HexA.top)
// finally attempting to build md input file
grompp -f simulate.mdp -c 2Hex_solv.gro -p 2Hex.top -o md.tpr
Any ideas what I might be doing wrong?
Thanks much,
Adam
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20071012/3c0d591a/attachment.html>
More information about the gromacs.org_gmx-users
mailing list