Re: Minimisation not converged

From: Kevin C Chan (
Date: Wed Apr 01 2015 - 22:26:19 CDT

Thanks Jerome,

Problem was solved by running more steps of sd. I have also loaded the
trajectory and it was flipping of a flexible loop.

However in fact with sd minimisation, the energy keeps decreasing instead
of converges. I understand sd is a relatively "slow" minimisation but how
can I tell when to stop? Now I have run more than 100,000 steps and I
choose to step as the decrease of energy between steps was <0.01 kcal/mol.


Sent from my iPhone

On 1 Apr, 2015, at 15:14, Jérôme Hénin <
<javascript:_e(%7B%7D,'cvml','');>> wrote:

Hi Kevin,

Two things you can try:
1) open that trajectory in VMD and see exactly what is oscillating, just in
case it indicates a deeper problem;
2) run a few steps of steepest descent.


On 1 April 2015 at 08:20, Kevin C Chan <
<javascript:_e(%7B%7D,'cvml','');>> wrote:

> Dear Users,
> I am currently doing minimisation for a protein in water simply in a
> stepwise manner. For the first minimisation in which I fixed the protein
> and only water and ion could move, it worked fine. However for the second
> one in which I fixed water and backbone of the protein instead and let the
> side chains to relax, when I look at the energies using NamdPlot (
>, it converges
> to two distinct energy values and jumps between them every few steps.
> I have tried adjusting parameters like minBabyStep or minLineGoal. It does
> not alter the fact of bouncing energy values but changes the higher energy
> value, i.e. it now bounces from -1.0087e6 to -1.00865e6 instead of
> -1.0084e6.
> I am wondering whether I am having problem with the protein structure or
> minimisation parameters.
> Thanks in advance,
> Kevin
> PhD Candidate
> Department of Physics and Material Science
> City University of Hong Kong
> <javascript:_e(%7B%7D,'cvml','');>

This archive was generated by hypermail 2.1.6 : Thu Dec 31 2015 - 23:21:47 CST