Re: AMBER: Sander / Anal

From: David A. Case <case.scripps.edu>
Date: Mon, 1 Dec 2003 12:58:32 -0800

On Mon, Dec 01, 2003, Emmanuel Giudice wrote:
>
> I'm comparing the energy output of Sander and anal. To do so I have
> performed with sander a "single point" energy calculation in vacuo on a
> very simple system...

>
> Here is my input file for Sander
>
> &cntrlHere is my input file for Sander
> imin=0, ntx=1, ntpr=1, ntwr=1, ntwx=1, ntwe=1,
> nscm=5000,
> ntf=2, ntc=2,
> ntb=0, igb=2
> nstlim=1, t=0.0, dt=0.002,
> cut=999.0,
> irest=0
> &end
> -----------------------------------------------------------------------------
>

All components of your energy are different, so my guess is that you are not
extracting the intial (step 0) results from sander. It's easier to set
imin=1 and maxcyc=1 to ensure you are using the same coordinates.

Also, you should turn off shake (set ntf=ntc=1), since anal also doesn't know
about shake.

...hope this helps...dac

-- 
==================================================================
David A. Case                     |  e-mail:      case.scripps.edu
Dept. of Molecular Biology, TPC15 |  fax:          +1-858-784-8896
The Scripps Research Institute    |  phone:        +1-858-784-9768
10550 N. Torrey Pines Rd.         |  home page:                   
La Jolla CA 92037  USA            |    http://www.scripps.edu/case
==================================================================
-----------------------------------------------------------------------
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:03 PST
Custom Search