From: Hadi (dinpajooh_at_gmail.com)
Date: Wed Aug 06 2014 - 13:57:26 CDT
This thread started with a question about rigid models which use bond
constraints in MD simulations, i.e., when the structure can be defined by
bond constraints the network of bonds. This is one way to enforce a rigid
body in a simulation, i.e. by imposing a sufficient number of bond
constraints between the atoms in the unit. Therefore, in this thread I am
using the term "rigid model" to pertain to these methods.
On Tue, Aug 5, 2014 at 12:59 PM, Axel Kohlmeyer <akohlmey_at_gmail.com> wrote:
> On Tue, Aug 5, 2014 at 3:49 PM, Hadi <dinpajooh_at_gmail.com> wrote:
> > On Tue, Aug 5, 2014 at 12:29 PM, Kenno Vanommeslaeghe
> > <kvanomme_at_rx.umaryland.edu> wrote:
> >> Two questions:
> >> (1) knowing that computers don't have infinite precision, there will
> >> *always* be small deviations. The question to ask is: are they
> >> (2) Did you plot the same distribution for an MC simulation?
> > Yes. My MC code does produce a distribution which is a delta function.
> As I
> > mentioned, the NAMD results have slight deviations from a delta function.
> > However, I quite agree with you: in addition to the previous factors
> > (algorithm, rigid tolerance, time step), precision should be considered
> > well.
> you have to be careful with generalizations in your statements. NAMD
> is *one* MD code and the accuracy of SHAKE depends very much on the
> length of the time step and convergence parameter. also, there are
> *plenty other* MD codes that do offer to integrate rigid particles as
> rigid bodies via separately integrating translation and rotation and
> for those there is *not any* difference in the way how positions are
> updated to what you do. as a person dealing with statistics and
> distributions, you should know that a sample size of 1 doesn't allow
> for much generalization, right?
> Dr. Axel Kohlmeyer akohlmey_at_gmail.com http://goo.gl/1wk0
> College of Science & Technology, Temple University, Philadelphia PA, USA
> International Centre for Theoretical Physics, Trieste. Italy.
This archive was generated by hypermail 2.1.6 : Thu Dec 31 2015 - 23:21:05 CST