AMBER Archive (2005)

Subject: AMBER: igb10 vlimit exceeded

From: Wells, David H (d.h.wells_at_imperial.ac.uk)
Date: Thu Mar 31 2005 - 03:08:55 CST


  
I'm trying to simulate a small organic solute in DMSO (83 molecules)
with igb=10. Previously there was trouble with solvent "penetrating" the
solute, apparently due to the way the non-bonded list is generated in
igb=10. That can be gotten around, but now simulations more than 20 ps
are blowing up. The RESTRAINT energy ramps rapidly, vlimits start to be
breached, and eventually temperature rises rapidly and solvent molecules
are forcibly ejected from the "cap" (even with FCAP=1000.0). Changing
igb=10 to igb=0 eliminates the problem.
Are the FCAP restraints only applied on atoms with indexes > NATCAP ?
Any further insights would be appreciated.
Many thanks, Dave
 



image001.gif
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu