Re: NAMD Energies

From: Peter Freddolino (petefred_at_umich.edu)
Date: Wed Dec 14 2016 - 13:54:52 CST

Quick followup on Brian’s point — the -debug flag or debug checkbox are what you want to keep the config file so that you can look and see if its settings match yours.
Best,
Peter

> On Dec 14, 2016, at 2:42 PM, Radak, Brian K <bradak_at_anl.gov> wrote:
>
> The log output is the energies used in the simulation - all terms and cutoffs. PME counts as ELECT and LJcorrections count as VDW.
>
> The energy plugin might make certain fixed assumptions on the cutoffs, etc. that don't match your simulation settings. If memory serves, you can save the NAMD input file that the plugin uses and compare the settings. I can't think of any particular reason why you shouldn't be able to match to at least 2-3 decimal places with the same settings.
>
> I have no idea what the conformational energy output does - I'd expect this is documented with the plugin? You might also be able to infer from the NAMD input file that is used if you can get your hands on it.
>
> Cheers,
> Brian
>
> Brian Radak
> Postdoctoral Appointee
> Leadership Computing Facility
> Argonne National Laboratory
>
> 9700 South Cass Avenue, Bldg. 240
> Argonne, IL 60439-4854
> (630) 252-8643
> brian.radak_at_anl.gov
>
> ________________________________________
> From: owner-namd-l_at_ks.uiuc.edu [owner-namd-l_at_ks.uiuc.edu] on behalf of Chris Goedde [chris.goedde_at_gmail.com]
> Sent: Wednesday, December 14, 2016 1:13 PM
> To: namd-l_at_ks.uiuc.edu
> Subject: namd-l: NAMD Energies
>
> Hi all,
>
> Couple of questions.
>
> How are the energies reported in the NAMD log file computed? The bonded energies are pretty obvious; I’m wondering more about the electrostatics and VdW. In particular, are the values for cutoff, pairlistdist, etc. used, or is the full energy of all interactions calculated? If PME is on, does that affect the electrostatic energy calculation?
>
> Second question. When I use the NAMD Energy plugin in VMD on a dcd file, the electrostatic energy is a bit different from that reported in the namd log file (about 10% different). What could be the cause of this? Also, what goes into the conformational energy reported by this plugin? It doesn’t seem to be a sum of other energies, that I can tell …
>
> Thanks.
>
> Chris
>
>
>

This archive was generated by hypermail 2.1.6 : Tue Dec 27 2016 - 23:22:44 CST