----- Original Message -----
From: Oleksandr <nablaobl...@yahoo.com>
Date: Monday, May 31, 2010 18:48
Subject: [gmx-users] water clusters MD
To: gmx-users@gromacs.org

> Hi,
> I'm new user of Gromacs. There are plenty examples how to do 
> solvation study of various proteins. But I'd like to "solvate" 
> highly ordered heavy water clusters in box of light water. At 
> the first step when I execute  "pdb2gmx -f watercluster.pdb" I 
> get an error no matter which model is chosen:
> -------------------------------------------------------
> 
> Fatal error:
> Atom H in residue HOH 4 not found in rtp entry with 3 atoms
>              while sorting atoms. Maybe different protonation state.
>              Remove this hydrogen or choose a different 
> protonation state.
>              Option -ignh will ignore all hydrogens in the input.
> -------------------------------------------------------
> 
> Can anybody help me how to solve this problem?

You're using a tool for a different job. pdb2gmx builds a topology file, mostly 
by constructing a [moleculetype] from a polymer of repeating units. Your 
moleculetypes are trivial, and at least one is already present in an .itp file 
for use with #include.

editconf and genbox make a box and fill it with generic solvent. It will be 
simplest to either write your .top by hand, or adapt an existing .top, 
depending what you mean by "heavy water". Either way, you'll need some fluency 
with GROMACS workflows and file types, so do all the general tutorial material 
you can find. If you'd done so, you might have realised that pdb2gmx is not 
really the tool for the job.

Mark
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to