Please do not send private emails asking for help. I've said this dozens of
times; it is much more beneficial to the community as a whole to post the
gmx-users list, where discussions are archived for future use. You also have a
better chance of reaching someone who knows how to solve your problem. I know I
respond a lot, but I don't always have an answer, nor do I advertise myself as a
private tutor. Contribution to the list is voluntary, when I've got a few moments.
Your system is blowing up, a fairly common problem. See here for general tips:
http://oldwiki.gromacs.org/index.php/blowing_up
rituraj purohit wrote:
Hello Justin..
I am getting a problem during dynamics(position restrain) of pr.mdp,
at first i given dt = 0.002 and for 200ps (because protein is very big
having 1480 aa = 11058 atoms) shows an error called table size at 1.6
ps. Then i changed it dt = 0.001 for 200ps, but again i had a same
error at ~45 ps.
So finally i giving dt = 0.0001 for 200 ps, now its moving fine..but
steps become 20 lacks :(
1) The correction which i done , is it ok?
You're probably just prolonging the inevitable; the system is not stable.
2) Is error occurred (with dt = 0.002 & 0.001) due to big size of protein?
No; it is due to some physically unrealistic aspect of the model physics
(unresolvable clashes, inappropriate parameters). Size should not be a problem;
people are simulating viral capsids nowadays.
3) To do dynamics for such a big protein is rel event?
That's up to you.
4) Is 200 ps enough for this protein? is their any relation between
protein size and pre-equilibrium period ?
There is no standard rule, really. To justify equilibration, your observables
have to have converged (i.e., temperature, pressure, density, RMSD, etc).
-Justin
Please help me in this..
Regards,
Rituraj
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
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
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
