On Thu, Jun 08, 2006, M. Scott Shell wrote:
>
> Given the fact that our minimization was erratic and not entirely
> reproducible and knowing that many of our initial pdb configurations had
> "perfect" dihedral angles, I decided to perturb all initial pdb phi/psi
> angles by a random amount by up to +- 1 degree. As it turns out, this
> fixed our problem and we no longer saw amber hanging during minimization.
Weird...I can't offhand think of any reason why amber9 would be different
than earlier versions in this respect. The basic routines that work with
angles and torsions haven't changed for many years.
> Might this be due to the minimization procedure itself when symmetries
> are perfect along the backbone, or does this sound like an internal
> coding problem perhaps akin to the acos example mentioned above? We
> haven't seen this problem in amber7 or amber8.
I suppose you could see if MD every hangs with the idealized geometry...that
might help decide if it involves minimization per se, or the way in which
forces are being evaluated.
...dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Jun 11 2006 - 06:07:18 PDT
