[gmx-users] simulated annealing

Justin A. Lemkul jalemkul at vt.edu
Wed Apr 21 13:33:32 CEST 2010

leila karami wrote:
> Dear Justin
> thanks for your attention
> I want to use simulated annealing to obtain global minimum for 
> protein-dna structure. I will do first heating and then slowly cooling.
> is my mdp file true, now?
> T-couple = brendsen
> ref-t = 300
> annealing: single
> annealing_npoints: 5
> annealing_time: 0  3  6  9  12 
> annealing_temp: 150  650  600  550  500  450  400  350  300

No.  Using ':' to separate lines will cause grompp to not read the parameters. 
Use '=' to separate parameters and values.  Further, the parameter 'T-couple' 
does not exist, it is tcoupl.  Also note that 'Berendsen' is spelled 
incorrectly, which will also cause grompp to give an error about an unrecognized 

I'm also assuming (based on having only a single value of annealing_npoints) 
that you are coupling the entire system to one thermostat.  Doing so is 
unreliable, and the consequences have been discussed numerous times across the 
list and in the references provided on the Gromacs wiki:


Further, you give 9 values of annealing_temp, but based on annealing_npoints, 
you should have 5 (if you only have this single group).  Please refer to the 
manual.  An example of proper usage is even provided online:



> Best regards.


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080


More information about the gromacs.org_gmx-users mailing list