Minimizing a structure using OPTIM and AMBER9

From Docswiki
Revision as of 15:07, 14 March 2011 by import>Jmc49
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

To optimize a geometry you need four (or five, see below) files with names specified in the following list:

1. coords.inpcrd -- file with coordinates in AMBER9 format

2. coords.prmtop -- AMBER9 topology file (see description here)

3. min.in - file with AMBER9 keywords describing the conditions for AMBER minimization of the studied system (only one cycle!), e.g.:

  Minimization
  &cntrl
     imin=1, maxcyc=1, ncyc=1,
     igb=2, saltcon=0.2,
     ntx=1, ntb=0,
     cut=15.0, rgbmax=12.0, ifswitch=1
   /

4. odata -- file with OPTIM keywords specifying method of minimization, convergence criterion, etc. e.g.:

  BFGSMIN 1.0D-6
  BFGSSTEPS 10000
  UPDATES 1000
  MAXERISE 1.0D-4
  AMBER9 coords.inpcrd inpcrd

If the file for the starting geometry contains three columns with x, y, z coordinates, respectively, name this file start and change the last line of mentioned above odata file into:

  AMBER9 start

You still need to have a coords.inpcrd file in the same directory or OPTIM won't run!