AMBER Archive (2004)Subject: Re: AMBER: Amber 7 (sander) use of PME
From: Carlos Simmerling (carlos_at_csb.sunysb.edu)
Date: Thu May 27 2004 - 14:47:50 CDT
I think that what you're seeing is reasonable and expected.
with use_pme=0, you basically have a cutoff on long-range
electrostatics and thus there are large fluctuations as
atoms move in and out of the cutoff. this "one factor" is
a really important one...
----- Original Message -----
From: "Linda Prengaman" <lindap_at_cmu.edu>
To: <amber_at_scripps.edu>
Sent: Thursday, May 27, 2004 3:42 PM
Subject: AMBER: Amber 7 (sander) use of PME
> Hello,
> I am running a system consisting of one amino acid solvated in a small
> water box (~300 TIP3 waters), cut=8.0, box size is~ 20x20x20. I have
> found that running MD simulations with use_pme=1 and use_pme=0 give me
> very different results for the total system energy: if pme is turned
> on, my total system energy is very nice and stable in a minimum energy
> state, but if pme is turned off, this energy appears to fluctuate
> wildly. I have been told that turning pme on is the "way to go"
> because it's more accurate, but I'm concerned about the accuracy of my
> simulation because of the differences between my results that come
> from changing just this one factor. Can anyone please give me some
> advice or explaination about using PME so I can resolve this? Thanks!
>
> Linda Prengaman
> Dept. of Chemistry
> Carnegie Mellon University
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|