From: Mert Gür (gurmert_at_gmail.com)
Date: Fri Oct 05 2012 - 15:01:03 CDT
Thanks a lot for your response.
The potential in TMD is
U= 0.5 k/N ( RMS(t)-RMS*(t) )^2
If you take the gradient of it then the force becomes (According to my
F= -k/N^2 ( 1-RMS*(t)/RMS(t) ) (R-R_target)
So there is an extra N. Does " k = K/it->second.size()" takes care of that ?
I use the following definition for RMS. Maybe the person who implemented it
used a different definition.
I am sending you a separate e-mail with my detailed derivation.
On Fri, Oct 5, 2012 at 2:26 PM, JC Gumbart <gumbart_at_ks.uiuc.edu> wrote:
> What about this line?
> k = K/it->second.size();
> On Oct 2, 2012, at 12:43 PM, Mert Gür wrote:
> Dear all,
> I wonder if the equation shown in the following link is correct?
> Is it k or k/N. If it is k/N (as written in the link) then there must be a
> second division by N in the source code (
> which I can not find.
> Can anyone comment on this?
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:22:09 CST