[gmx-users] minimization

Nuno R. L. Ferreira nunolf at ci.uc.pt
Mon Apr 5 13:17:01 CEST 2004


Your potential energy seems good, so my guess is that minimization went just fine.

That error you are receiving iis just an warnig. No problem at all.
Regards,
Nuno



######################################
Nuno Ricardo Santos Loureiro da Silva Ferreira
Departamento de Química
Faculdade de Ciências e Tecnologia
Universidade de Coimbra
3004-535 Coimbra - Portugal
Fax: +351 239 827703  - www.biolchem.qui.uc.pt
########################################
  ----- Original Message ----- 
  From: gold gunaseelan 
  To: spoel at xray.bmc.uu.se ; spoel at gromacs.org 
  Cc: gmx-users at gromacs.org 
  Sent: Monday, April 05, 2004 12:11 PM
  Subject: [gmx-users] minimization


  dear sir ,
  Iam using gromacs 3.2.1
  I changed the following in em.mdp  file as

  define = -DFLEXIBLE and ran energy minimisation but it converged as 

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

  Double precision normally gives you higher accuracy.

  writing lowest energy coordinates.

  Steepest Descents converged to machine precision in 202 steps,
  but did not reach the requested Fmax < 2000.
  Potential Energy  = -4.0429609e+05
  Maximum force    =  4.6900386e+03 on atom 1650
  Norm of force    =  1.4885711e+04

  Previously i had given define =-DFLEX_SPC it gave the same error but it converged in 57 steps.I looked in one of ur old reply,in which u have asked to see toplogy file 

  > ; Include water topology
  > #ifdef FLEX_SPC
  > #include "flexspc.itp"
  > #else
  > #include "spc.itp"
  > #endif
  > 
  > then you can try with this file:
  > 
  > 
  > title                    = xxxx
  > cpp                      = /lib/cpp
  > include                  = -I../top
  > define                  = -DFLEX_SPC

  I would like to point out that tis has been changed a bit in 3.2
  If you use spc water this still will work, but for other models the more
  generic way of writing is:

  #include "water.itp" (either tip3p, tip4p, spc, spce)

  define = -DFLEXIBLE

  But i didnot find this in top file

  ; Include water topology
  #include "spc.itp"

  #ifdef POSRES_WATER
  I found these terms missing as specified by u above.
  sir please suggest whats wrong .is protein minimized properly.thanks in advance -gsmith 



   
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20040405/41314dcd/attachment.html>


More information about the gromacs.org_gmx-users mailing list