[gmx-users] system not equilibrated
    Gmx Niki 
    gmx_gmx_12 at yahoo.com
       
    Wed Jan  9 09:41:09 CET 2013
    
    
  
Hi all,
I ran energy minimization (em) step, (protein-CNT in water) by these conditions:
title        = Minimization     
integrator    = steep         
emtol        = 1000.0       
emstep      = 0.01       
nsteps        = 150000  
energygrps    = Protein CNT     
 nstlist        = 1             
ns_type        = grid      
rlist        = 1.0         
coulombtype    = PME         
rcoulomb    = 1.0         
rvdw        = 1.0         
pbc            = xyz      
but when it recieves to 40000 step (nearly), it gives error:
Step=  125, Dmax= 7.8e-06 nm, Epot= -1.79432e+07 Fmax= 8.35837e+04, atom= 10656
Step=  129, Dmax= 1.2e-06 nm, Epot= -1.79432e+07 Fmax= 8.35979e+04, atom= 10656
Stepsize too small, or no change in energy.
Converged to machine precision,
but not to the requested precision Fmax < 1000
Double precision normally gives you higher accuracy.
writing lowest energy coordinates.
Back Off! I just backed up em2.gro to ./#em2.gro.1#
Steepest Descents converged to machine precision in 130 steps,
but did not reach the requested Fmax < 1000.
Potential Energy  = -1.7943236e+07
Maximum force     =  8.3583656e+04 on atom 10656
Norm of force     =  2.1333945e+03
then I ran my system with emtol=10 , but the error repeated too. then I ran it by double-precision , but was not successful too.
1- I changed fmax value only in .mdp file, is it correct? or I should change any parameter in topology file?
2- is the nsteps too much ?and problem is long time?what is its problem?
please guide me
tnx
Niki
    
    
More information about the gromacs.org_gmx-users
mailing list