Hi,
On Jun 19, 2007, at 5:26 PM, Liwei Li wrote:
I want to carry out normal mode analysis on a large protein (about
10,000 atoms) without using coarse grain approximations. Amber is
out of the question since it can only handle about 5,000 atoms for
NMA. I am wondering if GROMACS can deal with this assuming running
on a 64-bit machine with 8G memory. I truely appreciate it if you
can help.
If you use release 3.3 or later you won't even need lots of memory
since we're using sparse matrix representation and diagonalization as
soon as you have cutoffs, which should always be the case if you're
doing NMA in vacuo.
We regularly do such NMA on 25,000-atom systems.
The best interaction form is usually shifted interactions with a
cutoff of 1.2 to 1.4 nm. First minimize with steepest descent, then
_double_ precision l-bfgs. This is virtually the only case when you
absolutely need double precision. The rlist variable should be
~0.2-0.3nm larger to conserve energy with l-bfgs, since that helps
convergence.
Cheers,
Erik
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
