[gmx-users] Residue ‘XXX’ not found in residue topology database - SnCl6
Paolo Costa
paolo.costa85 at gmail.com
Fri Apr 17 22:50:20 CEST 2020
Dear GROMACS user,
I am aware that such error occurs frequently for new users as I am.
I am learning how to add new residue in Gromacs; I tried first following a
tutorial (https://www.svedruziclab.com/tutorials/gromacs/2-methane-in-water/ )
and later by creating a benzene, C6H6 molecule, using amber99.ff force
field. Everything goes fine without any problems.
For the next step I tried to create residue for a molecule, SnCl6, for
which there are no parameters (just for training). Thus, within amber99.ff,
I add in the atomtypes.atp the atom Sn with the corresponding mass (the Cl
(chlorine) was already present). Then, I add in ffnonbonded.itp the LJ
parameters for Sn (again for Cl there are already present). Finally I add
in ffbonded.itp the parameters for bond stretch and angle (I calculated
them by Gaussian and VFFDT software).
I created *my .pdb file* as following:
COMPND stannate
REMARK 1 File created by GaussView 6.0.16
HETATM 1 Sn SnCl6 0 0.000 0.000 0.000
Sn
HETATM 2 Cl1 SnCl6 0 -0.000 0.000 2.502
Cl
HETATM 3 Cl2 SnCl6 0 -0.000 2.502 -0.000
Cl
HETATM 4 Cl3 SnCl6 0 2.502 0.000 0.000
Cl
HETATM 5 Cl4 SnCl6 0 -0.000 -2.502 0.000
Cl
HETATM 6 Cl5 SnCl6 0 -2.502 0.000 -0.000
Cl
HETATM 7 Cl6 SnCl6 0 -0.000 -0.000 -2.502
Cl
END
CONECT 1 2 3 4 5
CONECT 1 6 7
CONECT 2 1
CONECT 3 1
CONECT 4 1
CONECT 5 1
CONECT 6 1
CONECT 7 1
and *the .rtp file* which I add in the amber99.ff folder:
[ SnCl6 ]
[ atoms ]
Sn Sn 1.2468 1
Cl1 Cl -0.5411 2
Cl2 Cl -0.5411 3
Cl3 Cl -0.5411 4
Cl4 Cl -0.5411 5
Cl5 Cl -0.5411 6
Cl6 Cl -0.5411 7
[ bonds ]
Sn Cl1
Sn Cl2
Sn Cl3
Sn Cl4
Sn Cl5
Sn Cl6
However, although to my eyes I did everything OK, I get such error.
*Please, can somebody help me to figure out where I am doing wrong? *
Thanks a lot in advance.
Paolo Costa
--
Paolo Costa, Ph.D.
Postdoctoral Researcher
Department of Chemistry and Biomolecular Sciences
University of Ottawa
10 Marie Curie, Ottawa, ON K1N 6N5, Canada
Room number: DRO 326 (D'Iorio Hall)
More information about the gromacs.org_gmx-users
mailing list