On 7/28/12 7:34 AM, Shima Arasteh wrote:
Hi all,
My system has got BLOWING UP . I followed the protocol in
http://www.gromacs.org/Documentation/Terminology/Blowing_Up and got as below:
Step 1: If the crash is happening relatively early (within a few steps), set
nstxout (or nstxtcout) to 1, capturing all possible frames. Watch the
resulting trajectory to see which atoms/residues/molecules become unstable
first.
I got this warning:
Warning: 1-4 interaction between 5 and 18 at distance 2.003 which is larger
than the 1-4 table size 2.000 nm
The output shown has nothing to do with the suggested diagnostic step. At what
time step does this occur? Have you done what "Step 1" suggests, and if so,
have you watched the resulting trajectory to see what happens?
Step 2: Simplify the problem:
If you have a new box of solvent, try minimizing and simulating a single
molecule to see if the instability is due to some inherent problem with the
molecule's topology or if instead there are clashes in your starting
configuration.
I got this:
Now there are 24 residues with 360 atoms
Making bonds...
Warning: Long Bond (1-18 = 0.357049 nm)
Again this output is not really relevant to the suggested "Step 2" - this is
from pdb2gmx. What are atoms 1 and 18? It seems that the starting geometry is
very poor. Long bonds indicate missing residues or atoms or just a bad
configuration.
Step 3: Monitor various components of the system's energy using g_energy. If
an intramolecular term is spiking, that may indicate improper bonded
parameters, for example.
The improper-dihedral plot is strange. It "goes up in first picoseconds" and
then goes down slightly.
I think there is something wrong in dihedrals which I defined in .rtp file in
the new residue. atom 18 is the H atom of formyl ( Formyl is as the N-terminus
and is connected to a Valine residue)
Now , what do I need to do? Going to the next step of diagnosing? or no? How
can I fix the problem?
I think you need to back up and do some relevant investigating in steps 1 and 2
first. You haven't really followed any of the advice given thus far so you're
sort of shooting in the dark. Something is clearly wrong either in the topology
or the starting configuration, but it is not apparent at this point which it is.
-Justin
--
========================================
Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/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/Support/Mailing_Lists