[gmx-users] atoms not part of any of the T-coupling groups
Andy Shelley
robert.shelley at gmail.com
Fri Oct 24 21:35:03 CEST 2008
Garry,
I believe you are trying to not temperature couple the 4 zincs. To do this
put the 4 zincs in one group and include that group in the listing of
tc_grps and then correspondingly assign tau_t. tau_t=0 means no temperature
coupling, but I think you will also need to give a temperature for the
tau_t=0 even though that group does not have temperature coupling. See
section 7.3.13 in the gromacs user manual.
Andy
garry wrote:
>
> I am a novice at Gromacs (although I used Gromos 20+ years ago!)....I
> have a dimer with two zincs per monomer, and I get the following
> complaint when running the grompp with pr.dmp, which I assume is
> referring to the 4 zincs.
>
>
> There are: 6789 OTHER residues
> There are: 262 PROTEIN residues
> There are: 0 DNA residues
> Analysing Protein...
> Analysing Other...
>
> -------------------------------------------------------
> Program grompp, VERSION 3.3.2
> Source code file: readir.c, line: 865
>
> Fatal error:
> 4 atoms are not part of any of the T-Coupling groups
> -------------------------------------------------------
>
> the zincs are in the pdb file as (properly formatted in the pdb file!):
>
> ATOM 2103 SG CYS B 133 -10.574 -5.438 18.502 1.00
> 51.42 S
> ATOM 2104 C CYS B 133 -10.084 -9.287 20.134 1.00
> 50.50 C
> ATOM 2105 O CYS B 133 -9.171 -9.545 20.923 1.00
> 50.69 O
> HETATM 2106 ZN ZN2+C 1 -12.350 1.941 1.066 1.00
> 43.80 ZN
> HETATM 2107 ZN ZN2+C 2 -21.690 2.493 -3.958 1.00
> 38.99 ZN
> HETATM 2109 ZN ZN2+C 3 -11.783 11.509 10.159 1.00
> 28.52 ZN
> HETATM 2108 ZN ZN2+C 4 -17.911 8.618 18.155 1.00
> 32.42 ZN
>
> What am I doing wrong ?
>
What does your .mdp file look like? Have you included the Zn2+ within one
of
your tc-grps?
-Justin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20081024/f164fb9f/attachment.html>
More information about the gromacs.org_gmx-users
mailing list