Dear Dr. Case,
Thank you very much for your help. It gives me the same results if I set
"vdwmeth = 0". I will check the TraPPE paper for the right option of these
flags.
All the best.
Lu
On 5/24/06, David A. Case <case.scripps.edu> wrote:
>
> On Wed, May 24, 2006, luckyang.gmail.com wrote:
> >
> > I have identified the source of the difference. The main difference
> comes
> > from the EVDW term. How does use_pme affect VDW energy? In attachment
> please
> > find the file contains the differences.
>
> My best guess is the the longrange vdw correction term is a part of the
> "PME"
> code, so that you are getting this in one calculation but not the other.
> Try setting vdwmeth=0 in both runs to see if you get the same results.
>
> I don't know if your force field was parameterized assuming such a long
> range
> correction or not. Leaving it out would lead to a slight change in the
> density of the system, but it would depend on the force field to know
> which
> one is "better".
>
> ...dac
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri May 26 2006 - 21:05:19 PDT
