AMBER: GBSA problem

From: <northrop.chem.ucla.edu>
Date: Thu, 11 Dec 2003 23:28:27 -0000

AMBER users,

I'm trying to run a simulation of a benzamide guest in a beta-cyclodextrin host with the implicit
GBSA water model. I've been using the gaff force field (with bugfix 7 "GBSA calculations using
the GAFF force field give a 'bad atom type' error" applied). In the input file I'm using igb = 1 and
gbsa = 1. When I try to do an initial minimization or if I try md, I get only question marks for all
my energy terms. The input files I'm using for minimization and the output file I get are below:

minimization of guest1, implicit water solvent
 &cntrl
  imin = 1, irest = 0, iwrap = 1, ntpr = 10, ntx = 1, igb = 1, gbsa = 1,

  maxcyc = 100, ncyc = 500, ntmin = 1, dx0 = 0.1, drms = 0.0001,

  ntb = 0, ntf = 1, cut = 10.0, nsnb = 10, ntc = 1,
 &end

representative output:
--------------------------------------------------------------------------------
   4. RESULTS
--------------------------------------------------------------------------------



   NSTEP ENERGY RMS GMAX NAME NUMBER
      1 6.1447E+01 3.7538E+01 1.9301E+02 S16 28

 BOND = 83.9909 ANGLE = 4.5039 DIHED = 0.0235
 VDWAALS = 1.6710 EEL = 19.2856 EGB = -16.0109
 1-4 VDW = 6.4571 1-4 EEL = -40.4250 RESTRAINT = 0.0000
 ESURF = 1.9505


   NSTEP ENERGY RMS GMAX NAME NUMBER
     10 ????????????? ????????????? 0.0000E+00 C2 1

 BOND = ????????????? ANGLE = ????????????? DIHED = ?????????????
 VDWAALS = ????????????? EEL = ????????????? EGB = ?????????????
 1-4 VDW = ????????????? 1-4 EEL = ????????????? RESTRAINT = 0.0000
 ESURF = 0.0000


   NSTEP ENERGY RMS GMAX NAME NUMBER
     20 ????????????? ????????????? 0.0000E+00 C2 1

 BOND = ????????????? ANGLE = ????????????? DIHED = ?????????????
 VDWAALS = ????????????? EEL = ????????????? EGB = ?????????????
 1-4 VDW = ????????????? 1-4 EEL = ????????????? RESTRAINT = 0.0000
 ESURF = 0.0000
.......

Has anyone come across this problem before? I've also tried to run the simulation using the 99
force field and have the same problem. Any suggestions would be greatly appreciated.

Thanks in advance,
 - Brian Northrop

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jan 14 2004 - 15:53:08 PST
Custom Search