[gmx-users] Problem with using all hydrogen force field

[Sandeep Somani]

Hi all,

I am trying to use the All Hydrogen forcefield (option 1 of pdb2gmx) for
simulation of a protein.  

The protein has a net charge of net charge is -1. The neutralization is
done using genion and the .top file is modified (reduce no of solvent by
one; add one Na molecule; add #include "ions.itp" line in ).

Next for energy minimization grompp is used for generating the .tpr file
but it gives an error: 
"Fatal error: atomtype K not found"

There is no error if Gromacs force field (option 0 of pdb2gmx) is used.

We need to use explicit hydrogens for our simulation but are stuck here
.. any pointers ?

Thnx in advance 

Sandeep

---
Sandeep Somani
Research Associate
BioInformatics Institute
Singapore
http://www.bii.a-star.edu.sg/~ssomani/

"Not everything that counts can be counted, and not everything that can
be counted counts." -- Einstein        


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20040113/46562f8b/attachment.html>