[gmx-users] Not enough ref_t and tau_t values

Ashutosh Jogalekar ajogale at emory.edu
Thu Mar 16 16:05:56 CET 2006


Thanks Tsjerk
Yes, there was a non-standard amino acid residue and it was not  
recognising it as part of the protein. I changed the total number of  
residues in the aminoacid.dat file (which I had not!) and then  
everything worked ok.
Ashutosh

On Mar 16, 2006, at 3:09 AM, Tsjerk Wassenaar wrote:

> Hi Ashutosh,
>
> If no index file is explicitly given to grompp, it will use default  
> index group (as if you did echo q | make_ndx -f test_b4pr.gro -o  
> index.ndx) Grompp does not know of definitions given in the .top  
> file with regards to the groups. Does your protein contain some non- 
> standard residues? Try to explicitly create an index file and see  
> if the groups "Protein" and "SOL" cover all atoms in your system.  
> ( N(Proteins) + N(SOL) = N(System) ). If not, make adjustments to  
> the index file by merging appropriate groups, i.e. merge non- 
> standard protein groups with the protein if necessary. Note that  
> this basically was the problem given in the last post about "Not  
> enough ref_t and tau_t values".
>
> Cheers,
>
> Tsjerk
>
>
> On 3/16/06, X.Periole <X.Periole at rug.nl> wrote:
> On Wed, 15 Mar 2006 17:21:52 -0600
>   Ashutosh Jogalekar <ajogale at emory.edu> wrote:
> > I know this question has been discussed a number of
> >times on the  list, but almost everytime it has been
> >discussed in reference to the  index file.
> > I am getting this error during grompp for the position
> >restrained MD  even though I am not using any index file.
> >The only groups I have are  solvent and protein and
> >accordingly, there are two columns for these  two in the
> >tau_t and ref_t fields in the pr.mdp file:
> >
> > tau_t=       0.1         0.1
> > tc_grps=   protein  sol
> > ref_t=        300       300
> >
> > my command is:
> >
> > grompp -f pr.mdp -c test_b4pr.gro -p test.top -o
> >test_em.tpr
> >
> > In spite of this, the error is showing up again and
> >again. The top  file also has a row only for protein and
> >solvent.
>
> If all the atoms of your protein or solvent are not
> included
> in the indexes defined by default (protein sol) I think
> grompp find them out of any index and then looks for more
> tau_t and ref_t ... I would check the topology of your
> system. It might be uncorrect.
>
> XAvier
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org .
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
>
>
> -- 
>
> Tsjerk A. Wassenaar, M.Sc.
> Groningen Biomolecular Sciences and Biotechnology Institute (GBB)
> Dept. of Biophysical Chemistry
> University of Groningen
> Nijenborgh 4
> 9747AG Groningen, The Netherlands
> +31 50 363 4336
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php

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


More information about the gromacs.org_gmx-users mailing list