[gmx-users] TIP5P and cg minimization

[rob yang]

Hi list,
I have solvated my protein in TIP5P box that's provided by 3.3.3. However, when I try to run an energy minimization using cg (at the grompp stage), I get:
ERROR: can not do Conjugate Gradients with constraints (218328)
 
I don't have any problem using the same em.mdp file on TIP4P, so I figured it has something to do with TIP5P. Steep descent works, but it bugs me that cg wouldn't work. As far as I can tell, there's no constraints in my em.mdp file and my *.top (*.itp files too). Has anybody else encounter the same problem? Thank you in advance.
 
em.mdp:
cpp                 =  /lib/cppdefine              =  -DFLEXIBLEconstraints         = noneunconstrained_start = yesintegrator          =  cgdt                  = 0.002nsteps              =  1000nstlist             =  10ns_type             =  gridrlist               =  1.0rcoulomb            =  1.0rvdw                =  1.0;;       Energy minimizing stuff;emtol               =  100.0emstep              =  0.01
-Rob
_________________________________________________________________

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