[gmx-users] minimization force criteria

tarak karmakar tarak20489 at gmail.com
Thu Jul 26 14:30:42 CEST 2012


Dear All,

While running minimization I imposed the the condition for the
minimization as to be converged only at Fmax < 10 . But I got the
following


Stepsize too small, or no change in energy.
Converged to machine precision,
but not to the requested precision Fmax < 10

Double precision normally gives you higher accuracy.
You might need to increase your constraint accuracy, or turn
off constraints alltogether (set constraints = none in mdp file)

Steepest Descents converged to machine precision in 15318 steps,
but did not reach the requested Fmax < 10.
Potential Energy  = -3.3397742e+06
Maximum force     =  9.8232233e+02 on atom 6057
Norm of force     =  3.6321018e+00


Then according to the above instruction I increased the stepsize up to
0.5 [ I didn't use any constraint during minimization] but got the
same. So do I need to increase the stepsize more or is there other
procedure to get to the desired convergence criteria ?
The minimization .mdp file is as follows


; 7.3.2 Preprocessing
;define                  = -DFLEXIBLE    ; defines to pass to the preprocessor

; 7.3.3 Run Control
integrator              = steep         ; steepest descents energy minimization
nsteps                  = 100000         ; maximum number of steps to integrate

; 7.3.5 Energy Minimization
emtol                   = 10            ; [kJ/mol/nm] minimization is
converged when max force is < emtol
emstep                  = 0.01           ; [nm] initial step-size
nstcgsteep              = 10            ;  freq of performing one
steepest descent step

; 7.3.8 Output Control
nstxout                 = 500           ; [steps] freq to write
coordinates to trajectory
nstvout                 = 500           ; [steps] freq to write
velocities to trajectory
nstfout                 = 500           ; [steps] freq to write forces
to trajectory
nstlog                  = 100           ; [steps] freq to write
energies to log file
nstenergy               = 100           ; [steps] freq to write
energies to energy file
energygrps              = System        ; group(s) to write to energy file

; 7.3.9 Neighbor Searching
nstlist                 = 1             ; [steps] freq to update neighbor list
ns_type                 = grid          ; method of updating neighbor list
pbc                     = xyz           ; periodic boundary conditions
in all directions
rlist                   = 0.8           ; [nm] cut-off distance for
the short-range neighbor list

; 7.3.10 Electrostatics
coulombtype             = PME           ; Particle-Mesh Ewald electrostatics
rcoulomb                = 0.8           ; [nm] distance for Coulomb cut-off

; 7.3.11 VdW
vdwtype                 = cut-off       ; twin-range cut-off with
rlist where rvdw >= rlist
rvdw                    = 0.8           ; [nm] distance for LJ cut-off
DispCorr                = Ener          ; apply long range dispersion
corrections for energy

; 7.3.13 Ewald
fourierspacing          = 0.12          ; [nm] grid spacing for FFT
grid when using PME
pme_order               = 4             ; interpolation order for PME, 4 = cubic
ewald_rtol              = 1e-5          ; relative strength of
Ewald-shifted potential at rcoulomb


Thanks,

Tarak



More information about the gromacs.org_gmx-users mailing list