From: Axel Kohlmeyer (akohlmey_at_gmail.com)
Date: Mon Dec 05 2011 - 17:36:58 CST
On Mon, Dec 5, 2011 at 3:54 PM, M.Raviprasad <M.Raviprasad_at_my.ndsu.edu> wrote:
> To NAMD users,
> I calculated distance traveled by an Oxygen atom (gas molecule) in the
> polymeric system using DCD file. In order to make sure the calculation, I
> used velocity to calculate the distance (distance = Velocity*time step).
> Velocity was taken from the velocity DCD file. However, these two values are
> not matching. I believe that results should be the same (not exact but
> close) irrespective of the methods.
i disagree. if it was possible to even approximately predict
the trajectory of atoms this way, we wouldn't have these
issues with requiring so much computer time to do MD
simulations of big systems.
> What am I doing wrong?
you are completely ignoring the fact how time integration works
and the fact that atomic velocities are significantly impacted
by high frequency vibrations and librations, that would integrate
out on long time scales.
> Is there any other way to calculate and verify the distance travelled by an
i don't see why you need to "verify" the coordinate data.
the coordinate dcd file is the perfect(!) record of how each
atom has traveled during your simulation and the distance
between any two reference points is as accurate as it gets
(outside of quantum effects that don't show up in classical
MD in the first place).
> Any suggestion is appreciated
-- Dr. Axel Kohlmeyer akohlmey_at_gmail.com http://goo.gl/1wk0 College of Science and Technology Temple University, Philadelphia PA, USA.
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:57:59 CST