On Thu, Jan 01, 2009, Paul Ledbetter wrote:
>
> I have been experimenting with changing force field parameters for a
> model of a relatively small molecule (22 atoms). For some choices of
> parameters, I see very strange results during minimization.
> Specifically, the energy of successive geometries during the
> minimization will sometimes jump to a higher energy than the starting
> geometry. This happens with greater probability as the number of
> cycles increases. I have set the non-bonded cutoff to larger than the
> system, so I am unable to think of a good reason for this to happen.
> I have included an mdout from one particular run. Any advice would be
> appreciated.
This is part of the nature of the conjugate gradient routine used. You
should set ntpr=1 to see the details of what is happening. Also, the
"xmin" (ntmin=3) will generally give a smoother and faster minimization
path.
...dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Sun Jan 04 2009 - 01:08:43 PST
