<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
Hi,<div><br><div><div>On Jun 19, 2007, at 5:26 PM, Liwei Li wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><div>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.</div><div><br></div></span></blockquote></div><br></div><div>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. </div><div><br class="webkit-block-placeholder"></div><div>We regularly do such NMA on 25,000-atom systems. </div><div><br class="webkit-block-placeholder"></div><div>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.</div><div><br class="webkit-block-placeholder"></div><div>Cheers,</div><div><br class="webkit-block-placeholder"></div><div>Erik</div><div><br class="webkit-block-placeholder"></div><div><br class="webkit-block-placeholder"></div><div><br class="webkit-block-placeholder"></div><div><br class="webkit-block-placeholder"></div><div><br class="webkit-block-placeholder"></div></body></html>